Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	522,844	G→A	intergenic (+14/‑417)	ybbP → / → rhsD	putative ABC transporter membrane subunit YbbP/protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	522,844	0	G	A	100.0%	33.6 / NA	13	intergenic (+14/‑417)	ybbP/rhsD	putative ABC transporter membrane subunit YbbP/protein RhsD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/5);  total (8/5)

GTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAGGGGCGTTCTGGTTGAATGGAACGCCTTGTTTAAATTAATTCTACAATATGTTTCATTGATGATGTTATTGAATTGGTGATTTCCTATCTTTCTATTGCTGATATTAATCTGAATCGGTGGAGTTTATTGATTTTAAATATAGCCCTACAAAAATAACTCAGACTAATAAAATCATAAATCATATGCGTTGAATGGAT  >  NC_000913/522604‑523041                                                                                                                                                                                                                                                 |                                                                                                                                                                                                      gtgAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTAGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAAGGGCGTTcc                                                                                                                                                                                              <  1:1457/250‑2 (MQ=255)                                          gCAATGGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAAGGGCGTTCTGGTTGAAtn                                                                                                                                                                                     <  1:48121/218‑2 (MQ=255)                                                     ngAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAAGGGCGTTCTGGTTGAATGGAACGCCTTGTTTAAATTAATTCTACAATATGTTTCATtgatg                                                                                                                                          >  1:32971/2‑250 (MQ=255)                                                               cTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAAGGGCGTTCTGGTTGAATGGAACGCCTTGTTTAAATTAATTCTACAATATGTTTCATTGATGATGTTATTGa                                                                                                                                >  1:149829/1‑250 (MQ=255)                                                                             ngCGAAAGCGTTTGATTTTCAGTGGGAGCCAGACTGGCGATTGTGGATTGTTATGCCGTGCAGCGGAGCGCTGCTGATGTCGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGGGATGAGAGTAAATAAAGGGCGTTCTGGTTGAATgg                                                                                                                                                                                    >  1:48320/2‑184 (MQ=255)                                                                                                        gCCAGACTGGCGATTGTGGATTGTTCGGCCGTGCAGCGGAGCGGTGCGGCTGTCGCTTTTCGTCGGCTGGCGGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTGCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAAGGGCGTTCTGGTGGAATg                                                                                                                                                                                     <  1:90524/156‑1 (MQ=255)                                                                                                                              gTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAaggg                                                                                                                                                                                                    <  1:9187/119‑1 (MQ=255)                                                                                                                                        gCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAAGGGCGTTCTGGTTGAATGGAACGCCTTGTTTAAATTAATTCTACAATATGTTTCATTGATGATGTTATTGAATTGGTGATTTCCTATCTTTCTATTGCTGATATTAATCTGAATCGGTGGAGTTTATTGATTTTAAATATAGccc                                                      >  1:58939/1‑251 (MQ=255)                                                                                                                                                             cGCTTTTCGGCGGCTGGCTGTGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAAGGGCGTTCTGGTTGAATGGAACGCCTTGTTTAAATTAATTCTACAATATGTTTCATTGATGATGTTATTGAATTTGTGATTTCCTATCGTTCTAGTGCTGATAGTATTCTGAATCGTTGGAGTTTATTGATTTTAAATATAGCCCTACAAAAATAACTCAGACta                                  <  1:157577/250‑1 (MQ=255)                                                                                                                                                             cGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAAGGGCGTTCTGGTTGAATGGAACGCCTTGTTTAAATTAATTCTACAATATGTTTCATTGATGATGTTATTGAATTGGTGATTTCCTATCTTTCTATTGCTGATATTAATCTGAATCGGTGGAGTTTATTGATTTTAAATATAGCCCTACAAAAATAACTCAGACta                                  >  1:2964/1‑250 (MQ=255)                                                                                                                                                                         gctggctgGGTGCGCGACTGGTTAAGGGTAAGTCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAATGGCGTTCTGGTTGAATGGAACGCCTTGTGTAAATTAATTCTACAATATTTTTCATTGACGATGTTATTGAATTGGTGATTGCCTATCTTTCTATTACTTATATTAATCTGAATCGGTGGAGTTTATTGATTTTAAAtatan                                                         <  1:14114/215‑2 (MQ=255)                                                                                                                                                                                            ggTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTACATAAAGGGCGTTCTGGTTGAATGGACCGCCTTGTTTAAATTAATTCTACAATATGTTTCATTGATGCTGTTATTGAATTGGTGATTTCCTATCTTTCTATTGCTGATATTAATCTGAATCGGTGGAGTTTATTGATTTTAAATATAGCCCTACAAAAATAACTCAGACTAATAAAATCATAAATCATATGCGTTGAATGGat  >  1:141052/1‑251 (MQ=255)                                                                                                                                                                                            ggTTAAGGGTAAGGCGCTGATCAGGCAGTTTGCGGGGTGATGAGAGTAAATCAAGGGCGTTCTGGTTGAATGGAACGCCTTGTTTAAATTAATTCTACAATATGTTTCATTGATGATGTTATTGAATTGGTGATTTCCTATCTTTCTATTGCTGATATTAATCTGAATCGGTGGAGTTTATTGATTTTAAATATAGCCCTACAAAAATAACTCAGACTAATAAAATCATAAATCATATGCGTTGAATGGat  >  1:138689/1‑251 (MQ=255)                                                                                                                                                                                                          ngcTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAAGGGCGTTCTGGTTGAATGGAACGCCTTGTTTAAATTAATTCTACAATATGTTTCATTGATGATGTTATTGAATTGGTGATTTCCTATCTTTCTATTGCTGATATTAATCTGAATCGGTGGAGTTTATTGATTTTAAATATAGCCCTACAAAAATAACTCAGACTAATAAAATCATAAATCATATGCGtt          >  1:40229/2‑229 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                      GTGAGTTCGCCATGCTTGGGTTTGTTTCCGGCCTGGTGGCCGCAATTGGTGCGGAAACGGCGCTGGCAGTGTTGCAGGCGAAAGTGTTTGATTTTCCGTGGGAGCCAGACTGGCGATTGTGGATTGTTCTGCCGTGCAGCGGAGCGCTGCTGCTGTCGCTTTTCGGCGGCTGGCTGGGTGCGCGACTGGTTAAGGGTAAGGCGCTGTTCAGGCAGTTTGCGGGGTGATGAGAGTAAATAAGGGGCGTTCTGGTTGAATGGAACGCCTTGTTTAAATTAATTCTACAATATGTTTCATTGATGATGTTATTGAATTGGTGATTTCCTATCTTTCTATTGCTGATATTAATCTGAATCGGTGGAGTTTATTGATTTTAAATATAGCCCTACAAAAATAACTCAGACTAATAAAATCATAAATCATATGCGTTGAATGGAT  >  NC_000913/522604‑523041

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