Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,253	C→T	V423V (GTG→GTA)	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,253	0	C	T	100.0%	156.7 / NA	49	V423V (GTG→GTA)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/17);  total (32/17)

CTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACA  >  NC_000913/73007‑73491                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                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)                           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)                                         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)                                                                                                                                         cGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑Ca                                                                                                     >  2:269170/1‑250 (MQ=255)                                                                                                                                                     cgcAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCAt                                                                                        >  2:33571/1‑251 (MQ=255)                                                                                                                                                            cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCa                                                                                   <  2:308006/247‑1 (MQ=255)                                                                                                                                                               ccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGcgcc                                                                              >  2:309716/1‑248 (MQ=255)                                                                                                                                                                  cTTTTAAGCGCGACCCGCCTTCAATGCCCAGCGACTGACACAAAATACTGTAGAGGGCGGTGATATCGCGCATCGGGTTTTCAAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATACCGTCAACAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACCATCACCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTg                                                                            >  1:182538/1‑250 (MQ=255)                                                                                                                                                                                        aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCcgcg                                                     >  1:336400/5‑251 (MQ=255)                                                                                                                                                                                        aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGAGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCGCATCACCAGCGCCTGACCCGCCAGCCTTTTCTGCcgcg                                                     >  1:356498/5‑251 (MQ=255)                                                                                                                                                                                              cAGCGACTGAACCAACATACTGGAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGCCAATCACAATGCCGTCCGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGGAGCCAGCACACTCCCCGGCAAGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCCGACCCGCCAGCCTTTTCTGCCGCGCCCgac                                               >  2:240693/1‑249 (MQ=255)                                                                                                                                                                                                  gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTc                                           >  1:262059/1‑251 (MQ=255)                                                                                                                                                                                                                         ggcggTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGcatc                     >  2:124723/1‑250 (MQ=255)                                                                                                                                                                                                                           cggTGATATCGCGCGTCGGGTTCTCCAGTACTTTCAGAGCATAAGGGATCGCCATTAAAGCATTGGTGAAAATCAAACTGCCGTCAGCATATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGACAGCACAATCCCCGGCATGGCGAGGATACACCA‑GCCGCTCATCTCCAGAGCCTGACCCGCCAGCATTTTCTTCCGCGCCCGAAGTTCAAGACTGCGCCATCGCCGCATCATg                  >  2:243883/1‑251 (MQ=255)                                                                                                                                                                                                                              tGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATACCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATACCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTc               >  2:185244/1‑251 (MQ=255)                                                                                                                                                                                                                                atatCGCCCATCGGGTTTGCGAGTACTTGCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGACAGCAGATTGCGGCAGCCCAAGGGTGTTGTTGACAAGAAAAAAAAAGCCGTTAGCCAGCACAACCCCCGGCTTGGCGAGGAT‑CAACATGCCGCGCATCTCCAGCCCCTGACCCGCCATCATTTGCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAg             <  1:32135/251‑1 (MQ=255)                                                                                                                                                                                                                                       gcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTAcc      >  2:230326/1‑250 (MQ=255)                                                                                                                                                                                                                                          tCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGcc                                                        >  1:109276/1‑198 (MQ=255)                                                                                                                                                                                                                                           cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTACCGCATTGGTGAACATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGCT‑AAACATGCCGCTCCTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCCGCAGAAGTTCGCGACTGCTCCCTAGCAGCATCATGGACAGCACAACGccc   >  2:176587/1‑248 (MQ=255)                                                                                                                                                                                                                                            gggTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGcaca  >  1:191797/1‑250 (MQ=255)                                                                                                                                                                                                                                                tttCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGcc                                                                     >  2:23811/1‑179 (MQ=255)                                                                                                                                                                                                                                                tttCCAGTACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCGCCTGACCCGcc                                                                     <  1:23811/179‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                CTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCATTAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTGTTGTTGAGCAGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGAT‑CAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCACA  >  NC_000913/73007‑73491

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