Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	309,549	A→C	K13Q (AAA→CAA)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	309,549	0	A	C	100.0%	83.0 / NA	30	K13Q (AAA→CAA)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (9/21);  total (9/21)

ATTGATAAGCGACTTCTGCTTCTGGGATCCCGTCGGTTTTGTGCGGGTCAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGAAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATGTGTTTGAGCTCCCTTGGCTTTCCAGTAATGCAAGAAAATTACGTATTTCATCATGTAACCCGTTTTACCAATAATCGGCGCTATGTTATAGCAATC  >  1/309304‑309786                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              atTGATAAGCGACTTCTGCTTCTGGGATCCCGTCGGTTTTGTGCGGGTCAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAAt                                                                                                                                                                                                                                          <  1:124583/251‑1 (MQ=255)             cttctgcttctgGGATCCCGTCGGTTTTGTGCGGGTCAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTa                                                                                                                                                                                                                              >  1:773813/1‑251 (MQ=255)              ttctgcttctgGGATCCCGTCGGTTTTGTGCGGGTCAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTaa                                                                                                                                                                                                                             <  1:218730/251‑1 (MQ=255)              ttctgcttctgGGATCCCGTCGGTTTTGTGCGGGTCAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTaa                                                                                                                                                                                                                             <  1:148448/251‑1 (MQ=255)                         ggATCCCGTCGGTTTTGTGCGGGTCAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTc                                                                                                                                                                                                                   <  1:185389/250‑1 (MQ=255)                             cccGTCGGTTTTGTGCGGGTCAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTAttt                                                                                                                                                                                                                                  >  1:554249/1‑231 (MQ=255)                                            cGGGTCAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGatttatt                                                                                                                                                                                               <  1:533354/251‑1 (MQ=255)                                                 cAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTActc                                                                                                                                                                                          <  1:15242/251‑1 (MQ=255)                                                             gCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTAc                                                                                                                                                                                            <  1:604314/237‑1 (MQ=255)                                                                           ttgATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAgg                                                                                                                                                                <  1:183143/251‑1 (MQ=255)                                                                                 aCGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGaaaa                                                                                                                                                          >  1:764349/1‑251 (MQ=255)                                                                                    atatTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATaaaaa                                                                                                                                                                               >  1:285668/1‑227 (MQ=255)                                                                                      atTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGaa                                                                                                                                                                                                        <  1:873607/200‑1 (MQ=255)                                                                                          tgtgTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGttt                                                                                                                                                 <  1:883859/251‑1 (MQ=255)                                                                                                                                 gcATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTc                                                                                                            >  1:199767/1‑249 (MQ=255)                                                                                                                                       aCGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTaa                                                                                                    <  1:761703/251‑1 (MQ=255)                                                                                                                                            gCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATaa                                                                                                                       <  1:99316/227‑1 (MQ=255)                                                                                                                                               aGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGt                                                                                                                                                   <  1:223930/196‑1 (MQ=255)                                                                                                                                                             cAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTAttt                                                                                                                                                                                                                                  <  1:367133/103‑1 (MQ=255)                                                                                                                                                              aaTGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATGTGTTTGAGCTCCCTTGGCttt                                                                             <  1:155464/251‑1 (MQ=255)                                                                                                                                                                 ggTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATaa                                                                                                                                                                                  >  1:742647/1‑147 (MQ=255)                                                                                                                                                                           aGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTaaaaa                                                                                                                                       <  1:527666/180‑1 (MQ=255)                                                                                                                                                                                             tGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATGTGTTTGAGCTCCCTTGGCTTTCCAGTaa                                                                      >  1:36187/1‑227 (MQ=255)                                                                                                                                                                                                                  aTGAATTATCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATGTGTTTGAGCTCCCTTGGCTTTCCAGTAATGCAAGAAAATTACGTATTTCATCATGTAACCCGTTTTACCaata                         <  1:71375/251‑1 (MQ=255)                                                                                                                                                                                                                          tCGTAAGAGTCAGAAAAATGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATATGTTTGAGCTCCCTTGGCTTTCCAGTAATGCAAGAAAATTACGTATTTCATCATGTAACCCGTTTTACCAATAATCGgcg                  >  1:451548/1‑250 (MQ=255)                                                                                                                                                                                                                                       aaaaaTGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATGTGTTTGAGCTCCCTTGGCTTTCCAGTAATGCAAGAAAATTACGTATTTCATCATGTAACCCGTTTTACCAATAATCGGCGCTATGTTATAg       <  1:47836/248‑1 (MQ=255)                                                                                                                                                                                                                                       aaaaaTGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATGTGTTTGAGCTCCCTTGGCTTTCCAGTAATGCAAGAAAATTACGTATTTCATCATGTAACCCGTTTTACCAATAATCGGCGCTATGTTATAGc      <  1:70549/249‑1 (MQ=255)                                                                                                                                                                                                                                       aaaaaTGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATGTGTTTGAGCTCCCTTGGCTTTCCAGTAATGCAAGAAAATTACGTATTTCATCATGTAACCCGTTTTACCAATAATCGGCGCTATGTTATAGc      <  1:340151/249‑1 (MQ=255)                                                                                                                                                                                                                                         aaaTGATCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATGTGTTTGAGCTCCCTTGGCTTTCCAGTAATGCAAGAAAATTACGTATTTCATCATGTAACCCGTTTTACCAATAATCGGCGCTATGTTATAGCAATc  <  1:106488/251‑1 (MQ=255)                                                                                                                                                                                                                                              aTCTATCGCAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCtg                                                                                                         >  1:285964/1‑143 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              ATTGATAAGCGACTTCTGCTTCTGGGATCCCGTCGGTTTTGTGCGGGTCAAAACTTTGTGGCTCGGAACCGTTGTTGATCACGATATTCTGTGTTTCTGCGAGCTGTGTTCCTTCAGGTACTTTGGCGGCATAAACGCTGCAAGATAAGGCAAGCGCAATGGTTGAGACAAGTAATTTTTGTTGCATTTGAATCTCCTTGTTATGAGCAATGAATTATCGTAAGAGTCAGAAAAATGATCTATCGAAAAATAGTTATTTTGTAAAGATCAATTCACAAAAAAGAACGATTTATTTACTCTTCTATAAAAAAGGAAAAAATAAAGGCGAAAAATTATGTTTTATTTAAAAATTTCGCCTCTAAATAAGGACCGATTTCCTGTGTAATGTGTTTGAGCTCCCTTGGCTTTCCAGTAATGCAAGAAAATTACGTATTTCATCATGTAACCCGTTTTACCAATAATCGGCGCTATGTTATAGCAATC  >  1/309304‑309786

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑