Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,001	A→G	G507G (GGT→GGC)	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,001	0	A	G	100.0%	139.5 / NA	42	G507G (GGT→GGC)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/22);  total (20/22)

GATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCG  >  NC_000913/72754‑73242                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  gATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCg                                                                                                                                                                                                                                                >  2:171722/1‑251 (MQ=255) gATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCg                                                                                                                                                                                                                                                >  2:204824/1‑251 (MQ=255)            aaCCGCTGGCTGGCGGCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGGTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCg                                                                                                                                                                                                                                     <  2:118615/251‑1 (MQ=255)                              gAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTg                                                                                                                                                                                                                   <  2:259123/251‑1 (MQ=255)                                           caagttCAGCAGGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAAc                                                                                                                                                                                                      <  1:231267/247‑1 (MQ=255)                                                 cagcagGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAgg                                                                                                                                                                                                >  2:106989/1‑251 (MQ=255)                                                     agGGTACTTTTACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGc                                                                                                                                                                                            <  1:115097/251‑1 (MQ=255)                                                              ttACCTGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCCAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGGAAACAGCAGAAAACAGAGGAGCAGCAGAATTAACACGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTAGCTGGTAGCGATAAAACGGCAGGGGGCGGAAatcat                                                                                                                                                                                   >  2:261897/1‑251 (MQ=255)                                                                    gcgcCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTAc                                                                                                                                                                             <  2:227780/251‑1 (MQ=255)                                                                     cgccgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTACAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTAcc                                                                                                                                                                            >  2:156624/1‑251 (MQ=255)                                                                             ggCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGACAGGGTGCGGAAATCATCGTTACCGaacaac                                                                                                                                                                     <  2:81413/250‑1 (MQ=255)                                                                               cccGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAAAAGAGTAGCAGCAGAATTAACGCGGTGAACGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGCAATCATCGTTACCGAACAACGcc                                                                                                                                                                  >  1:178708/1‑251 (MQ=255)                                                                                       tCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGcc                                                                                                                                                          >  1:280116/1‑251 (MQ=255)                                                                                          ccACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATAGCAGACAGTGTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCGGGCTGCGATAGGAGCCAAGTTGCTGGTAGAGATAAAACGGCAGGGGGCGGAAATCATCGTTTGCGAGCAACGCCACCACCCCaaa                                                                                                                                                       <  2:286097/251‑1 (MQ=255)                                                                                                   cgccgcGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCa                                                                                                                                                <  1:39920/249‑1 (MQ=255)                                                                                                   cgccgcGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCa                                                                                                                                                <  2:275202/249‑1 (MQ=255)                                                                                                    gccgcGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATc                                                                                                                                             <  2:187143/251‑1 (MQ=255)                                                                                                    gccgcGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAGTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATc                                                                                                                                             <  2:157731/251‑1 (MQ=255)                                                                                                                cGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGc                                                                                                                                   >  2:186609/1‑249 (MQ=255)                                                                                                                      gCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCaaa                                                                                                                            >  2:224638/1‑250 (MQ=255)                                                                                                                           aaCGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCac                                                                                                                      >  2:97801/1‑250 (MQ=255)                                                                                                                           aaCGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAg                                                                                                                      <  2:25141/251‑1 (MQ=255)                                                                                                                            aCGCATCGGCAAATGTTGGTAAAGCCAGGTGATAGCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAgc                                                                                                                     <  2:59185/251‑1 (MQ=255)                                                                                                                                           ttggtAAAGCCAGGTGATATAAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTgg                                                                                                       >  1:111770/2‑250 (MQ=255)                                                                                                                                              gtAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGt                                                                                                   <  2:257909/251‑1 (MQ=255)                                                                                                                                                     cAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAggg                                                                                            <  1:194837/251‑1 (MQ=255)                                                                                                                                                      aGGTGATCTCAGTCAGTTTTAACATCTCGACCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGgc                                                                                           >  2:18887/1‑251 (MQ=255)                                                                                                                                                         ggATATAAGTCAGTTTTAACATCCCGCCCCGGTAGTTTTTCAATCACGGAAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGTAGGGTGCGGAAATCATCGTTACCGAACAAAGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCAAGTGGACGTTTCAGGgcgcg                                                                                        <  2:273238/250‑1 (MQ=255)                                                                                                                                                          gATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAACCAGCAGAAAACAGAGTAGCAGCAGACTTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGGGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCACCAGTGGACGTTTCAGGgcgcgc                                                                                       >  2:129031/1‑251 (MQ=255)                                                                                                                                                             atCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGc                                                                                    <  1:156624/251‑1 (MQ=255)                                                                                                                                                               cagtcagtTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTc                                                                                  <  1:110873/251‑1 (MQ=255)                                                                                                                                                                                 cGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCg                                                                   <  1:37769/248‑1 (MQ=255)                                                                                                                                                                                    cccGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCa                                                             <  2:23178/251‑1 (MQ=255)                                                                                                                                                                                                  aaTCACGGGAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATgccc                                                >  2:120917/1‑250 (MQ=255)                                                                                                                                                                                                   aTCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACTGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCCGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATgccc                                                >  2:224262/1‑249 (MQ=255)                                                                                                                                                                                                     cACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCg                                            >  1:248334/1‑251 (MQ=255)                                                                                                                                                                                                           aaaCAGCAGAACATAGAGTAGCAGCAGAATTAAGGCGGTGACCGCGCCGTCCTGGCTGCGAAAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTCCGGACATCATCGTTACCGAACAACTCCACCACGCCCAAATAAACAATCGACAGCACGCAGGCAAAAGCCAGTGCCTGCGCCAGTGGACGTTGCAGAGCGCGCCGCTCGACCACTTTAAAGCGCGCCCAGCCTTAAATGACCAGCGCCTGa                                       >  2:218778/1‑250 (MQ=255)                                                                                                                                                                                                            aaCAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGAc                                      >  1:207881/1‑250 (MQ=255)                                                                                                                                                                                                                           agtagcagcagAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAgcg                      <  2:207881/248‑1 (MQ=255)                                                                                                                                                                                                                           agtagcagcagAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAg                        >  2:233117/4‑249 (MQ=255)                                                                                                                                                                                                                                              cgGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGGGATATCGCGCATc   <  1:321057/251‑1 (MQ=255)                                                                                                                                                                                                                                               ggTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCg  >  1:117139/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  GATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCGCGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCG  >  NC_000913/72754‑73242

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