Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,806,897	C→A	R296R (CGC→CGA)	proW →	glycine betaine ABC transporter membrane subunit ProW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,806,897	0	C	A	92.3%	56.6 / ‑2.3	27	R296R (CGC→CGA)	proW	glycine betaine ABC transporter membrane subunit ProW
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (8/16);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

TCTTTGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGCGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGAC  >  NC_000913/2806656‑2807138                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   ccTTTTCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTACGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCAAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGG‑TCAGATGGTACTTCGCGGTATTGgtc                                                                                                                                                                                                                                          <  2:103949/250‑1 (MQ=255)                cgATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGCGGTATCGGACGTCTGGATAT‑gggg                                                                                                                                                                                                                           <  1:343769/251‑1 (MQ=255)                         gTCTGACCATTCTGGGGATTAACCTGGTTCCGGCGGATCTGATTGAAGCCGCGCGCTCAGTCGGTGCCAGCCCGCGCCATATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCTGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑GTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCAcc                                                                                                                                                                                                                  <  2:477437/251‑1 (MQ=255)                                cATTCTGGGGAGTAACCAGGTTCCGGCGGTTCTGATTGAAGCCTCGCGCTCATTCGGTGGCAGCCCGCGCCAGAGGCTGTTCAAAGTTCAGTTGCCGCTGGCGAGGCCGACCATGATGGCGGGCGGTAACCAGACGCTGATGCTGGCCCTTGCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGGATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTggcg                                                                                                                                                                                                           <  2:8764/251‑1 (MQ=255)                                              aCCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATtgtg                                                                                                                                                                                             <  1:232565/251‑1 (MQ=255)                                               ccAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTTGCGGCGTCGGGATTGTGa                                                                                                                                                                                            <  1:187038/251‑1 (MQ=255)                                                           ggATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGGTGCTGGCCCTTTCTTTGGGGGTCATCGCCTCGATGATTGCCGTCGGCGGGGTGGGGTC‑GAAGGTACTTCGGGGGATTGGGCGGCTGGGGAT‑TGGGGTTGTCCCCCTTGGCGGCGGCGGGATTGTGGTCCttccccttt                                                                                                                                                                                >  1:130511/1‑243 (MQ=255)                                                                                  cATTCGGGGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGc                                                                                                                                                          <  2:379974/250‑1 (MQ=255)                                                                                                                 cAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAAcc                                                                                                                          <  1:161126/251‑1 (MQ=255)                                                                                                                              ccgtgggCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTacac                                                                                                             <  1:170602/246‑1 (MQ=255)                                                                                                                                gCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑GTTGGTCAGGTGGTAATTCGAGGTATCGGGCGTCTGGGTATTGGGGCT‑GCCACCCTTTGCGGGGGCCGGGTTGTGATCCTCCCCCCTATCCTCCGTCCTCTGGCGCCGGCCCTTTGGCGCGACTCCCCCCGTCGCCGCAACcgccgcccgtaccccc                                                                                                           >  1:254274/1‑241 (MQ=255)                                                                                                                                          ccGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTATACCACTGGCCCTGt                                                                                                 <  1:394128/251‑1 (MQ=255)                                                                                                                                                                         ggATGCTGGCCCTTTCTAGGGGGGTCATCGCCTCGAGGATTGCCGTCGGCGGG‑TGGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGAGCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTa                                                                   <  2:22354/249‑1 (MQ=255)                                                                                                                                                                          gATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTTGGTCAGGTGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCCCCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGCCGCAGGCCGTTTGGCGCGACCCCCGCAGTCGCGGCAACCGCCGCTGGTACACCACTGGGCCTTTTGGGCTGCTGACCCCCCCCTTCATTacataac                                                                 >  2:10390/1‑245 (MQ=255)                                                                                                                                                                           aTGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAAct                                                                <  1:376459/251‑1 (MQ=255)                                                                                                                                                                           aTGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑CTGGGTCAGATGGTACTTCGAGGTATCGGTCGGCTGGGTAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCCCCATTATCCTCGATCGTCTGACGCAAGCCGTTTGGCGCGGCTCACGCAGTCGCGGCAACCGCCGCTTGTACACCCCTTGCCCTGTTGGTCTGCTGGCCCCCCCATTCCTTAAATAAcc                                                                >  1:59548/1‑250 (MQ=255)                                                                                                                                                                             gcgtgCCCTTGCTCTGGTGGTCATCGCCTCGGTTATTGCCGTCGGCGGG‑TGGGGTCAGATGGTACTTCGAGGTACCGGTCGTCTGGATAT‑GGGGCTTGCCCCCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAActct                                                              <  2:344886/247‑1 (MQ=255)                                                                                                                                                                                           ggtgggtCATCGCCTTGATGATTGCCGTCGGCGGG‑TTGGGTCAGAGGGTCCTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCCCCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCACCCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTg                                                <  2:327883/247‑1 (MQ=255)                                                                                                                                                                                                gtCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAATTAACTCTGCACTTTCCCGGTGACGcc                                           >  1:155192/1‑251 (MQ=255)                                                                                                                                                                                                  cATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGGATCGGTCGTCTGGGTAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGGCCCCCCCAGTCGCCGCCACCCGCGCCGGGACACCCCTTGCCCCGTTGGGGTGCTGAACCCCCCCTTCATTAAGAAACCCTTCCCTTGCCCGGTGACACCgg                                         >  2:262629/1‑251 (MQ=255)                                                                                                                                                                                                          cGATGATTGCCGGCGGCGGG‑TTGGGTCAGATGGTACGTCGAGGTAGCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATc                                 <  2:380161/251‑1 (MQ=255)                                                                                                                                                                                                                 tGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCCTTTTCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGCCTCACGCAGTCGCGGCAACCGTCCCTGGTACCCCACCGGCCCTGTTGGTCTGCTGACCCGCCCATTCCTTTTTTAACTCTGCACCTGCCCGGTGACGCCGGGCCTTTTCAACCTGc                          >  2:377014/1‑251 (MQ=255)                                                                                                                                                                                                                 tGCCGGCGGCGGG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTg                           <  1:423829/250‑1 (MQ=255)                                                                                                                                                                                                                          cggG‑TTGGGTCAGATGGTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaaaa                   >  1:374416/1‑249 (MQ=255)                                                                                                                                                                                                                            gg‑TTGGGTCAGGTGGTACTTCGAGGTATCGGTCGTCGGGTTAT‑GGGGCTTGCCGCCGTTGGCGGCGTCGGGATTGTGACCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCTCCCTGCCAAAAAAAGGa               <  2:203143/251‑1 (MQ=255)                                                                                                                                                                                                                                      aGATGGTCGTTCGAGGTCTCGGTCGGTTGGATAT‑GGGGCTGGCCACCGTGGGGGGCGTCGGGATTGTGAGCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATg      <  2:248871/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggTACTTCGAGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGAc  >  2:160500/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   TCTTTGCTCTGCCGCCGATTATCCGTCTGACCATTCTGGGGATTAACCAGGTTCCGGCGGATCTGATTGAAGCCTCGCGCTCATTCGGTGCCAGCCCGCGCCAGATGCTGTTCAAAGTTCAGTTACCGCTGGCGATGCCGACCATTATGGCGGGCGTTAACCAGACGCTGATGCTGGCCCTTTCTATGGTGGTCATCGCCTCGATGATTGCCGTCGGCGGG‑TTGGGTCAGATGGTACTTCGCGGTATCGGTCGTCTGGATAT‑GGGGCTTGCCACCGTTGGCGGCGTCGGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGAC  >  NC_000913/2806656‑2807138

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