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%	230.3 / NA	60	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/27);  total (33/27)

ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCC  >  NC_000913/1652087‑1652571                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  aTGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCCCCCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTTATGTTCACCTCag                                                                                                                                                                                                                                               >  1:261098/1‑249 (MQ=255)   ggTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagaga                                                                                                                                                                                                                                            <  1:432984/250‑1 (MQ=255)            cTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGc                                                                                                                                                                                                                                     >  1:181674/1‑248 (MQ=255)                ttATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAAGTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGATCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcaga‑aca                                                                                                                                                                                                                              <  2:584774/250‑1 (MQ=255)                  ataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGAGCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCCGCAATCGACGCTCTCCGGATTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGCATGGcaga‑acagaa                                                                                                                                                                                                                           <  2:81295/251‑1 (MQ=255)                                        gCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGGCATTCATCTTATTAAACCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATTGCAGA‑ACAGAAAAAAAAAAATGCACCTTTGttt                                                                                                                                                                                                     >  1:338842/1‑251 (MQ=255)                                           gtgtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGttttt                                                                                                                                                                                                   <  2:345420/250‑1 (MQ=255)                                                cACTTGGATGGGCGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAA‑CAATCGGGTAATTCATCTTATTAAGCCAACAATCGACGCCCTCCGGAGTCAGATGACATTAACGAGAATGAGCAAAAAGCATATCATTTATTTTCACCTCAGAGAGCATGGCAGA‑AAAAAAAAACAAAAATTAACCTTTGATTTTCAAcac                                                                                                                                                                                             >  2:393679/1‑249 (MQ=255)                                                 aCTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGAAGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                                                                                                                                                                            <  2:61928/251‑1 (MQ=255)                                                         gggAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCa                                                                                                                                                                                    <  2:261098/251‑1 (MQ=255)                                                          ggaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAg                                                                                                                                                                                   <  1:510580/251‑1 (MQ=255)                                                           gaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCAGTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc                                                                                                                                                                                  <  2:408503/251‑1 (MQ=255)                                                            aggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAA‑CAATCGGGTCATTCATCTTATTAAGTCAGCAAACGACGATTTCCGGAGTCAGAAGACATTAAAGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGGAGA‑ACAGAAAAAAAAAAATGAACCTTTGTTTTTCAAACTGAAgtgttag                                                                                                                                                                                   >  1:531089/1‑246 (MQ=255)                                                              gaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTAACGGTACCAAAAACACAGGCGGGCACACAA‑AGGGTCAGTTATATTATTAAACCAGCAAAAGAGGCTCTCCCGAGACAGATGGCATTAACAAAAATGAAGAAAGAACATATCATTGATGGTCAACTCAGAAAGTATGGAAGG‑ACAGGAAAAAAAAAAAACACCCTTGTTTTTCAAACAGAAgttgacgcgcg                                                                                                                                                                               >  2:519565/1‑242 (MQ=255)                                                                     aTCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCCGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGGGTCAGCGAGAGTaaa                                                                                                                                                                         >  1:555604/1‑250 (MQ=255)                                                                         gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                                                                                                                                                                         <  2:18931/246‑1 (MQ=255)                                                                            aaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGGATGGCAGA‑AAAGAAAAACAAAAATGCACCTTTGTTTTTCAAACTGAAGTGTCAGCGAGAGTAAAAAATTATg                                                                                                                                                                 >  1:576409/1‑251 (MQ=255)                                                                            aaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGGCAGATGACATTAACGAGACTGAGCAAAGAGCATATCCTTTATGTTCACCTCAGAGGGTATGGCAGA‑AAAGAAAAACAAAAATGCACCTTTGTTTTTTAACCTGAAGTGTCAGCGAGAGTAAAAAATTATg                                                                                                                                                                 >  1:563538/1‑251 (MQ=255)                                                                                ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTg                                                                                                                                                              <  2:181674/250‑1 (MQ=255)                                                                                 gAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGGATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc                                                                                                                                                            >  1:506551/1‑251 (MQ=255)                                                                                  aaCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAcc                                                                                                                                                           <  2:102428/251‑1 (MQ=255)                                                                                     aaTGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTTACCTCAGAGAGAATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGGGAGGAAAAAATTATGGTGACCAtt                                                                                                                                                        >  1:518673/1‑251 (MQ=255)                                                                                      atgatgATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTCCAAAAAACACAGGCGGGAAC‑CAAACGGGGCATTCATCTTATTAAACCAGCAATCGACGCCCTCCCCAGGAAGATGACATTAAACAGACTGCGCAAAGAGAAAATCATTGTTGTTCACCTCAGAGAGAAAGGCAAACACAAAAAAAAAAAAAACCAC‑TTTGTTTTTAAACCTGAAGTGTCAGAGAGAGCAAAAAATTATGGTGaaccttt                                                                                                                                                       >  2:533004/1‑245 (MQ=255)                                                                                      atgatgATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCCTTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATTACATTAACGAGAATGAGCAACGAGCATATCATTGATGTTCACCTCAAAGAGTATGGGAGA‑ACAGAAAAAAAAAAATGCACATTTTTTTTTTAAAATGAAGTGTACGCGAGAGTAAAAAATTAAGTTGACCAttt                                                                                                                                                       >  1:436020/1‑251 (MQ=255)                                                                                        gatgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAtttt                                                                                                                                                      <  1:257798/250‑1 (MQ=255)                                                                                         atgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAA‑CAAACGGGGCATTCATCTTATTAAGCCAGCAAACGCCGCCCTCCGGAGTCAGATGACATTAAAGAGAATGAGCAAAGAGCATATCATTGATGTTCACCTCCGAGAGGATGGCAGA‑ACAGAAAAACAAAAATACACCTTTGTTTTTCAACCCTAAGTGTCAGAGAAAGTAAAAAAATATTGGGAACATTTTAc                                                                                                                                                    >  2:15973/1‑251 (MQ=255)                                                                                           gatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCTCCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAcc                                                                                                                                                   <  2:608291/250‑1 (MQ=255)                                                                                            atCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATAGGGCCATTCCTCTTCTTTAGCCAGCAATCAACGCTCTCCGGAGTGACATGAAATTAACGAGACTGAGCAAAAAGCCTATAATTGATGTTCACCTCCGCGAGAATGGCAGA‑ACAGCAAAACAAAAATGCACCTTTGTTTTTTAAACTGCAGTGTCAGCGAGAGAAAAAAAATTTGGTGAACCTTTTACCgc                                                                                                                                                 >  1:236501/1‑250 (MQ=255)                                                                                                     aaTTTAGCAAAAGACCGTTTCACGGTACAAAAAACACAGGCGGGAAC‑CAATCGGGGCATTCAACTTATTAAGCCAGCAATCGACGCTCTCCGGGGTCCGATAACATTAACGAGACTGAGCAAAAAGCATAACATTGATGTTCAACTCAGGGAGTATGGCAGA‑ACAGAAAAAAAAAAAAGCACCTTTGTTTTTCAACCTTAAGTGGTAGCGGGAGTAAAAAAATAAGGGGAACATTTTACAGTTTGCTCAAt                                                                                                                                        >  2:466540/1‑251 (MQ=255)                                                                                                         tAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAg                                                                                                                                     <  2:142542/250‑1 (MQ=255)                                                                                                           gCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGCCATTCCTCTTATTAAGCCAGACATCGATGCTCCCCGGAATAAGATGACAATAACGCGACCGAGCAAAGAGCAAATAATTGATGTTTACCTCCGAGAGTATGGCAGA‑ACTGAAAAACAAAAATGCACCTTTGTTTTTTAAACTGAATGGTCAGCGAGAGGAAAAAATTATGGTGACCCTTTTAccagttaattaataaggc                                                                                                                                   >  1:475930/1‑234 (MQ=255)                                                                                                             aaaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGa                                                                                                                                <  2:133834/251‑1 (MQ=255)                                                                                                                     gTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGAc                                                                                                                        <  2:576409/251‑1 (MQ=255)                                                                                                                                  aaaaaCACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGTACGCAGCGATAAAAc                                                                                                           <  2:259266/251‑1 (MQ=255)                                                                                                                                       cacaGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGcc                                                                                                      >  1:278862/1‑251 (MQ=255)                                                                                                                                         caGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCgg                                                                                                    >  2:101966/1‑251 (MQ=255)                                                                                                                                          aGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGa                                                                                                   <  2:506551/251‑1 (MQ=255)                                                                                                                                               ggAAC‑CAATCGGGTCATTCCTCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTTACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGAATCCg                                                                                              >  2:40989/1‑251 (MQ=255)                                                                                                                                                 aaC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTAAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGcc                                                                                            >  1:121847/1‑251 (MQ=255)                                                                                                                                                         cGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaa                                                                                      >  2:549741/1‑250 (MQ=255)                                                                                                                                                                 tcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCatat                                                                              <  1:483367/250‑1 (MQ=255)                                                                                                                                                                   atCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATc                                                                             >  2:571646/1‑249 (MQ=255)                                                                                                                                                                      ttattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACTTGCCGGACTCCGCCATCGAAAATCATATCAGTCg                                                                        >  2:172058/1‑251 (MQ=255)                                                                                                                                                                           aaGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGAc                                                                    <  1:120511/250‑1 (MQ=255)                                                                                                                                                                           aaGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACa                                                                   <  2:518673/251‑1 (MQ=255)                                                                                                                                                                                       gACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctg                                                       >  1:113555/1‑251 (MQ=255)                                                                                                                                                                                       gACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctg                                                       >  1:130478/1‑251 (MQ=255)                                                                                                                                                                                          gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgctg                                                    >  1:78037/1‑251 (MQ=255)                                                                                                                                                                                          gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCct                                                        >  2:478584/1‑247 (MQ=255)                                                                                                                                                                                          gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCct                                                        <  1:478584/247‑1 (MQ=255)                                                                                                                                                                                             ctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTc                                                 >  1:545028/1‑251 (MQ=255)                                                                                                                                                                                                       aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCt                                       <  1:15973/251‑1 (MQ=255)                                                                                                                                                                                                                tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAc                              >  1:3127/1‑251 (MQ=255)                                                                                                                                                                                                                 aaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAcc                             <  1:201460/251‑1 (MQ=255)                                                                                                                                                                                                                    gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTACCCCgg                          >  2:430571/1‑251 (MQ=255)                                                                                                                                                                                                                    gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCgg                          >  1:130012/1‑251 (MQ=255)                                                                                                                                                                                                                    gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCgg                          <  1:533004/251‑1 (MQ=255)                                                                                                                                                                                                                           aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat                   <  1:101966/251‑1 (MQ=255)                                                                                                                                                                                                                                  agCATATCATTGATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAc            >  2:378642/1‑251 (MQ=255)                                                                                                                                                                                                                                             gATGTTCACCTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGcc  <  2:465431/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAC‑CAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGA‑ACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCC  >  NC_000913/1652087‑1652571

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: ‑