Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,300 C→T intergenic (+506/+136) 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_0009134,296,3000CT85.0% 30.6 / 0.5 20intergenic (+506/+136)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (7/10);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

GACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTC  >  NC_000913/4296052‑4296534
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
gACATCAACGCCTGATGCTACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                                                                                                                                                                                                            <  1:124041/251‑1 (MQ=25)
                                                                                                                                                             aCGCCAGACAGTGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTGAACCGTTGGTTAAttttt                                                                                 <  1:135999/249‑1 (MQ=255)
                                                                                                                                                                ccAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATAAGGCGTTTATGCCGCATTCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                                            >  1:235877/1‑251 (MQ=255)
                                                                                                                                                                  aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGATATCAACGCCTGATGAGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATTCGACATCAATGCCTGATGCGTCGCTTTCCGCGTATTATAAGGCCTATCTTAACCGTTGTTTAATTATTCAAActc                                                                          >  1:453509/1‑251 (MQ=11)
                                                                                                                                                                  aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATACGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATAAGGACTACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTTACGCCGAATCAGACATCAATGCATGATGCGACGCTTGCCGTGTCTTATCAGGCATATATTAACCGTTGGTAAAGTTTTCAAActc                                                                          >  1:453504/1‑251 (MQ=11)
                                                                                                                                                                                ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                                            >  1:280187/1‑251 (MQ=255)
                                                                                                                                                                                   gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                                          >  1:101985/1‑250 (MQ=255)
                                                                                                                                                                                          cggacgtttctaggctggATAAGGCGTTTACGCCGCACCCGACATCAACGCCTGATGCGACGTTTAACGCGCCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGATTTGGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                                  <  2:14837/235‑1 (MQ=17)
                                                                                                                                                                                                  tGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAg                                           <  1:110419/250‑1 (MQ=255)
                                                                                                                                                                                                  tGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAGCAAACCCTTCCATACAg                                           <  1:489522/250‑1 (MQ=255)
                                                                                                                                                                                                         cGGATAAGGCGTTTCTGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTCGGTCGGATACGGCGTTTACGCCTCATCCGACATCAATGCCTGAAGCGACGCTTCCCACGTCTTATCAGGCCTATCTTAACCGTTCTTTATTTTTTCACCCTCTTCCCACACGGTTTCAAAACCTTCCATACACATCAGGtt                                   >  2:43044/1‑251 (MQ=11)
                                                                                                                                                                                                                 gCGTTTAGTCCGCTTCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt                           <  2:193089/251‑1 (MQ=18)
                                                                                                                                                                                                                  cGTTTATGCCACATCCGACATAAACACCTGATTCGACGTTTTACGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTg                          <  2:110923/251‑1 (MQ=18)
                                                                                                                                                                                                                       aTGCCGCATCCGACATAAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGAGCAgcg                     >  1:464912/1‑251 (MQ=255)
                                                                                                                                                                                                                             cATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtttt               <  1:285638/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGATTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATAAGGCATATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCAGGTATCACACCCTTCCATACAGCTCAGGttaacccaattcaccgccctttact            >  2:273922/1‑228 (MQ=17)
                                                                                                                                                                                                                                    cATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGtt        <  1:411174/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tCAACTCCTGATGCTACGTTTAACTCGTCTTATCAGTCCTACGCCAGACATCGCAATATCCTGATTTAGCGTGATTTTGTAGGACGGATAAGGCGTTTACGCCGCAGCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttc       <  2:487594/250‑1 (MQ=21)
                                                                                                                                                                                                                                       caacgccgtgtGCGACGCTTAACGCGTCTTACAATGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGGATGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCTAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttctt     <  2:453093/241‑1 (MQ=21)
                                                                                                                                                                                                                                            ccTGATGCAACGTTTAACGCGCTTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTGAGCGTGATTTTGTAGGTCCGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc  <  2:110865/249‑1 (MQ=21)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            
GACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTC  >  NC_000913/4296052‑4296534

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