Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,329,743	G→T	D261Y (GAT→TAT)	dacB →	peptidoglycan DD‑endopeptidase DacB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,329,743	0	G	T	100.0%	257.1 / NA	74	D261Y (GAT→TAT)	dacB	peptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (31/43);  total (31/43)

GCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGT  >  NC_000913/3329503‑3329969                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     gCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTGTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTaaa                                                                                                                                                                                                                            <  2:222955/251‑1 (MQ=255)   cccAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCGTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATGCTGAAATATGAGTTAAAAc                                                                                                                                                                                                                          <  2:54418/251‑1 (MQ=255)        aGCCTGGTGATATGGCTTTTATACGCGTGGCAACTTATTACCCCGTTACGATGGTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGCCCTGAACCTCTTTACGCTGACGGGATGCCTGCCACAACTTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCg                                                                                                                                                                                                                     <  1:462286/251‑1 (MQ=255)         gCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACAGCGTTACGATGTTCAGCCAGGTAAGCAACCTCCCCCGTGGCTCTACCGAAGCGCAATACTGCGAACTGGATGAGGCGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGCTCTGGGCCGCTCCCGATGGTTTTTGCCTTGCAGGCAGGAGCCAGCTATGACGGGGCAATTCTGAACTATGAGTTAAAAAAggggg                                                                                                                                                                                                                    >  2:342067/1‑248 (MQ=255)                  tatGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCt                                                                                                                                                                                                           >  1:488783/1‑251 (MQ=255)                   atGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCt                                                                                                                                                                                                           <  2:443500/250‑1 (MQ=255)                    tGGCGTTTATACGCGTGGCATCTTATTTCCCCGTTACGAGGTGCAGCCAGGTACTCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATCCTGCGAACTGGAGGTGGTGCCAGGCGACCTGATCCTCTTTACGCTGGCGGGATGTCTGCCACAACGTTCGGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGTGCCCGCTATGCCGGTGCAATGCTGAAATATGAGTTAAAACAGGCGGGTATCACCTg                                                                                                                                                                                                          <  2:342458/250‑1 (MQ=255)                     ggCTTTTATATGCGTGGTTCCTTATTTCCCCGTTACTATGTTCAGCCAGGTAAGCACCCGCCCCCGTGGGTCTGCCGAAGCGAAATTCTGCAAACTGGATGTTTTGCAAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCGATGCCGGTTCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTgg                                                                                                                                                                                                         <  2:315487/250‑1 (MQ=255)                     ggCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGa                                                                                                                                                                                                                                  >  1:61308/1‑225 (MQ=255)                     ggCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTgg                                                                                                                                                                                                         <  1:305469/250‑1 (MQ=255)                     ggCTTTGATACGCGTGGCATCGTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGTCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTGTGCCGTGCAGGATGGAGCCAGCTAGGCCGGTTCAATTCTGAAATATGa                                                                                                                                                                                                                                  <  2:61308/225‑1 (MQ=255)                           tataCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATTACCTGGAGCGGa                                                                                                                                                                                                   >  2:52794/1‑250 (MQ=255)                              ccgcgTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑acac                                                                                                                                                                                              <  1:258923/250‑1 (MQ=255)                               cgcgTGCCATCTTAGTACCCGGGTACGATGTTCAGCCAGGTACGCACCCGCCCCCGTGGTGCGGCGGAAGCGCAATACTGCGAACTGGATGGGGTGTTAGGCGACCGGACCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGATCCGCTCCCGTTGGCGTGGACCGTGCAGGATGTGGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACACTCGGGTaggacgtcgagcgcaacact                                                                                                                                                                                              <  2:540678/251‑20 (MQ=255)                                 cgTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGTGTGGA‑ACActgc                                                                                                                                                                                           >  1:146988/1‑251 (MQ=255)                                 cgTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACActgc                                                                                                                                                                                           <  1:472972/251‑1 (MQ=255)                                           tattaCCCCGTTACGATGTTCAGCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAg                                                                                                                                                                                   <  2:200858/249‑1 (MQ=255)                                               aCCCCGTTCCGATGGTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCa                                                                                                                                                                              <  2:437017/250‑1 (MQ=255)                                                 cccGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGCTGTGGCGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTACCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCCGACTCAGGt                                                                                                                                                                           >  1:265237/1‑251 (MQ=255)                                                    gTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGACGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGG‑ACACTGCTGCGCCAGACTCAGGTTa                                                                                                                                                                         >  1:303296/1‑250 (MQ=255)                                                    gTTACGATGTTCAGCCAGGAACGCACCCTCCCCCGTGGTTCTGCCGACGCGCAATACTGCGAACTGGATGTGGTGCCCGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTTTCCCCTGGGGCGGA‑ACCCTGCTGCGCCAGACTccaggta                                                                                                                                                                         >  2:546740/1‑244 (MQ=255)                                                     ttACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGG‑AAActgcggggccctaaccagggtaac                                                                                                                                                                       >  2:452770/1‑231 (MQ=255)                                                             ttCAGCCAGGTAAGCACCCTCCCCCGTGGATCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATAACCTGGAGCGGA‑ACACTGCTGCGCCCGACTCAGGGTAACGAACCTgg                                                                                                                                                               >  1:173185/1‑251 (MQ=255)                                                                aGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAAc                                                                                                                                                            <  1:292195/251‑1 (MQ=255)                                                                aGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAAc                                                                                                                                                            >  1:208170/1‑251 (MQ=255)                                                                aGCCAGGTACGCACCCTCCCCCGGGGTTCTGCCGAAGCGCTATTCTGCGAACTTGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGAGGCCTGACACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGAGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAAc                                                                                                                                                            <  2:112315/251‑1 (MQ=255)                                                                 gCCAGGTACGCCCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGTTGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACg                                                                                                                                                           <  1:63881/251‑1 (MQ=255)                                                                 gCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACg                                                                                                                                                           <  1:4616/251‑1 (MQ=255)                                                                  ccAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGAGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGG‑ACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACgg                                                                                                                                                          >  2:83320/1‑251 (MQ=255)                                                                     gtgactcACCCTCCCCCGTGGTTCTGCCGAAGTGAAATACTGCGAACTGGATGTGGTGCCAGGCGACCTTAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAg                                                                                                                                                       <  2:11136/245‑1 (MQ=255)                                                                         cGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGc                                                                                                                                                  <  2:473529/251‑1 (MQ=255)                                                                              ccTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTa                                                                                                                                             <  2:16150/251‑1 (MQ=255)                                                                                 cccccGTGGTTATGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGATTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAAc                                                                                                                                          <  2:367553/251‑1 (MQ=255)                                                                                      ggggTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCg                                                                                                                                     <  2:98348/249‑1 (MQ=255)                                                                                      gTGGTTCTGCAGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGATGCAATTCTGAAATATGAGTTAAAACAGGCGGGAATCACCTGGAGCAGATAC‑CTGCTCCGCCCGACCCCGGGTAAAGAACCTGGGACGGAAG‑ttacccggcaaaaggcg                                                                                                                                     >  1:433849/1‑236 (MQ=255)                                                                                        ggTTCTGCCGAAGCGCAATACTTCGACCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGTTATCACCTGGAGCGGA‑ACACTGCTGCGCCATACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGc                                                                                                                                   <  1:55107/251‑1 (MQ=255)                                                                                        ggTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGc                                                                                                                                   >  1:106147/1‑251 (MQ=255)                                                                                         gTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGCGAGGA‑AAACGGCTGGGCCAGAATCAGGTTAACGAAACGGGAACGGAAG‑TTGACAGAAAACAGTCggac                                                                                                                                  >  1:476964/1‑249 (MQ=255)                                                                                              gccgacggccaatagtccGAACGGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACTTTCTGATCCGCTCCCGGTGGCTTTTGCCGTGCAGTATGGAGCCAGCTATGCCGGTGCAAGTCTGAAATATGAGTTAAAACAGTCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCg                                                                                                                               <  1:524921/232‑1 (MQ=255)                                                                                              gCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGcccc                                                                                                                                <  1:256672/248‑1 (MQ=255)                                                                                              gCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCt                                                                                                                             <  1:406514/251‑1 (MQ=255)                                                                                                  aagcgaaatccggCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGCGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCAc                                                                                                                         <  2:488268/239‑1 (MQ=255)                                                                                                     cgcAATACTGCGAACTGGAGGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGAt                                                                                                                      <  2:109072/251‑1 (MQ=255)                                                                                                           aCTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAAc                                                                                                                                                            <  2:475968/208‑1 (MQ=255)                                                                                                              gCGAACTGGATTTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGACCGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTaa                                                                                                              <  1:374577/250‑1 (MQ=255)                                                                                                               cGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGTGCGGA‑ACACTGCTGCGCCAAACCCAGGTTTACGAACCTGGAACGGGAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCCCGATCTGCTTaaa                                                                                                             >  1:153312/1‑249 (MQ=255)                                                                                                                gAACTGGATGTGGGGCCAGGCGCCCTAAACCGCTTTACGCTGCCGGGATGCATGCCACAACGTTCTGAGCCGCTCCCGATGGCTTTTGCAGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGAGGGGATCACATGGGGCGGA‑ACAATGCTGCGCCAGACGCCGGGTAAAGAACCTGGAACGGGCG‑TTGGCCGTAAAAAGTCGGGCCCGCTGGACGagctcgttatgat                                                                                                           >  2:345625/1‑239 (MQ=255)                                                                                                                     ggATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATg                                                                                                       <  1:26041/250‑1 (MQ=255)                                                                                                                     ggATGTGGGGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGACACAACGTTCGGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTAt                                                                                                        <  1:546740/249‑1 (MQ=255)                                                                                                                         gtggtgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGGGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCCGCCCGATCCGCTTAAAATTTTGCTgg                                                                                                   >  1:35915/1‑249 (MQ=255)                                                                                                                         gtggtgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGA‑ACAATGCTGCGCCAGACTCAGGTTAACGAAACAGGGACGG‑AGAGTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGa                                                                                                   >  1:120107/1‑250 (MQ=255)                                                                                                                                cAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGCTCACCTGGGGCGGG‑AAACTGCTGCGCCAGACTAAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGCACGCCACGCCGCACGGACTGCTTAAGAGTATGCTGAAAAAga                                                                                             >  2:460542/1‑248 (MQ=255)                                                                                                                                    cGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGaca                                                                                       <  2:208170/251‑1 (MQ=255)                                                                                                                                             ccGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAACGGA‑ACCCTTCTGCGCCAGACTCAGGTTAAAGAACCTGGGACGGGAG‑GTTCCAGTAAACAGTCGGCCCCGCTGCACCATCTGCTTAAGATTATGCTGGAAAAGTCGGACAACATGATCg                                                                              >  2:542632/1‑251 (MQ=255)                                                                                                                                                cTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑GTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCg                                                                           >  1:239637/1‑251 (MQ=255)                                                                                                                                                 tttACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGAGCTGCTTAAGATTATGCt                                                                                                     <  1:128316/224‑1 (MQ=255)                                                                                                                                                 tttACGCTGACGGGATGCCGGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACGGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGaaa                                                                                                 <  1:78441/228‑1 (MQ=255)                                                                                                                                                   tACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGAt                                                                         <  1:161144/250‑1 (MQ=255)                                                                                                                                                     cGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGG‑ACACTGCTGCGCCAGACTCAGGTTAAAGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACg                                                                      >  1:206448/1‑251 (MQ=255)                                                                                                                                                      gCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATAACCCGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGGACGGGAG‑ATGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACCTGATTGCCCATaagg                                                                     >  1:59372/1‑248 (MQ=255)                                                                                                                                                          aCGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGCGCGGA‑ACACTGCCGCGCCAGAATCAGGTTAACGAACCTGGAACGGCAG‑TTGCCAGAAAACAGTAGGCCCCGCTGCACGATCTGCTTAAGATTATGGTGAAAAAGTCGGACAACATGATCGCCGCTACGGtttt                                                                 >  1:463495/1‑251 (MQ=255)                                                                                                                                                           cGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑GTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGGCGGACAACATGATCGACGATACGGTTTTc                                                                >  1:393839/1‑251 (MQ=255)                                                                                                                                                                             gTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATgcg                                              <  1:248152/251‑1 (MQ=255)                                                                                                                                                                              ttCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGGTTTCCGCATGATAAGCCATgcgc                                             >  1:141333/1‑251 (MQ=255)                                                                                                                                                                                  gAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCt                                          <  2:106147/250‑1 (MQ=255)                                                                                                                                                                                            cGTTGGCTTTTGCCGGGCAGGATGGAGCCAGCTATGCCGGAGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGG‑ACACTGCTGCGCCCGACTCAGGGTAACGAAACTGGAACGGTAG‑TTGCCAGGAAACAGACGGCCCCGCTGCAAGATCTGCTTAAGATTATGCTGAAAAAGACGGACAACATGATCGCCTCAACGGTTTTCCGCATGAAAGGCCACGCGCGCTTCAAAGTGCCt                               >  2:292050/1‑251 (MQ=255)                                                                                                                                                                                              ttGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATTCTGAAAAAGGCCGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCCTGCGCGCTTCAATGTGCCTgg                             >  1:134867/1‑251 (MQ=255)                                                                                                                                                                                                     ttGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCGCCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATg                       <  2:516885/250‑1 (MQ=255)                                                                                                                                                                                                           ggCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGcc                <  2:120502/250‑1 (MQ=255)                                                                                                                                                                                                                 aTGGAGCCAGCTAGGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCCCCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGGCGGCCCCGCTGCACGATCTGCTGAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTGCAATGTGCCTGGAACATGGCGGTCCGGGTCg          <  2:433849/251‑1 (MQ=255)                                                                                                                                                                                                                       ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGcc    <  2:134867/251‑1 (MQ=255)                                                                                                                                                                                                                         aGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGCTACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGGCCGGCCGGACGCCGt  >  1:558146/1‑251 (MQ=255)                                                                                                                                                                                                                                       aaTTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCAt                                                 <  1:193238/190‑1 (MQ=255)                                                                                                                                                                                                                                       aaTTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCAt                                                 >  2:193238/1‑190 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     GCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGA‑ACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAG‑TTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGT  >  NC_000913/3329503‑3329969

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