Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,412,466	G→A	T23I (ACT→ATT)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,412,466	0	G	A	100.0%	111.6 / NA	40	T23I (ACT→ATT)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/16);  total (24/16)

TAATAGATGAAAATGCTTTCTTTATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAGTAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTTCGCTCATGTAACTCGTAAACAATCCGCTTTCTACTCCAATCCTTTTTTTTCATAGTTCCTCG  >  1/1412248‑1412709                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    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)                                                              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)                                                                                                     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)                                                                                                                                                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)                                                                                                                                                                          aTTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTTCGCTCATGTAACTCGTAAAc                                            <  1:613601/251‑1 (MQ=255)                                                                                                                                                                              ttcCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACt                                                                                             <  1:613128/198‑1 (MQ=255)                                                                                                                                                                                   gAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCaatgtaatgt                                                                  >  1:626132/1‑220 (MQ=255)                                                                                                                                                                                        gCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACt                                                                                             <  1:221044/188‑1 (MQ=255)                                                                                                                                                                                              tACCAATTCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGttt                                                                <  1:102694/211‑1 (MQ=255)                                                                                                                                                                                                    ttCCCCATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTTCGCTCATGTAACt                                                   >  1:145754/1‑218 (MQ=255)                                                                                                                                                                                                        ccATTCACCTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTTCGCTCATGTAACTCGTAAACAATCCGCTTTCTACTCCAATCCtttttttt              >  1:38287/1‑251 (MQ=255)                                                                                                                                                                                                               ccTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTTCGCTCATGTAACTCGTAAACAATCCGCTTTCTACTCCAATCCTTTTTTTTCATAGtt       >  1:680530/1‑251 (MQ=255)                                                                                                                                                                                                               ccTGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTTCGCTCATGTAACTCGTAAACAATCCGCTTTCTACTCCAATCCTTTTTTTTCATAGtt       >  1:510180/1‑251 (MQ=255)                                                                                                                                                                                                                 tGTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTtaat                                                                                                                                                                                                                  >  1:872130/1‑46 (MQ=255)                                                                                                                                                                                                                  gTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAAttt                                                                                                                                                                                                             <  1:638159/50‑1 (MQ=255)                                                                                                                                                                                                                  gTAAGTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTTCGCTCATGTAACTCGTAAACAATCCGCTTTCTACTCCAATCCTTTTTTTTCATAGTTc      >  1:626668/1‑249 (MQ=255)                                                                                                                                                                                                                      gTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGttt                                                                >  1:680172/1‑187 (MQ=255)                                                                                                                                                                                                                      gTGCAATAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTTCGCTCATGTAACTCGTAAACAATCCGCTTTCTACTCCAATCCTTTTTTTTCATAGTTCCTCg  >  1:827800/1‑249 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    TAATAGATGAAAATGCTTTCTTTATATGTAATAAATAATCTTTTTTCTGCTTTGTATCATGAGTGTTTATAGCAATTTGATATTTAAATGATTCTAAGGAATCAATACCAAATAACGAAACATAAGAATACTTATTCATTTTGAAATTATTATGATTCTTTTCATAGAAATTCTTCCAGAAATGCGTTTTACCAATTCCCCATTCACCTGTAAGTGCAGTAGCAAAAGAGTCTGCTTGATTGATAAGGTTTAATAATTTTTCTTCAATTTTTTCAGTAGTTGCCATAATTGGTTACTCTAAATGATCCATTCTCCCTTCGGATAACTAACTCGCAATGCATTTTTTAATGTGGAGGGAGCCAATCCAGACTTAACTGAAAGCGAATGCAATGTAATGTTTCGCTCATGTAACTCGTAAACAATCCGCTTTCTACTCCAATCCTTTTTTTTCATAGTTCCTCG  >  1/1412248‑1412709

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