Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,031,075	T→C	K43K (AAA→AAG)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,031,075	0	T	C	100.0%	92.0 / NA	33	K43K (AAA→AAG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/15);  total (18/15)

ATATTATGCTGTTTATCATACTTTGAGTAACTGAAAGACTAAAACCTTTGGCTTCTAGTTTTGGGAAGTCAATGTTACTTATAATAGATGAAAATGCTTTCTTTATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTTTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTT  >  1/1030837‑1031317                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   atatTATGCTGTTTATCATACTTTGAGTAACTGAAAGACTAAAACCTTTGGCTTCTAGTTTTGGGAAGTCAATGTTACTTATAATAGATGAAAATGCTTTCTTTATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGaaa                                                                                                                                                                                                                                        >  1:471045/1‑251 (MQ=255)             ttATCATACTTTGAGTAACTGAAAGACTAAAACCTTTGGCTTCTAGTTTTGGGAAGTCAATGTTACTTATAATAGATGAAAATGCTTTCTTTATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGaaa                                                                                                                                                                                                                            >  1:845713/1‑251 (MQ=255)                                         aaaaCCTTTGGCTTCTAGTTTTGGGAAGTCAATGTTACTTATAATAGATGAAAATGCTTTCTTTATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCGATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTg                                                                                                                                                                                                <  1:687239/251‑1 (MQ=255)                                              cTTTGGCTTCTAGTTTTGGGAAGTCAATGTTACTTATAATAGATGAAAATGCTTTCTTTATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGt                                                                                                                                                                                           <  1:416036/251‑1 (MQ=255)                                                                                  taataGATGAAAATGCTTTCTTTATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGttt                                                                                                                                                       >  1:158180/1‑251 (MQ=255)                                                                                                     ctttATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAAttt                                                                                                                                    >  1:558871/1‑251 (MQ=255)                                                                                                                    ataatCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCttgat                                                                                                                                                                  <  1:729162/206‑1 (MQ=255)                                                                                                                       atCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACACAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATaa                                                                                                                  >  1:545540/1‑251 (MQ=255)                                                                                                                           tttttCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCAcc                                                                                                                                                                                                  >  1:207358/1‑167 (MQ=255)                                                                                                                                               aGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCagt                                                                                                                             >  1:553330/1‑216 (MQ=255)                                                                                                                                                  gtTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATtctc                                                                                        >  1:737754/1‑250 (MQ=255)                                                                                                                                                      ataGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTc                                                                                   <  1:391973/251‑1 (MQ=255)                                                                                                                                                          cAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGa                                                                                <  1:169220/250‑1 (MQ=255)                                                                                                                                                                            aTTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCttgatt                                                                                                                                                                 >  1:488062/1‑151 (MQ=255)                                                                                                                                                                            aTTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATtctc                                                                                        <  1:162474/224‑1 (MQ=255)                                                                                                                                                                            aTTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTCCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCAt                                                             >  1:821184/1‑251 (MQ=255)                                                                                                                                                                                    gAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAAt                                                                                                                                                                                      <  1:761475/122‑1 (MQ=255)                                                                                                                                                                                      aTCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAAt                                                                                                   <  1:845244/203‑1 (MQ=255)                                                                                                                                                                                      aTCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGAt                                                                                                >  1:382603/1‑206 (MQ=255)                                                                                                                                                                                          aTACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATg                                                                                                  <  1:641552/200‑1 (MQ=255)                                                                                                                                                                                            aCCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTggaggg                                             >  1:527454/1‑251 (MQ=255)                                                                                                                                                                                                  tAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCaa                                       >  1:516325/1‑251 (MQ=255)                                                                                                                                                                                                       aaaCATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTtaata                                                                                                                                                   >  1:870926/1‑138 (MQ=255)                                                                                                                                                                                                          cATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGa                                 >  1:114073/1‑249 (MQ=255)                                                                                                                                                                                                          cATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGACCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGa                                 >  1:103011/1‑249 (MQ=255)                                                                                                                                                                                                            tAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCagt                                                                                                                             <  1:464768/155‑1 (MQ=255)                                                                                                                                                                                                            tAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATtctc                                                                                        <  1:828964/192‑1 (MQ=255)                                                                                                                                                                                                                  tACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGc                                                                                                                                                                                         <  1:105935/89‑1 (MQ=255)                                                                                                                                                                                                                                 aaaTTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCaatgt         <  1:311202/250‑1 (MQ=255)                                                                                                                                                                                                                                   attattATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATtctc                                                                                        <  1:677812/169‑1 (MQ=255)                                                                                                                                                                                                                                       ttatGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATg                                                                                                  >  1:879286/1‑155 (MQ=255)                                                                                                                                                                                                                                       ttatGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGttt  >  1:773902/1‑251 (MQ=255)                                                                                                                                                                                                                                       ttatGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGttt  <  1:184885/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   ATATTATGCTGTTTATCATACTTTGAGTAACTGAAAGACTAAAACCTTTGGCTTCTAGTTTTGGGAAGTCAATGTTACTTATAATAGATGAAAATGCTTTCTTTATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTTTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTT  >  1/1030837‑1031317

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: ‑