Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,214	G→A	S436S (AGC→AGT)	thiP ←	thiamine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,214	0	G	A	100.0%	179.4 / NA	53	S436S (AGC→AGT)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/20);  total (33/20)

ACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTC  >  NC_000913/72968‑73450                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             aCAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGt                                                                                                                                                                                                                                          <  1:355001/251‑1 (MQ=255)   agagTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCCATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGt                                                                                                                                                                                                                                          >  1:214758/1‑249 (MQ=255)           agcagAATTAACGCGGTGACCGCGCCGTCCTGGCTTCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCCGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCACTGCCCAGCGACTGACCCAACCTACTGTAGCGGGCGGt                                                                                                                                                                                                                                >  2:264140/1‑251 (MQ=255)           agcagAATTAACGCGGTGACCGCGCCGTCCTGGCTGCTATAGTAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCTTTACCGACCAACTCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTGCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGCAgcgggcagt                                                                                                                                                                                                                                <  2:315965/251‑4 (MQ=255)           agcagAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCCGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACAAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAgcgggcggg                                                                                                                                                                                                                                >  1:188369/1‑250 (MQ=255)             cagAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTg                                                                                                                                                                                                                               >  1:101522/1‑250 (MQ=255)             cagAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGAAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTg                                                                                                                                                                                                                               >  1:39881/1‑250 (MQ=255)                          gTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCgg                                                                                                                                                                                                                 >  1:185336/1‑251 (MQ=255)                           tGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCggg                                                                                                                                                                                                                >  1:295797/1‑251 (MQ=255)                                       ccTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAgtac                                                                                                                                                                                                     <  1:354627/250‑4 (MQ=255)                                        cTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACtt                                                                                                                                                                                                   >  1:331302/1‑251 (MQ=255)                                           gCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCa                                                                                                                                                                                                <  2:265623/251‑1 (MQ=255)                                            cTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTc                                                                                                                                                                                                 >  2:64040/1‑249 (MQ=255)                                                   aGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATa                                                                                                                                                                                         >  1:309100/1‑250 (MQ=255)                                                                 tGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAgcgg                                                                                                                                                                                                                                     <  1:77122/192‑1 (MQ=255)                                                                 tGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAgcgg                                                                                                                                                                                                                                     >  2:77122/1‑192 (MQ=255)                                                                  ggcagagATAAAACGGCAGGGTGCGGAAATCAGCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAAc                                                                                                                                                                         <  2:150956/248‑1 (MQ=255)                                                                   gTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCCTAAGGGATCGCCATTAAcc                                                                                                                                                                        >  2:103707/1‑250 (MQ=255)                                                                      gagaTAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCAt                                                                                                                                                                     <  1:111346/251‑1 (MQ=255)                                                                       agaTAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCAtt                                                                                                                                                                    >  1:341432/1‑251 (MQ=255)                                                                       agaTAAAACGGCAGGGTGCGGAAATCATCGTGACCGAACAACGCCACCACGCAAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTTATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGCATCGTCATTAACGCact                                                                                                                                                                    <  2:362141/251‑3 (MQ=255)                                                                            aaaCGGCAGGGTGCGGAAAGCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCGGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGa                                                                                                                                                               <  2:133993/251‑1 (MQ=255)                                                                             aaCGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATcgc                                                                                                                                                                                                                        <  1:248453/193‑1 (MQ=255)                                                                             aaCGGCAGGGGGCGGAAATCATCGTTACCGAACAACGCCCCCACGCCCAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCCGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATcgc                                                                                                                                                                                                                        >  2:248453/1‑193 (MQ=255)                                                                                ggCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCAAGCCAAAATCACCAATCGCCAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAAt                                                                                                                                                           >  1:166012/1‑251 (MQ=255)                                                                                 gCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATc                                                                                                                                                          >  2:277722/1‑251 (MQ=255)                                                                                   aGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATcac                                                                                                                                                        <  1:264140/251‑1 (MQ=255)                                                                                       tGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAACAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATg                                                                                                                                                    >  1:154546/1‑251 (MQ=255)                                                                                          ggAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGATCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCg                                                                                                                                                 >  2:256612/1‑251 (MQ=255)                                                                                          ggAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCg                                                                                                                                                 >  1:281759/1‑251 (MQ=255)                                                                                           gAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGt                                                                                                                                                <  2:331302/251‑1 (MQ=255)                                                                                                           aacaacGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGc                                                                                                                                  >  1:353100/1‑249 (MQ=255)                                                                                                            acaacGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAgg                                                                                                                               >  1:19180/1‑251 (MQ=255)                                                                                                            acaacCCCACCATCCCAGAATCACCAATCGACAACACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCAATTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAAAATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATtgcgaccgg                                                                                                                               >  2:359611/1‑244 (MQ=255)                                                                                                                 gCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGccaat                                                                                                                          >  2:221812/1‑248 (MQ=255)                                                                                                                            aaaaTCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGttgttg                                                                                                                >  1:345307/1‑250 (MQ=255)                                                                                                                             aaaTCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCTCTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCATGCCAATGGTGTTGTTGAg                                                                                                              <  2:326340/251‑1 (MQ=255)                                                                                                                                    cAATCGACAGAACGCAGGAAAAAGCCAGCGCCTGCGCCAGTGGACGTGTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTTACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGTTGGTGTTGAGCAGTaaa                                                                                                       <  2:166215/251‑1 (MQ=255)                                                                                                                                     aaTCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTaaaa                                                                                                      <  2:117182/251‑1 (MQ=255)                                                                                                                                        cGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGGTGTTGAGCAGTAAAaag                                                                                                    <  2:166012/250‑1 (MQ=255)                                                                                                                                         gACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAaaga                                                                                                   >  2:374877/1‑250 (MQ=255)                                                                                                                                                    ggCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCa                                                                                       <  2:321321/250‑1 (MQ=255)                                                                                                                                                          agCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACaa                                                                                 <  1:174333/250‑1 (MQ=255)                                                                                                                                                                             ggACGTTTCAGGGCGCGCAGCGCGACCACTGTTAAGCGCGACCAGCCTCCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCTCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTGGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCAGCATGGCGAGGAt                                                              <  2:270638/251‑1 (MQ=255)                                                                                                                                                                                cGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCa                                                            >  2:269170/1‑250 (MQ=255)                                                                                                                                                                                            cgcAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCAt                                               >  2:33571/1‑251 (MQ=255)                                                                                                                                                                                                   cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCa                                          <  2:308006/247‑1 (MQ=255)                                                                                                                                                                                                      ccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGcgcc                                     >  2:309716/1‑248 (MQ=255)                                                                                                                                                                                                         cTTTTAAGCGCGACCCGCCTTCAATGCCCAGCGACTGACACAAAATACTGTAGAGGGCGGTGATATCGCGCATCGGGTTTTCAAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATACCGTCAACAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACCATCACCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTg                                   >  1:182538/1‑250 (MQ=255)                                                                                                                                                                                                                               aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCcgcg            >  1:336400/5‑251 (MQ=255)                                                                                                                                                                                                                               aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGAGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCGCATCACCAGCGCCTGACCCGCCAGCCTTTTCTGCcgcg            >  1:356498/5‑251 (MQ=255)                                                                                                                                                                                                                                     cAGCGACTGAACCAACATACTGGAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGCCAATCACAATGCCGTCCGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGGAGCCAGCACACTCCCCGGCAAGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCCGACCCGCCAGCCTTTTCTGCCGCGCCCgac      >  2:240693/1‑249 (MQ=255)                                                                                                                                                                                                                                         gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTc  >  1:262059/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             ACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTC  >  NC_000913/72968‑73450

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