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

TGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACG‑GTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCC  >  NC_000913/3329502‑3329986                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      tGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGGGCAATTCTGAAATATGAGTTa                                                                                                                                                                                                                                               >  2:30841/1‑250 (MQ=255)        aaGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑c                                                                                                                                                                                                                                      >  1:189657/1‑251 (MQ=255)         agcctgaggATATGACTTTTATACGCGTGGCATCTTATTACCCCGTTCCGATGTTGAGGCAGGTACGCACCCTCCCCCGTGGTTCTGACGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGATGAAGGGATGCCTGCCACAACGTTCTGAGCCGCTCCAGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑Cg                                                                                                                                                                                                                                     >  2:47606/9‑251 (MQ=255)               gTGATATGGCTTTTATACGCGTGGCATCTTATTACCCAGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTCTGCCGGTGCAATTCTGATATATGAGTTAAAACAGG‑CGGGTATc                                                                                                                                                                                                                               >  2:392123/1‑251 (MQ=255)                tGATATGGCTTTTATACACGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGCCCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCCGCTATGCCGGTGCAATTCTGAAATATGAGTTATA‑CAGGGCGGGTATcc                                                                                                                                                                                                                              >  2:171938/1‑250 (MQ=255)                             ataCGCGTTGCATCTTAGTCCCCCGTTACGTTGTTCAGCCATGTACGCACCCTCCCCCGTGGTTCGGCGGAAGCGCAATACGGCGAACTGGATGTGGGGCCAGGCGACCTGAGCCGCTTGACGCTGACGTGATGCCTTCCACAACGTGCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGTATGTAGCCAGGTATGCCGGTGCAATTCTGAAATATAAGTGAAAACACG‑CGGGTATGTCCTGGAGCGGAac                                                                                                                                                                                                                 <  2:367970/251‑1 (MQ=255)                             ataCGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGTTACGCACCCTCCCCCGTGGTTCTGCCGACGCGCAATACTGCGAACTGGATGTGGTGCCCGGCGAGCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTTCCATGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAATCAGG‑CGGGTATCACCTGGAGCGGaa                                                                                                                                                                                                                  <  1:158215/250‑1 (MQ=255)                                       cATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTTGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGCTGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGCGCGGGACACAGCTgcgc                                                                                                                                                                                                       >  2:106192/1‑251 (MQ=255)                                             agtcccacGTTACGATGTTCAGCCAGGTACGCACCCTCCCCTGTGTTTCTGCCGAAGCGTATTACTGCGATCTGGATGTGGTGCCATCCGACCTTAACCGCTTTACGTTGACTGGATGCCTGCCACAACGTTCTGAGCCTCTCCCGTTGGCTTTTGCCGTGCAGGATGGCGCCAGCGCTGCCGGTGCAATTCTGAAATATGAGCTATAACAGG‑CGGGTACCTTCTGGATCGGAACACTGCTGCGCCAGAcg                                                                                                                                                                                                 <  2:353707/244‑2 (MQ=255)                                                acaccGTTACGATGTTCAGACAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCCGGGGAGGAAACCTGCTGCGCCAGACTccg                                                                                                                                                                                              >  2:284822/4‑249 (MQ=255)                                                 ccccGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTATAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCCGACTCAgg                                                                                                                                                                                             >  1:172879/1‑251 (MQ=255)                                                          gATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTTCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGACTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGATAAAACAGG‑CGGGGATCCCCCGGCGCGGAACCATGCTGCGCCAGACTCAGGTTaaagacc                                                                                                                                                                                    >  2:53819/1‑246 (MQ=255)                                                          gATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGGATCACCTGGAGCGGAACAATGCTGCGCCAGACTCAGGTTAACGAAc                                                                                                                                                                                    >  1:190608/1‑251 (MQ=255)                                                          gATGTTCAGCCAGGTACGCACCCCCCCCCGTGGTCCTGCCGAAGCGCAAGACTGCGAACTTGATGGGTTGCCAGGCGTCCTGGGCGGCTTTACGCTGACGGGATGCCTGCCACGACGTGCTGAGCCGTTCCCGGTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGGCGTTGCAATTCTGATATATGAGTTAAAACAGG‑CGGGTAGCACCTGGAGTGTAACACTGCTGCGCCAGACTCAGGTTAACGAAc                                                                                                                                                                                    <  1:242547/251‑1 (MQ=255)                                                            tGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGATGCGACAGACTCAGGTTAAAGAACCt                                                                                                                                                                                  >  1:216643/1‑251 (MQ=255)                                                                         aCGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTgggaagggagttg                                                                                                                                                                     >  1:331717/1‑240 (MQ=255)                                                                            cACCCTCCCCCGTGGTTCTGCCGAAGCGAAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCa                                                                                                                                                                  <  2:331717/251‑1 (MQ=255)                                                                                                       gcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCCGTAAACAGTCGGCCCCGCTGccagctc                                                                                                                                       >  1:301925/1‑245 (MQ=255)                                                                                                       gcAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATc                                                                                                                                       >  1:123320/1‑251 (MQ=255)                                                                                                                  aaCTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGAtt                                                                                                                            <  1:30841/251‑1 (MQ=255)                                                                                                                         tgtgGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGCGTTAAAACAGG‑CGGGTATCACCTGGGGCGGGACACTGCTGCGCCAGACTCCGGTTAACGAACCTGGAACGGTAGTTGCCCGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTg                                                                                                                      >  1:58810/1‑250 (MQ=255)                                                                                                                          gtggtgCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGaa                                                                                                                    <  1:101751/251‑1 (MQ=255)                                                                                                                          gtggagCCAGGCTACCGTAACCGCTTTACGCTGACTGGATGCCTGCCACAACGTTCCGAGCCGCTTCCGTTGGCTTTTGCCGTGCAGGATGGATCCAGCTATTCCGGTGCAATTCTGAAATATTAGTTAAAACAGG‑CGGGTATCACCTGTAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGTCCCGCTGCACGATCTGCTTAAGATTATGCTGaa                                                                                                                    <  1:291988/251‑1 (MQ=255)                                                                                                                                      gACCTGACCCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGacaa                                                                                                        <  2:89592/251‑1 (MQ=255)                                                                                                                                        ccTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTTCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGacaaca                                                                                                      <  1:47606/251‑1 (MQ=255)                                                                                                                                                gCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGcc                                                                                              <  2:123320/251‑1 (MQ=255)                                                                                                                                                gCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGcc                                                                                              >  1:111009/1‑251 (MQ=255)                                                                                                                                                 cTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGccc                                                                                             >  1:13152/1‑250 (MQ=255)                                                                                                                                                          gACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACG‑Gttt                                                                                   <  1:320857/251‑1 (MQ=255)                                                                                                                                                          gACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGAGCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACG‑Gt                                                                                     <  2:174764/249‑1 (MQ=255)                                                                                                                                                                    ctgccACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAAAACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAATTGCCAGTAAACAGTCGGCCCCCCTGCACGCTCTGCTTAAGATTATGCTGAAAAACTCGGACAACATGATCGCCGATACG‑GTTTTCCGCATGAt                                                                         >  1:42897/1‑251 (MQ=255)                                                                                                                                                                       ccatagcGTTCTGAGGCGATCCCGTTGGCATTGGCAGTGCAGGATGGAGCCTGCTATGCCGGTGCAATTTTGAAATATGAGTTAAAACAGG‑CGGGTATCACCGGGAGCGGAACACTGCTGTGCCAGAATCAGGTTAACGAACCTGGAATGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAGCATGATCTCCGATACG‑GTTTTCCGGATTATAga                                                                      <  2:223285/245‑2 (MQ=255)                                                                                                                                                                          caACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACG‑GTTTTCCGCATGATAGGCCa                                                                   <  1:319552/251‑1 (MQ=255)                                                                                                                                                                               ttCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTCTGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGGATTACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGCTAACGAACCTGGAAAGGCAGTTGCCAGTAAACAGGCGGCCCCGCTACACGATCTGCTTATGATTATGCTGAAAAAGTCGGACAACATGATCGCCGCTACG‑GTTTTCCGCATGATAGGCCCTgcgc                                                              >  1:266685/1‑251 (MQ=255)                                                                                                                                                                                tCTGAGCAGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCACTTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGTACACAGGTGCGACAGACCAAAGTTAACGAACATGGCACGGGAGTTGCCCGTAAACAGGCGGCCCCGATGCACGACCTGGTTAAGATTATTCTTAAAAAGGAGGACAACAGAATCGCCGAACCGAGTTTTC‑GCATGATaagacaagcgcg                                                             >  2:130516/1‑240 (MQ=255)                                                                                                                                                                                                                       gCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGGGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGGAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTTTGCTGGAAAAGTCGGACAACAt                                                                                                     >  2:203485/1‑173 (MQ=255)                                                                                                                                                                                                                       gCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACAt                                                                                                     <  1:203485/173‑1 (MQ=255)                                                                                                                                                                                                                                    gTGCAATTCTGAAATATGAGTTAAAACAGG‑AGGGGATCACCTAGAGCGGAACACTGCTGCGGCAGAATCAGGTTAACGAACCTGGAACGGAAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACG‑GTTTTCCGCATGATAGGCCATGCTCGCTTCAATGTGCCTGGAAAATGGCGGGCCGGGGCGGAAGCCGTGCGTc              >  1:348942/1‑246 (MQ=255)                                                                                                                                                                                                                                       cAATTCTGAAATATGAGTTAAAACAGG‑CGGGGATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTCTGCTGAAAAAGTCGGACAACATGATCGCCGCTACG‑GTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGGCGTTCGTTAGAACCTg      >  1:270897/1‑251 (MQ=255)                                                                                                                                                                                                                                              gAAATATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACG‑GTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGcc  <  1:126309/248‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      TGCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGG‑CGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACG‑GTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCC  >  NC_000913/3329502‑3329986

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