Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,953,546	G→A	Q35Q (CAG→CAA)	cmoB →	tRNA U34 carboxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,953,546	0	G	A	100.0%	171.9 / NA	47	Q35Q (CAG→CAA)	cmoB	tRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/26);  total (21/26)

GCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCC  >  NC_000913/1953304‑1953788                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   gCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGgtt                                                                                                                                                                                                                                            <  1:216578/251‑1 (MQ=255)              aTGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTgg                                                                                                                                                                                                                               <  1:229877/251‑1 (MQ=255)                gCTGACCTATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAGCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGTTAACTTTTATTCTCTGATTGCCTAAAATCATCTTTCTCACTGGCTCGAAACGCTGCTCGTGCAAATTGCTAACTTGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTc                                                                                                                                                                                                                             <  1:380494/251‑1 (MQ=255)                gCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTc                                                                                                                                                                                                                             >  2:8149/1‑251 (MQ=255)                gCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTc                                                                                                                                                                                                                             <  2:138213/251‑1 (MQ=255)                     ccGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATg                                                                                                                                                                                                                        <  2:103536/251‑1 (MQ=255)                       gATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCg                                                                                                                                                                                                                      >  1:359278/1‑251 (MQ=255)                        aTTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATgagg                                                                                                                                                                                                                     >  2:230554/1‑248 (MQ=255)                               ggAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGtt                                                                                                                                                                                                              <  2:277447/251‑1 (MQ=255)                                     ccATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTAcc                                                                                                                                                                                                        >  2:333650/1‑251 (MQ=255)                                      cATAAAGCACGCCTGCATAACTCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCt                                                                                                                                                                                                       <  2:133003/251‑1 (MQ=255)                                      cATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCt                                                                                                                                                                                                       >  1:61572/1‑251 (MQ=255)                                                gCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAAc                                                                                                                                                                                             <  1:184957/251‑1 (MQ=255)                                                    gCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGt                                                                                                                                                                                          >  1:22688/1‑250 (MQ=255)                                                     cATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTAt                                                                                                                                                                                        >  2:149186/1‑251 (MQ=255)                                                                 tttGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGGTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTGCCTGAAATTAAACCGTATCGTCTTGAttt                                                                                                                                                                             <  1:333958/250‑1 (MQ=255)                                                                    gAGCATAGTAATCTGTGGTTCCAATTCTTTAACTTTGGTTCACTGATGGCATTAAAAGCAGAGGACGCTGCATTGGCGGCTTTGTTAACTTTTATTCTCTTATTGCCAAAAATCATCTTTCTCACTGGCGCGAAACGCTACCCGCGCAGATTGCTAACTGGCGGCGCGAGCAGCAACGCGGGTTGTTTAAGCAGTGGTCAAATGCAGTGGACTTTCTACCGGAAATTAAACCGTATCGTCTTGATTTAtt                                                                                                                                                                          <  1:388528/250‑1 (MQ=255)                                                                     aGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGc                                                                                                                                                                        >  1:220799/1‑251 (MQ=255)                                                                     aGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGc                                                                                                                                                                        <  2:137657/251‑1 (MQ=255)                                                                     aGCATAGTGAGCTGAAGTACCAGTGCTTTAACTTTGGTTCACTGGAGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTCACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGCGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGc                                                                                                                                                                        >  2:51008/1‑251 (MQ=255)                                                                       cATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCAt                                                                                                                                                                      >  1:294358/1‑251 (MQ=255)                                                                                       tCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGaa                                                                                                                                                       >  1:52292/1‑250 (MQ=255)                                                                                                      ttGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTg                                                                                                                                        >  2:346416/1‑250 (MQ=255)                                                                                                        ggTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAAGGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAg                                                                                                                                      <  2:193534/250‑1 (MQ=255)                                                                                                           tCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGGTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAATCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGCAGAGCCACTGAgcgc                                                                                                                                   >  2:227246/1‑250 (MQ=255)                                                                                                            cACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAg                                                                                                                                 <  1:367516/251‑1 (MQ=255)                                                                                                             aCTGGTGGCAGTAAAAGCAGATGACTCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTTCATAGCGTTACGGCAGAAAGCGAAGAGCCACGGAGCGCAgg                                                                                                                                <  2:225195/251‑1 (MQ=255)                                                                                                                ggtggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAggg                                                                                                                               <  2:125817/249‑1 (MQ=255)                                                                                                                   ggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGGGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAAt                                                                                                                          >  1:313709/1‑251 (MQ=255)                                                                                                                          aaaGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGc                                                                                                                   >  1:135205/1‑251 (MQ=255)                                                                                                                                gagGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGa                                                                                                              >  1:315244/1‑250 (MQ=255)                                                                                                                                gagGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGGTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGa                                                                                                              >  1:64204/1‑250 (MQ=255)                                                                                                                                         gCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATg                                                                                                    <  2:22688/251‑1 (MQ=255)                                                                                                                                              aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGCGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCa                                                                                                >  2:187722/1‑250 (MQ=255)                                                                                                                                                       ggTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGcc                                                                                      >  1:323949/1‑251 (MQ=255)                                                                                                                                                                                                  ggCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTGCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCTTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCg                                             <  2:35534/249‑1 (MQ=255)                                                                                                                                                                                                         aaCGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGAGTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCGCACTGTATGGCGTCAACATCGATACCGAAt                                    <  2:295700/251‑1 (MQ=255)                                                                                                                                                                                                            gCTGCCCGCGCAGATTGCTAACTGGCAGCTCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGc                                 >  2:129106/1‑251 (MQ=255)                                                                                                                                                                                                                     gcAGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACt                        <  1:8149/251‑1 (MQ=255)                                                                                                                                                                                                                       aGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTgg                      <  1:227246/251‑1 (MQ=255)                                                                                                                                                                                                                       aGATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGGTCCGACTgg                      <  1:396421/251‑1 (MQ=255)                                                                                                                                                                                                                        gATTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTgg                      >  1:334183/1‑250 (MQ=255)                                                                                                                                                                                                                         aTTGCTAACTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGaa                    <  2:323949/251‑1 (MQ=255)                                                                                                                                                                                                                               aaCTGGCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATggg               <  1:64156/250‑1 (MQ=255)                                                                                                                                                                                                                                   ggCAGCGCGAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCgt          <  1:79908/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggCAGCGCGAGCAACAACACGGGTTGTTTAAGCAGGTGTCAAATGCGGTGGAGGTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGTAAATGGGATc            <  2:52333/249‑1 (MQ=255)                                                                                                                                                                                                                                           gAGCAGCAACACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGcc  <  1:230554/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   GCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGGCGTCAACATCGATACCGAATGGCGTTCCGACTGGAAATGGGATCGTGTTCTGCC  >  NC_000913/1953304‑1953788

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