Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	72,736	T→G	T64T (ACA→ACC)	thiQ ←	thiamine ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	72,736	0	T	G	100.0%	183.1 / NA	51	T64T (ACA→ACC)	thiQ	thiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (31/20);  total (31/20)

GGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGAT  >  NC_000913/72443‑73017                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                          ggATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGa                                                                                                                                                                                                                                                                                     >  1:577694/1‑300 (MQ=255)   aTCGAGCGCAGAGAAGGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATc                                                                                                                                                                                                                                                                                   <  1:61547/300‑1 (MQ=255)    tCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGTCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATc                                                                                                                                                                                                                                                                                   <  2:18003/299‑1 (MQ=255)    tCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATca                                                                                                                                                                                                                                                                                  <  1:598937/300‑2 (MQ=255)          gcAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGCCAACGCGCTAACGCCACGCGCTGTCACTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCTGAGTGAGCCCTAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGc                                                                                                                                                                                                                                                                              <  2:257427/298‑1 (MQ=255)                    gTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATTTGGCGGGCGATAGCGTGCATTTTCTACTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAAGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAAGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTc                                                                                                                                                                                                                                                                   >  2:338287/1‑299 (MQ=255)                                aTAAAGTCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAAcagat                                                                                                                                                                                                                                                        >  2:135554/1‑294 (MQ=255)                                aTAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATGGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCAGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCgct                                                                                                                                                                                                                                                        >  2:136141/1‑298 (MQ=255)                                aTAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTACTGCTGTGTCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACAGGACGGCGTGACGGCGGCGTCGTGGTGTGGTCAACGCCGTCGATAGTCAGCGAACCgct                                                                                                                                                                                                                                                        >  2:136916/1‑298 (MQ=255)                                            gTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTTGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGa                                                                                                                                                                                                                                           <  2:283462/299‑1 (MQ=255)                                             ttCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGNTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGaaa                                                                                                                                                                                                                                         <  2:194574/300‑1 (MQ=255)                                                   caCCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGc                                                                                                                                                                                                                                   >  1:396623/2‑300 (MQ=255)                                                     ccAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGa                                                                                                                                                                                                                                 >  1:424554/1‑300 (MQ=255)                                                               cgcTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTca                                                                                                                                                                                                                        >  1:112901/1‑299 (MQ=255)                                                                  tAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCACCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTcagcag                                                                                                                                                                                                                    <  2:365077/300‑1 (MQ=255)                                                                         acgcgctgtcgctgACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACtt                                                                                                                                                                                                             <  1:182264/297‑1 (MQ=255)                                                                                gtcgctgACCGCCGGAAAGCTCGCCAGGTAACCGCGCAATTAAATTATCAATCCACATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCAGGATTTAGCACCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCAGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGCAAACCGGCGATCAAGTTCAGCAGGGTACTTTTAcct                                                                                                                                                                                                       >  2:325921/1‑298 (MQ=255)                                                                                                       ccaggTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGTGTGCATTTTCTACTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGcc                                                                                                                                                                               >  2:454705/4‑300 (MQ=255)                                                                                                          ggTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCAcc                                                                                                                                                                            <  1:158124/300‑1 (MQ=255)                                                                                                                                  aTCCCCATCTGGCGGGCGATGGCGTGCATTTTCTACTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCCCCTGCTCGCCGCGTTCTACCGTTaa                                                                                                                                                     >  2:548394/1‑299 (MQ=255)                                                                                                                                  aTCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTaa                                                                                                                                                     >  2:548387/1‑299 (MQ=255)                                                                                                                                           tGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAAc                                                                                                                                             <  1:609574/298‑1 (MQ=255)                                                                                                                                                 ggggtAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCAAGCCCGATGTTCTGCGCGACCGTCAGGTGGCTGAACAGGTTGTTTTCACGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCGGCGAACCGCTGGCTGGCGTCAGAACAACGGTGATCAAGTTCAGCAGGGTACCTTTCCCTGCGCCGATTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTGCAGCTAAGCTAAAACGCATGGGc                                                                                                                                     >  2:320260/5‑300 (MQ=255)                                                                                                                                                 gcgATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCAGGATTTAGCACCAGCCCGATGTTCTGCGCGACCGCCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCCGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCACGAGGATCGCCAACCGCTCGCCACGTTCTACCGTTAAGCTaaa                                                                                                                                               >  2:320479/1‑290 (MQ=255)                                                                                                                                                                  ctccTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCa                                                                                                                    >  2:330224/3‑300 (MQ=255)                                                                                                                                                                     ctgctgTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCAAGGt                                                                                                                 >  2:579180/1‑300 (MQ=255)                                                                                                                                                                           tgccGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCAATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAtc                                                                                                           <  2:317021/298‑1 (MQ=255)                                                                                                                                                                                 gTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCAAGGTGATAtcagtca                                                                                                      >  2:345892/1‑299 (MQ=255)                                                                                                                                                                                  ttCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATcagtcagt                                                                                                    >  1:452948/1‑300 (MQ=255)                                                                                                                                                                                        ttCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGNGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAAc                                                                                              >  2:195832/1‑300 (MQ=255)                                                                                                                                                                                         tCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACa                                                                                             <  1:602278/300‑1 (MQ=255)                                                                                                                                                                                                  gATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGccc                                                                                     >  1:452209/1‑299 (MQ=255)                                                                                                                                                                                                  gATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGNGTGATCAACGCCGTCGATAGTCAGCGAACCGCNGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGcccc                                                                                    >  2:482796/1‑300 (MQ=255)                                                                                                                                                                                                   aTTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTTGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGcccc                                                                                    >  1:601224/1‑299 (MQ=255)                                                                                                                                                                                                              gCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCa                                                                        >  2:531290/1‑300 (MQ=255)                                                                                                                                                                                                                ccGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACCTCTAGCCCAGGTAGTTTTTAAAt                                                                      >  2:190088/1‑300 (MQ=255)                                                                                                                                                                                                                                  gTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGTTAGTTTTTCAATCACGGTAAACAGCAGaaa                                                    <  1:63012/300‑1 (MQ=255)                                                                                                                                                                                                                                   tCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCACGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGaaaa                                                   >  2:358231/1‑300 (MQ=255)                                                                                                                                                                                                                                      ggTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACag                                                <  2:424554/300‑1 (MQ=255)                                                                                                                                                                                                                                        tgacTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTAACGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTAACCCCGGTCGTTTTTCAATCACGGTAAACAGCAGAAAACaga                                               >  2:530990/4‑299 (MQ=255)                                                                                                                                                                                                                                             gAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCTAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagtagca                                         >  1:401640/1‑300 (MQ=255)                                                                                                                                                                                                                                              aaCAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagtagcag                                        <  1:507627/300‑1 (MQ=255)                                                                                                                                                                                                                                                  ggTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGa                                    >  1:397072/1‑300 (MQ=255)                                                                                                                                                                                                                                                   gttgttTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGagtagcagc                                       <  1:266419/296‑1 (MQ=255)                                                                                                                                                                                                                                                         ttcCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAAcg                             <  1:491675/299‑1 (MQ=255)                                                                                                                                                                                                                                                                aaCAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGTTGATATCAGTCAGTTTTAACATCTCGCCCCGTTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGAc                      <  2:350867/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          acacCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTcc            <  2:577694/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          acacCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTAAGTTTTAACATCTCGCCCCGGTAGTTTTTAAATCACGGTAAACAGCAGAGAACAGAGTAGAAGCAGAAt                                  >  2:101025/1‑278 (MQ=255)                                                                                                                                                                                                                                                                               ggACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCt       <  1:568292/300‑1 (MQ=255)                                                                                                                                                                                                                                                                               ggACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCt       >  2:640887/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                    gCGTGACAGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTCACGCGGTGACCGCGCCGTCCTGGCTGcaat  >  2:306377/1‑297 (MQ=255)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                          GGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGAT  >  NC_000913/72443‑73017

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑