Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,296,381	+GC	intergenic (+587/+55)	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,380	1	.	C	96.4%	79.1 / ‑3.8	28	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (14/13);  minor base G (1/0);  total (15/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.76e-01
*	NC_000913	4,296,380	2	.	G	96.4%	82.5 / ‑4.1	28	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (14/13);  total (15/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.76e-01

GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGC  >  NC_000913/4296193‑4296611                                                                                                                                                                                              ||                                                                                                                                                                                                                                         gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACTCTTGCCGCGTC‑TTATCAGGCCTATCTTAAcc                                                                                                                                                                            <  1:217675/251‑1 (MQ=2)          cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTa                                                                                                                                                                    <  1:69544/250‑1 (MQ=19)             gcccaCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAAAAGGCGTTTACGCCCCATCCGACATCAACGCCTGATGAGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                                                                                                                <  2:345906/247‑1 (MQ=11)                         cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                                                                                                                    >  1:91590/1‑251 (MQ=255)                             gcAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCGGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTGCAAACTCTTCacac                                                                                                                                                <  2:431500/251‑1 (MQ=14)                              cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                                                               >  1:341272/1‑251 (MQ=255)                                aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGTCTGATTACGCTTGATTTTCGAGGTCGGATCAGGCG‑TTCACGCCGCATCCCACATCAATCCCTGCTTCTTCCCTTGCCGCGTC‑TTATCAACCCTATCTTCACCGTTAGTTAATCTTTCAATCCCTTCACACCCg                                                                                                                                             >  2:177871/1‑251 (MQ=2)                                aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                                                                                                             >  1:170606/1‑251 (MQ=255)                                      tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                                                                                                                        >  1:9336/1‑250 (MQ=255)                                       gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCTTGGTATATCCTGAATTCGCTTGATTCTGTATCTCGGATACGGCGATACAG‑CCGCATCCGACATCGATGCCTGATTCGTCGCTGGCCTCGTCGTTA‑CAGGCCTCTCCTCACCTTAGGTTCATTTTTCCCACCCTTCACACCGAGTATCa                                                                                                                                       >  2:37273/1‑250 (MQ=2)                                           tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc                                                                                                                                   <  1:364766/250‑1 (MQ=255)                                                    tttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                          <  1:356382/250‑1 (MQ=255)                                                     ttGTAGGTCGGATAAGACGTTTATGCCGCATCCGACATCAACGCATGATGCGGCGCTTAACGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGTATCCGACATCAATGCCTGATGCGACGCATGCCGCGTC‑TTCTCAGTCCTATCTTAACCGTTGGTTAATTTTTCAACCTCTTCACCCCCGGTATCAAACCCTTCCATACa                                                                                                                         >  2:24306/1‑250 (MQ=17)                                                        tAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                                                                                                                     >  2:127048/1‑251 (MQ=255)                                                             cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGACGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGtt                                                                                                                <  1:30768/251‑1 (MQ=255)                                                                     gCGTTTATGCCGCATCCGTCATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt                                                                                                        <  2:288300/251‑1 (MQ=255)                                                                       gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc                                                                                                      >  1:314860/1‑251 (MQ=255)                                                                         ttATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc                                                                                                      >  1:225206/1‑249 (MQ=255)                                                                                  aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTg                                                                                           <  2:151345/251‑1 (MQ=255)                                                                                  aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTg                                                                                           >  1:170157/1‑251 (MQ=255)                                                                                           cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttctt                                                                                  >  2:163859/1‑251 (MQ=255)                                                                                             aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTtcttct                                                                                <  1:180106/251‑1 (MQ=255)                                                                                                                         aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGa                                                     <  1:127048/250‑1 (MQ=255)                                                                                                                               ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAGCCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatcat                                              <  2:225206/251‑1 (MQ=255)                                                                                                                                  aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAtt                                             <  1:16087/249‑1 (MQ=255)                                                                                                                                       aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGccc                                      >  1:447690/1‑251 (MQ=255)                                                                                                                                                                           ggTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGCTACAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGACTAACCTGTTCGGGAAATTGCAGAGc  >  2:296491/1‑251 (MQ=38)                                                                                                                                                                           ggTCGGATAAGGCG‑TTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGc  >  1:130036/1‑251 (MQ=255)                                                                                                                                                                                              ||                                                                                                                                                                                                                                         GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTC‑TTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGC  >  NC_000913/4296193‑4296611

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