Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,331	T→C	intergenic (+2346/+397)	ydfD → / ← ynfP	lysis protein/protein YnfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,652,331	0	T	C	100.0%	169.3 / NA	46	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/17);  total (29/17)

TATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACTTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652103‑1652575                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      tataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCATTCATCTTATTTAGCCAGCAATA‑CA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCCCCTCAGAGA‑GTATgggcgaac                                                                                                                                                                                                                                   >  2:337344/1‑242 (MQ=255)                       aGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAAAGAGACTGAGCAAAGAGCATATCAAT‑GATGTTCAACTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                                                                                                                                                                                                          >  2:285080/1‑251 (MQ=255)                           gtgtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACATAAAAACAAAAATGCACCTTTGTTTTTc                                                                                                                                                                                                      <  1:392828/251‑1 (MQ=255)                           gtgtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTc                                                                                                                                                                                                      <  1:337344/251‑1 (MQ=255)                                cACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                                                                                                                                                                                 <  2:32603/251‑1 (MQ=255)                                                cATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGTAACCAAACGGGACATTCATCCTATTAAACCAGCAAAC‑GA‑CGCCCTCCGGAGGCAGATGAAATTAACAAGACTGAACAAAAAGCATACCATTTGAT‑TTCAACTCA‑AGAAGTAATGCAGAACAGAAAAACAAAAATTCACCTTTGTTTTTCAACATGAAGGGTCAACgagg                                                                                                                                                                                  >  2:29164/1‑249 (MQ=255)                                                   aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa                                                                                                                                                                               <  2:405388/250‑1 (MQ=255)                                                    gATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                                                                                                                                                                              >  1:127488/1‑250 (MQ=255)                                                     aTCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGTACCAATAGGGGAATTCATATTAATAAGCCCACAATCAGAGCGCTC‑CCGG‑GCCAGACGACATTTACGAGACTTAGCAAAGAACATATCCTT‑GATGTTTACCTCAGGGA‑GTTTGGCAGAAAAGAAACACAAAAATGCAACTTTGTTTTTCAAACTGCAATGTCAGCGAGAGTaaa                                                                                                                                                                             >  2:234473/1‑250 (MQ=255)                                                                ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCCCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑TATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGGAAAAAATTATGGTg                                                                                                                                                                  >  1:375549/1‑250 (MQ=255)                                                                ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGAACCAAAAACACAGGCGGGAAACAATCGGGTCATTCATCTTATTAAGACAACAATC‑GG‑CCCTCTCCGGAGGCAGATGACATTAAAGAGAATGAGCAAAAAGCATATCCTT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAAAGCACCTTTGTTTTTCAACCCTAAATGTCCGCGAGAGGAAAAAAatattggt                                                                                                                                                                  >  2:147604/1‑246 (MQ=255)                                                                  aaCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAcc                                                                                                                                                               <  2:201836/251‑1 (MQ=255)                                                                   aCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCATCGGGGCATTCATCTTTTTTAGCCAGCAATC‑GA‑CGCTCCCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATT‑GATGTTTACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGGGTCAGCGAGAGTAAAAAATTATGGTGAccc                                                                                                                                                              >  1:135119/1‑250 (MQ=255)                                                                         atgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCCGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAAAAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAACGAGAGTAAAAAATTATGGTGACCATTTTAc                                                                                                                                                        >  1:434433/1‑251 (MQ=255)                                                                         atgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGcagaacaga                                                                                                                                                                                                                                <  1:108250/179‑1 (MQ=255)                                                                           gatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCg                                                                                                                                                      <  2:135119/251‑1 (MQ=255)                                                                                   gaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCa                                                                                                                                              <  2:74125/251‑1 (MQ=255)                                                                                             aaaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATAGGGGCATTCATCTTATTAAGCCCGCAATC‑GA‑AGCCCTCCGGAGTCAGATGACATTAACGAGAATGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAAAAGAAAAACAAAAATGCACCTTTGTTTTTTAACCTGAAATGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAACGTTGACTCAATAAGGc                                                                                                                                       >  1:224904/1‑248 (MQ=255)                                                                                                 gACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGGGACCATTTTACCGTTGACTCAATAAGGCAGAtgtg                                                                                                                                >  1:23212/1‑251 (MQ=255)                                                                                                      tttCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACCTTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GGATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGgggggcg                                                                                                                           >  1:327776/1‑248 (MQ=255)                                                                                                      tttCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGCGACTGAGCAAAGAGCCTATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTTACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACg                                                                                                                           >  1:27112/1‑251 (MQ=255)                                                                                                                          aGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCgg                                                                                                        >  1:409513/1‑250 (MQ=255)                                                                                                                          aGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGAATATCATT‑GATGTTAACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACATGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCgg                                                                                                        >  1:366169/1‑250 (MQ=255)                                                                                                                          aGGCGGGAAACAATCGGGGCAATCATCTTACTACGCCAGCAATC‑GC‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GCTGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGGTTTTCAACCTGAAGTGACAGCGAGAGTAAAAAATTATGGGGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGAAGCAGCGATAAAACGAGCCgg                                                                                                        >  2:394891/1‑250 (MQ=255)                                                                                                                           ggCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTT‑AACGTTGACTCAATTAGGCAGATGGGGGACGCAGCGATAAAACCTGCCGGACt                                                                                                     >  1:368450/1‑251 (MQ=255)                                                                                                                                  aCCAATCGGGTCATTCATCTTATTAAGCCAGAAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGcc                                                                                                >  1:141845/1‑250 (MQ=255)                                                                                                                                        cGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACACAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaa                                                                                          >  2:283871/1‑250 (MQ=255)                                                                                                                                                      tattaAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCagag                                                                                                                                                                                                                                                 >  2:23027/1‑86 (MQ=255)                                                                                                                                                      tattaAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCagag                                                                                                                                                                                                                                                 <  1:23027/86‑1 (MQ=255)                                                                                                                                                       attaAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGCATCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATAAGTCg                                                                            >  2:389656/1‑249 (MQ=255)                                                                                                                                                           aGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACAt                                                                      <  2:182059/251‑1 (MQ=255)                                                                                                                                                             ccAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATAc                                                                    <  1:56383/251‑1 (MQ=255)                                                                                                                                                                     c‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCct                                                            >  1:85583/1‑251 (MQ=255)                                                                                                                                                                           gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgctg                                                        >  1:414443/1‑251 (MQ=255)                                                                                                                                                                                   gAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGa                                                <  1:33766/251‑1 (MQ=255)                                                                                                                                                                                    aGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGAc                                               <  1:384027/251‑1 (MQ=255)                                                                                                                                                                                              cATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGGGCt                                     >  1:132096/1‑251 (MQ=255)                                                                                                                                                                                                   aCGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGCGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGCCTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGACCTAAccc                                >  2:71276/1‑251 (MQ=255)                                                                                                                                                                                                      agaCTGAGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGc                             >  2:414898/1‑251 (MQ=255)                                                                                                                                                                                                            aGCAAAGAGCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat                       <  2:23212/251‑1 (MQ=255)                                                                                                                                                                                                               aaaGACCATAGCATT‑GATGTTCCCCCCAGATG‑GTATGGCAGAACAGAAAACCAAAAATCCACCTTTGTTTTTCACCCTAAAGTGTCAGCAAGAGAAAAAAATTATGGGGCCCATTT‑ACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTTGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCa                   <  2:368450/251‑1 (MQ=255)                                                                                                                                                                                                                    gCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcc               >  1:157833/1‑251 (MQ=255)                                                                                                                                                                                                                    gCATATCATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCATTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcc               >  2:144203/1‑251 (MQ=255)                                                                                                                                                                                                                          cATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATgg         <  1:29164/251‑1 (MQ=255)                                                                                                                                                                                                                          cATT‑GATGTTCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATg          <  1:34293/250‑1 (MQ=255)                                                                                                                                                                                                                                   ttCACCTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGc  >  1:404467/1‑250 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      TATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATC‑GA‑CGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATT‑GATGTTCACTTCAGAGA‑GTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652103‑1652575

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