Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,143,377	C→T	S283L (TCA→TTA)	frwC →	putative PTS enzyme IIC component FrwC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,143,377	0	C	T	95.5%	69.7 / NA	22	S283L (TCA→TTA)	frwC	putative PTS enzyme IIC component FrwC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/12);  minor base G (0/1);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAATTAATTGAGGAACG  >  NC_000913/4143135‑4143619                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   cTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTtgtg                                                                                                                                                                                                                                            >  1:285683/1‑251 (MQ=255)       gCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGAGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCa                                                                                                                                                                                                                                       <  1:186188/250‑1 (MQ=255)               gCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTgcg                                                                                                                                                                                                                              <  1:309199/251‑1 (MQ=255)                     cTCAGGGTGTTTATACCGTGGGGGCCATCGCTGCCGTGGGGATTTGTATCCCGCCGCTGGGGCTGGGGCTGGCGACGCTGATTGGTCGTAAAACTTTCTCCGCAGCAGAGCGCGAAACCGGCAAAGCGGCAATGGAGATGGGCTGTGTTGGCGT‑TCCTGAGAGGGCGATTCCTTTCGCCGCTGCCGAAACGCTGCGTGTCATTCCTTCCATCATGGCCGGTTTGGTTTGTGGTGCAGGAACTGCGGCGCTg                                                                                                                                                                                                                        >  1:322659/1‑251 (MQ=255)                              tttATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTgcg                                                                                                                                                                                                               >  1:8401/1‑251 (MQ=255)                                   cccggggtggtCATCGCTGCCGTGGTTATTTGTTTCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGTCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGt                                                                                                                                                                                                           <  2:228944/245‑1 (MQ=255)                                                                 tGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGtgct                                                                                                                                                                            >  1:98371/1‑251 (MQ=255)                                                                             cTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTg                                                                                                                                                                 >  1:76549/1‑250 (MQ=255)                                                                                              gactcttaTTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCCCTGGTGATGGGCTTTGTTGGCGTGTA‑GGAAGGGGCGATTCCGTTCGCCGCTGTCGATCCGCTGCGTTTCATTCCTTCCAGCATGGTCGGTTTAGCTTGTGGTGCAGTACCTGCTGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGAGTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTtata                                                                                                                                               <  1:18719/244‑1 (MQ=255)                                                                                                cGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATc                                                                                                                                             <  1:162055/251‑1 (MQ=255)                                                                                                          tCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCGGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTg                                                                                                                                    <  2:285683/250‑1 (MQ=255)                                                                                                              aaaaaTTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTg                                                                                                                                <  2:98371/250‑1 (MQ=255)                                                                                                                                       aaaCCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGAGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGGGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGGCGGAGCAGTGGTGACGGCTGGTtgg                                                                                                      >  2:211990/1‑250 (MQ=255)                                                                                                                                               aaaGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTAGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGACTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGGGGTGACGGCAGTTtgtgggtacgg                                                                                              >  2:230452/1‑244 (MQ=255)                                                                                                                                                              aTGGGCTGTGTTGGCGG‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGAGCCGCTTCGTGTCATTCCTTCCATCATGGTCGGTTGAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCGGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGgc                                                                               <  2:8401/251‑1 (MQ=255)                                                                                                                                                                gggCTGTGTTGGCGG‑TACTGAAGGGGCGATTCCTTTCGCCGATGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGTTGTTTGTGTGAACGTGCTGAAAAGTCTGgcg                                                                              <  1:282151/250‑1 (MQ=255)                                                                                                                                                                        ttGGCGA‑TACTGAAGGGGTGAATACTTTCGCCGCTGCCGCTCCGCTGCGTGTCATTCATTCCATCATGGTAGGTTTAGTTTGAGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGTAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGGGCAGTGGTGACGGCTGTTTTTGTGAACATGCTGAAAAGTCACGCGCGTAAAAAt                                                                     >  2:224628/1‑251 (MQ=255)                                                                                                                                                                                                   ttCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAaagaa                                             <  1:211990/249‑1 (MQ=255)                                                                                                                                                                                                     cgccgcTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGGCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAgac                                          <  2:325150/250‑1 (MQ=255)                                                                                                                                                                                                     cgccgcTGCCGATCCGCTGCGTGTCATGCCTTCCATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGGTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAgacg                                         <  2:109808/251‑1 (MQ=255)                                                                                                                                                                                                                                 ccttccATCATGGTCGGTTTAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAaattaatt              >  1:91298/1‑250 (MQ=255)                                                                                                                                                                                                                                            ggTCGGGTTAGTTTGTGGTGCAGGAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAATTAATTGAGGAACg  <  2:69263/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   CTTCATGCTGATTTGCGTTGCTCAGGGTGTTTATACCGTGGTGGCCATCGCTGCCGTGGGTATTTGTATCCCGCCGCTGGGGATGGGGCTGGCGACGCTGATTGGTCGTAAAAATTTCTCCGCAGAAGAGCGCGAAACCGGCAAAGCGGCACTGGTGATGGGCTGTGTTGGCGT‑TACTGAAGGGGCGATTCCTTTCGCCGCTGCCGATCCGCTGCGTGTCATTCCTTCCATCATGGTCGGTTCAGTTTGTGGTGCAGTAACTGCGGCGCTGGTCGGTGCGCAGTGCTATGCAGGCTGGGGTGGTCTGATTGTGCTGCCGGTGGTTGAAGGCAAGCTGGGTTATATCGCAGCAGTGGCTGTCGGAGCAGTGGTGACGGCTGTTTGTGTGAACGTGCTGAAAAGTCTGGCGCGTAAAAATGGGTCTTCGACTGATGAAAAAGAAGACGACCTGGATTTGGATTTTGAAATTAATTAATTGAGGAACG  >  NC_000913/4143135‑4143619

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