Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,381,378	T→C	G212G (GGT→GGC)	degQ →	periplasmic serine endoprotease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,381,378	0	T	C	95.2%	66.0 / NA	22	G212G (GGT→GGC)	degQ	periplasmic serine endoprotease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (2/18);  total (3/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTAACCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGTAACTCCGGCGGTGCACTATTAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAAC  >  NC_000913/3381254‑3381500                                                                                                                             |                                                                                                                           tcggtAACCCATTTGGCCTTGGGCGACCCGCCACCTCTGGCCTTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCcgg                                                                                                                    <  2:223151/133‑1 (MQ=255) tcggtAACCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCcgg                                                                                                                    <  1:475740/133‑1 (MQ=255)   ggtAACCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCcggcg                                                                                                                  <  2:59880/133‑1 (MQ=255)      aaCCCATTTGGCGTTGGGCAAACCGCCACATCTGGCCTTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTg                                                                                                               <  1:486674/133‑1 (MQ=255)      aaCCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTg                                                                                                               <  1:303331/133‑1 (MQ=255)      aaCCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTg                                                                                                               <  1:445755/133‑1 (MQ=255)      aaCCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTg                                                                                                               <  1:484506/133‑1 (MQ=255)      aaCCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTg                                                                                                               <  1:68608/133‑1 (MQ=255)       aCCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGc                                                                                                              <  2:399434/133‑1 (MQ=255)       aCCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGc                                                                                                              <  2:277397/133‑1 (MQ=255)       aCCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGc                                                                                                              <  2:414775/133‑1 (MQ=255)               ggCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGCACTATTaa                                                                                                      <  2:378001/133‑1 (MQ=255)                 ccTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGCACTATTAAAc                                                                                                    <  2:335449/133‑1 (MQ=255)                     gggCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTACATTAACCCCGGGAACTCCGGCGGTTCACTGTTAAACCTTa                                                                                                >  2:130467/1‑133 (MQ=255)                      ggCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGCACTATTAAACCTTaa                                                                                               <  2:472305/133‑1 (MQ=255)                                 aCCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGCACTATTAAACCTTAACGGTGAGTTaa                                                                                    <  2:127220/133‑1 (MQ=255)                                   ctctGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGCACTATTAAACCTTAACGGTGAGTTAAtt                                                                                  <  2:86747/133‑1 (MQ=255)                                        gCATTGTTTCCGCATCAGGCCGCAACGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGCACTATTAAACCTTAACGGTGAGTTAATTGGCAt                                                                             <  2:197623/133‑1 (MQ=255)                                                  cGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGCACTATTAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCa                                                                   >  1:67179/1‑133 (MQ=255)                                                        aGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGCAACTCCGGCGGTGCACTATTAAACCTTAACGGGGAGTTAAATGGGATGAACACTTCAATCCtt                                                             >  1:146243/1‑133 (MQ=255)                                                                                                 tCCAGACAGATGCTTCCATTAACCACGGTAACCCCGGCGGTACACTATTAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATccc                    <  2:146243/133‑1 (MQ=255)                                                                                                                   ttAACCGCGGCAACTCCGGCGGTGCACTATTAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAac  <  2:67179/133‑1 (MQ=255)                                                                                                                             |                                                                                                                           TCGGTAACCCATTTGGCCTTGGGCAAACCGCCACCTCTGGCATTGTTTCCGCATTAGGCCGCAGCGGGTTGAATCTTGAAGGTCTGGAAAACTTTATCCAGACAGATGCTTCCATTAACCGCGGTAACTCCGGCGGTGCACTATTAAACCTTAACGGTGAGTTAATTGGCATCAACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAAC  >  NC_000913/3381254‑3381500

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑