Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,296,187	C→T	intergenic (+393/+249)	gltP → / ← yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,187	0	C	T	95.7%	54.1 / NA	23	intergenic (+393/+249)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/12);  minor base A (1/0);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTT  >  NC_000913/4295944‑4296302                                                                                                                                                                                                                                                     |                                                                                                                    caccccctgttgataCGCTTAACGGGTCTTATCAGGCCTACGCCAGAAAGCGCAATAGCCTTTTTTAGCGTGATGGTGTAGGCCGGATAAGCCGTTTTCGCCGCAGC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAAc                                                                                                                <  2:61659/235‑1 (MQ=11) cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCGACGTCAGACGGCGTAATGGCGTGATTTAGCGTGATTTTTTAGGCCGGATAAGGCGTTTACGCCGCATC‑GTATATCACCGCCTGATGCGACGCTTGCCGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGACGCATCCGACATCAACGCCTGATGCGACGTTTAAcg                                                                                                               <  1:323364/250‑1 (MQ=14)            gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAg                                                                                                    <  1:549613/250‑1 (MQ=39)               aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCt                                                                                                <  1:187764/251‑1 (MQ=255)                         cgTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGATTTAGCGTGATTCTGTACGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCAAGGCAGCGCCATAGCATGCTTAACCGTGGTTTTGTCGGTCGGATGCGGCGTTTAGGCCGCAGCCGACATCAACTCCTTATGCGAGGATTACCGCGTCTTATCAGTCCTACACCacgc                                                                                       >  2:287110/1‑247 (MQ=11)                                           cctgtatgctcAATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATATGCCGTTTACTCCGC‑TCTCGACAGCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTGAGCTTGATTTTGTAGTTCGGATAAGGCGTTTATGCCGCACCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                    <  2:131383/241‑1 (MQ=11)                                             aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTAc                                                                                              <  1:173157/223‑1 (MQ=25)                                                   cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCa                                                                                          <  1:153015/221‑1 (MQ=25)                                                    gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGAttt                                                           >  1:193089/1‑251 (MQ=21)                                                        tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                        >  1:567531/1‑250 (MQ=25)                                                         aGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCTCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                      <  1:565582/251‑1 (MQ=255)                                                           ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTATGCCGCATCCGACATACACCCCTGCTGCGACGTTTAACGCGTCTTCTCCGGCCTACGCCAGCCATCCCAATAGCCTTATTTaccgattggaggagcgt                                                    >  2:38160/1‑235 (MQ=11)                                                           ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg                                                                   <  1:38160/236‑1 (MQ=25)                                                              gATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCg                                                  <  1:32854/250‑1 (MQ=25)                                                                tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAt                                               >  1:476038/1‑251 (MQ=21)                                                                      gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTAc                                                                                              >  1:480922/1‑198 (MQ=37)                                                                            ttGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAATGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAAAGCCTGATTTCGCGTTATTCTGTAGATCAGATACGGCGtata                                     >  2:117864/1‑246 (MQ=11)                                                                                     ggATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTAc                                                                                              >  1:201615/1‑183 (MQ=37)                                                                                          aGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc                     >  1:83222/1‑251 (MQ=255)                                                                                                     ccgcATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg          <  1:412320/251‑1 (MQ=18)                                                                                                       gcATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTCGGTCGAATAAGGCGTTTACGCCGCATCCGACATAAATGCCTGATGCg        >  2:484383/1‑251 (MQ=11)                                                                                                         aTC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATACGACATCAATGCCTGATGCGAc      >  1:14837/1‑251 (MQ=255)                                                                                                              gACATCAACGCCTGATGCTACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCtt  <  1:124041/251‑1 (MQ=25)                                                                                                                                                                                                                                                     |                                                                                                                    CAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATC‑CGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTT  >  NC_000913/4295944‑4296302

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: ‑