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	80.8%	43.7 / 5.5	26	intergenic (+393/+249)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (8/13);  total (10/16)
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

TTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCG‑TGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTT  >  NC_000913/4295895‑4296464                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                          ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAAcgc                                                                                                                                                                                                                                                                                    <  1:135868/300‑1 (MQ=39)                   tcggctaggttgTTTATTCCGCATCCGACATCAAAGCCTGATTCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGc                                                                                                                                                                                                                                                                 <  2:712330/289‑1 (MQ=255)                      gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCa                                                                                                                                                                                                                                                               >  1:347250/1‑299 (MQ=255)                        tagggCGTTTATGCCGCGTCTGACATCAACGCCTGATGCGACGCTG‑AACGCGTCTTATCAGGTCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGAc                                                                                                                                                                                                                                                            <  2:454597/297‑1 (MQ=255)                        tAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGAc                                                                                                                                                                                                                                                            <  1:727093/300‑1 (MQ=39)                           ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTCTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAgc                                                                                                                                                                                                                                                         <  1:203263/300‑1 (MQ=255)                            gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCATGATGCGAAGTTTAACGCGTCTTAT‑CAGGCCTACGCcagccagc                                                                                                                                                                                                                                                         >  2:553043/1‑299 (MQ=255)                                  aTGCCGCATCCGACATCAACGCCTGATTCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAg                                                                                                                                                                                                                                                  <  1:296435/300‑1 (MQ=255)                                            cGACATCAACGCCTGATGCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                                                                                                                                                                         >  1:15264/1‑299 (MQ=37)                                                                aCGCTTGAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCATTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATGCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                            <  1:286464/262‑1 (MQ=255)                                                                           cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCAGCATCCGACATCAACGCCTGATGCGACGTTTACCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGTACTAGCCTGAGCTAGCGGTG‑TTTTGTAGGTCGGATAAggagctta                                                                                                                                                                                                         >  2:544138/1‑294 (MQ=11)                                                                           cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCGTCGACTGATGTGCCGTTTACCGTGTCTTAT‑CAGGCCTATGCCAGACAGCGAAATAGCATGATTTAGTG‑GGATTTTGTAGATCgggtcaggaagtta                                                                                                                                                                                                         >  2:544110/1‑288 (MQ=11)                                                                                             ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAAACAGCGCACTATCCTGACCTAGCGTAACCTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCAACATCAACGCCTGCTGCGACGTTTCACGCATCTCAT‑CAGACCTACGCCAGACAGCGAAATAGCATGACCTAGCG‑TGATTCTGTAGACCGGATAAGGCGTGTCCGCAGCCTCCGACATCaa                                                                                                                                                                                       >  1:567530/1‑300 (MQ=11)                                                                                                    gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGc                                                                                                                                                                                                                                                                 <  1:806738/220‑1 (MQ=25)                                                                                                    gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGc                                                                                                                                                                                                                                                                 >  2:806738/1‑220 (MQ=25)                                                                                                                         tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCG‑TTATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                                                           >  1:47103/1‑300 (MQ=255)                                                                                                                            ttttGTAGGTCGGATAAGGCGTTGACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCG‑TGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                                                        <  1:246798/300‑1 (MQ=255)                                                                                                                                  aggttgggtagggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATGCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                            <  1:286152/185‑1 (MQ=11)                                                                                                                                     tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCG‑TGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                                                               <  1:779689/300‑4 (MQ=9)                                                                                                                                                               gacgtaaaCGCCTGATGCGACGCTTGCCGCNTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGAc                                                                                                                                                                                                                                                            <  2:600828/160‑1 (MQ=9)                                                                                                                                                               gacgtaaaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTAT‑CAGGCCTACGCCAGAc                                                                                                                                                                                                                                                            <  2:678242/160‑1 (MQ=9)                                                                                                                                                                                                                                                           taatgCGTTTATGCCGCAGCCGATATCAACGCCTGATTCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCG‑TGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTCAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                         <  2:92305/296‑1 (MQ=11)                                                                                                                                                                                                                                                           tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATTCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCG‑TGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                         <  1:325299/300‑1 (MQ=11)                                                                                                                                                                                                                                                            aaGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCTTGTCTTAT‑CATGCCTACGCTAGACAGCGCAATAGCCTGATTTAGCG‑TGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                         <  2:48066/299‑1 (MQ=11)                                                                                                                                                                                                                                                                        ccGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCG‑TGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAACGTTGGTTAAtttt            >  2:311065/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                  cATCAACGCCTGATGCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCG‑TGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCANCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt  >  1:166482/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                          TTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTT‑AACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCG‑TGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTT  >  NC_000913/4295895‑4296464

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