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%	201.8 / NA	57	D261Y (GAT→TAT)	dacB	peptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (28/29);  total (28/29)

ACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCG  >  NC_000913/3329498‑3329984                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  aCAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCCGGTACG‑AACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGTAGCCAGCTATGCCGGTGCAATTCTGAAATATg                                                                                                                                                                                                                                                >  1:184097/1‑249 (MQ=255)           aaaGCCTGGAGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAgg                                                                                                                                                                                                                                    >  2:118926/1‑251 (MQ=255)               ccTGGTGATATGGCTGTTATACGCGTGGCATCTTATTACCCCGTTACGATGTCCAGCCAGGTACG‑CGCCCTCCCCCGTGGTTCTGCCGAAGCGAAATACTGCGAACTGGATGTGTTGCCAGGCGACCTGTACCGGTTTATGCTGACGGGTGGCCTGCGACATCGTTCTGAGCCGCTCCCTTTTGCTTGTGGCGTTCATGATGGAGCCAGCTATGCCGGGGCAATTCTGAAATAGGAGTTAAAACATGCggg                                                                                                                                                                                                                                <  2:455031/251‑1 (MQ=255)                               ttATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCgg                                                                                                                                                                                                                 <  2:21403/250‑1 (MQ=255)                                tataCGCGTGGCATCTAATTACCCCGTTCCGATGTTCAGCCAGGTACG‑CACCCCCCCCCGGGGTTCGGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCTCCTTGAGCGGa                                                                                                                                                                                                                <  2:201998/250‑1 (MQ=255)                                     gcgtgttaTCTTATTACCCCGTTACTATGTTCAGCCAGGTACGTCTCCGTCCCC‑GGGGTTCTGCCGAAGCGCAAGAATGTGAACTGGATGTGGTGCTAGGCGACCGGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGTCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActg                                                                                                                                                                                                          <  2:379905/244‑1 (MQ=255)                                      cgTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgc                                                                                                                                                                                                         <  1:341134/251‑1 (MQ=255)                                        tGGCATCTTATTACCCCGTTACGTTGTTCAGCCAGCCACG‑CACCCTCCCCCGCTGTTCTGCCGGCGCGAACTACTGCGAACTGGATGTGGTGCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCATGCTACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGACTCAGATGGAGCGGAACActgctg                                                                                                                                                                                                       >  2:171031/1‑251 (MQ=255)                                         ggCATCTTATTACCCCGTTACGATGTTCACCCAGGTGCG‑CTCCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCTAACGGGATGTGGTGCCAGGCGACCTGTACCGCTTTACGCTGGCGGGATGCCTGCCACAACGTTCTTAGCCTCTCCCGTTGGCTTTTGCCTTGCATGATGGAGCCAGCTATGCCGGTGCAATTCTGATATTTGAGTTAAAACAGTCGGGTATCACCTGGAGCGGAACACtgctgc                                                                                                                                                                                                      <  1:399450/251‑1 (MQ=255)                                                    cccccGTGACGATGTTCAGCCAGGTACG‑CACCCTCCCCCGTGGGTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAg                                                                                                                                                                                           <  1:17642/250‑1 (MQ=255)                                                      cccGTTACGATGTTCAGCCAGGTCCG‑CACCCGCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGTAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGt                                                                                                                                                                                         <  2:234527/251‑1 (MQ=255)                                                          ttACGATGTTCAGCCAGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGGACACTGCTGCGCCCGACCCAGGTTAAc                                                                                                                                                                                     >  1:87496/1‑251 (MQ=255)                                                            aCGATGATCAGACAGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGTGTTAAAACAGCCGGGTTTAACCTCGGGCGGAACACAGCTGCGcccgaacacagataa                                                                                                                                                                                      >  2:345799/1‑235 (MQ=255)                                                              gATGTTAAGCCAGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAACGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGGTGCCTGCCACAACGTTCTGAACCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCAGCGCCAGAccccggattaccaa                                                                                                                                                                                  >  2:41810/1‑237 (MQ=255)                                                                 gTTCAGCCAGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTTCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCTCCAGACTCAGGTTATCGAACCTg                                                                                                                                                                              <  2:271497/251‑1 (MQ=255)                                                                         aGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCgg                                                                                                                                                                                                                                 >  2:122317/1‑192 (MQ=255)                                                                         aGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCgg                                                                                                                                                                                                                                 <  1:122317/192‑1 (MQ=255)                                                                             ccG‑CACCCGCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTg                                                                                                                                                                  <  2:97495/250‑1 (MQ=255)                                                                                gcccccccccccGTGGTGCTGCCGAAGCGCAATACTGCGAACTGGAGGGGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGcc                                                                                                                                                                <  2:231727/244‑1 (MQ=255)                                                                                gccaccTCCCCCGGGGTTCTGCCGAAGCTCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTGTACGCTGACGGGATGCCTGCCACAACGGTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATGCTGAAATATTAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACGCAGGTGAAGGAACTTGGAACGGTAGTTGcc                                                                                                                                                                <  2:216599/247‑1 (MQ=255)                                                                                 cACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAAACTGGGACGGTAGTTGccc                                                                                                                                                               >  1:177920/1‑250 (MQ=255)                                                                                      tCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGACAACTGCTGGGCCCAACTCAGGTTAACGAACCCGGAACGGTAGTTGcccggaa                                                                                                                                                           >  2:400934/1‑245 (MQ=255)                                                                                                 tCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGGGCGGAACACTGCTGCGCCAGACTCAGGTTAAAGAACCTGGAAAGGTAGTTGCCAGTAAACAGGCGGcccc                                                                                                                                               >  1:411640/1‑251 (MQ=255)                                                                                                      cGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCCGGAGCGGAACACTGCTGCGCCAGACTCCGGGTAACGAACCTGGGACGGTAGTTGCCCGGAAACAGTCGGCCCCGccgc                                                                                                                                          >  2:164434/1‑248 (MQ=255)                                                                                                         aGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGCACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACg                                                                                                                                       >  1:262478/1‑251 (MQ=255)                                                                                                              aaTACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGGGCGGGACACTGCTGCGCCAGACTCAGGTTCACGAAACCGGGACGGGAGTTTCCAGGAAACAGTCTGCCCCGCTGCACggcctg                                                                                                                                  >  2:288268/1‑246 (MQ=255)                                                                                                               aTACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGCGTTAAAACAGGCGGGTATCACCTGGGGCGGACAACCGCTGCGCCCGACTCAGGGTAAGCAACCTGGGCCGGTAGTTTCCCGGAAACAGTCGGCCCCCCTGcccggtatgc                                                                                                                                 >  2:180377/1‑242 (MQ=255)                                                                                                                    gCGAACTGGAGGTGGTGCCAGGCGCCCTGAACTGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGGTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAg                                                                                                                            <  2:179536/251‑1 (MQ=255)                                                                                                                         cTGGATGTGGTGCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACCCAGGATAACGCACCGGGAACGGTAGTTGCCCGTAAAAAGTCGGCCACGCCGCACGATCTGCTTTAGATTAt                                                                                                                       >  2:381644/1‑251 (MQ=255)                                                                                                                            gATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCCGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCCCCGCTGCACGATCTGCTTAAGATTAttct                                                                                                                    >  1:399831/1‑248 (MQ=255)                                                                                                                                     ccAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATGCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGATCCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCGGCACGATCTGCTTAAGATTATGCTGAAAAAGTc                                                                                                           <  2:414347/251‑1 (MQ=255)                                                                                                                                            aCCTGCACCGCTTGACGCTTACGGGAGGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGacaac                                                                                                    <  2:262478/251‑1 (MQ=255)                                                                                                                                                    cGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCCGTAAACAGTCGGCCCCGCCGCACGATCTGCTTAAGATTATGCTGGAAAAGTCGGACAACATGATCGc                                                                                            >  1:97408/1‑251 (MQ=255)                                                                                                                                                              tGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGACTCAGGTTTAAGAACCCGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGATTAAGATTATGCTGCAAAAGTCGGGCAACATGATCGCCGATACGGtt                                                                                  >  2:268236/1‑251 (MQ=255)                                                                                                                                                                 cGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCGGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTc                                                                               <  1:268236/251‑1 (MQ=255)                                                                                                                                                                 cGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTAGGACAACATGATCGCCGATAAGGTTTTc                                                                               >  2:353913/1‑251 (MQ=255)                                                                                                                                                                  gggATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTcc                                                                              <  1:128736/251‑1 (MQ=255)                                                                                                                                                                  gggAGGCCTGCCACAACGTTCGGAGCCGCGCCCGTAGGCTTTTGCCGGGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCGGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAGTTGCCAGAAGACAGTCGGCCCCGCTGCACGATCTGGTTAATATGATGGGGAAGAAGGCGTCCATCATGATCACCGATACAGTTTTca                                                                              <  2:300722/251‑2 (MQ=255)                                                                                                                                                                      ggcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCAt                                                                          <  2:353204/250‑1 (MQ=255)                                                                                                                                                                       gcctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCAt                                                                          <  1:232637/250‑1 (MQ=255)                                                                                                                                                                         ctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGGTTTCCGCCTGAt                                                                       >  1:129201/1‑251 (MQ=255)                                                                                                                                                                         ctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGAAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATAAGAGTTAAAACAGGCGGGTATCAACTGGGGAGGAACACTGCTGCGCCAGACCCAGGATAACGAAACTGGAACGGTAGGTGCCAGAAAACAGACGGCCCCGCTGCACGATATGCTTAAGGTTATGATGAAAAAGGCGGGCAACAAGATCGCCGACACCGTTTTCCGCCTGAt                                                                       >  1:291585/1‑251 (MQ=255)                                                                                                                                                                           gccACAACGTTCTGAGCCTCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCTGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTGGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATGATGCTTAAAAAGTCGGACAACATGATCGCCGTTACGGTTTTCCGCATGaga                                                                      <  1:229170/250‑3 (MQ=255)                                                                                                                                                                                                                  gCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACAATGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGGTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCCATGTGCCTGGAACATGGCGGGccc                              >  2:215763/1‑250 (MQ=255)                                                                                                                                                                                                                   cAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCGGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCgg                             <  1:41810/251‑1 (MQ=255)                                                                                                                                                                                                                    aGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCggg                            >  1:146432/1‑251 (MQ=255)                                                                                                                                                                                                                     ggaaggAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACTTGGGGCGGAACACTGCTGCTCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGTACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGt                           <  1:71030/247‑1 (MQ=255)                                                                                                                                                                                                                       aTGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGATCCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAATATTATGCTGAAAAATTCGGACAACATGATCGCCGATACGGt                                                                                   <  2:422132/193‑1 (MQ=255)                                                                                                                                                                                                                       aTGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATAAGGt                                                                                   >  1:422132/1‑193 (MQ=255)                                                                                                                                                                                                                          gAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGAc                      <  1:107976/251‑1 (MQ=255)                                                                                                                                                                                                                            gCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGt                                                                                                                                                             <  2:228291/114‑1 (MQ=255)                                                                                                                                                                                                                            gCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGt                                                                                                                                                             >  1:228291/1‑114 (MQ=255)                                                                                                                                                                                                                             ccAGCTATGCCGGTGCAATTCTGAACTCTGAGTTACAACAGGCGGGTATCCCCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTTCCAGTAAACAGTCGGCCCCGCTGCACGGTCTGCTTAAGATTATGCTGCAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACCAGGCGGGCCCGGGCGGGCGcc                   >  2:26651/1‑251 (MQ=255)                                                                                                                                                                                                                                     gCCGGTGCAAGTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCCCCAGAGTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCGGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACACCAGGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGTCCGGGTCAGACGCCGTGCGTc            <  2:385386/250‑1 (MQ=255)                                                                                                                                                                                                                                        ggTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGGTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGGCGGGCGCCGGGCGTccggt        >  1:1546/1‑247 (MQ=255)                                                                                                                                                                                                                                             aaTTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACCCTGCTGCGCCAGACTCAGGTTAACGAACCTGGCACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGGAAAAGTCGGGCAACATGATCGCCGGTACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCCGGCGTCAAGTCCTgc   >  1:382925/1‑251 (MQ=255)                                                                                                                                                                                                                                              atgcTGAAATATGAGGTAAAACAGGCGGGTAACCCCTGGAGCGGAACACGGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTATACAGTCGGCCCCGCTGCACGATCTGCTTAATATTATGCTGAAAAAGTCGGACAACTTGATCGCCGATACGGTTTTCCGCATGATAGGGCATGCGCGCTTCAATGTGCCTGGAACAGTGCGGGCCGGGTGGGACGCCGTGCGTCAGATCCTgcg  <  2:232598/248‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  ACAGTGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACG‑CACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCG  >  NC_000913/3329498‑3329984

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