Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,191 A→C 25.5% intergenic (+397/+245) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,1910AC25.5% 13.8 / 11.0 24intergenic (+397/+245)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (13/4);  new base C (3/3);  total (16/8)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

GGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACT  >  NC_000913/4295911‑4296461
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                 
ggttgggtaggcggTTTATGCCGCATCATACATGACCGGCTGATGCGACGCTTACCGGGTGTTATCAGGCCTACGGCAGACAAGGGAAGATCCTAATTTAGGGTGATGTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCTACGCTTGCCGTGTCTTATCAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                                                                                                                                                                                                                                                                                   <  2:784410/268‑1 (MQ=11)
           gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACTCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCTGCATCCGACATCAACACCTGATGCGACGCTTAACGCGTCTTac                                                                                                                                                                                                                                                                         >  1:801869/1‑279 (MQ=255)
                                               aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATGCCAGACAGCGGAATAACCTGATTTATCGTCATTTTTTATGCCGGATCAGCAGTTTACGCCCCATCCGACATCAATGCCTGATGCCACGCTTAACCCGTCTTATCAGGCCTACGCAAGACCTCGGAATAGTCCGATTTAgt                                                                                                                                                                                                                                    >  2:1084596/1‑280 (MQ=11)
                                                                   ggCCGACGCCATGCAGCGCGATAGCCTGCTTTAGCGTGATTTTGTGGTTCGGATAAGGCGTTTACACCTCATCCGACATCAACTCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAt                                                                                                                                                                                                                 <  1:528113/280‑1 (MQ=11)
                                                                     ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGACAAGGCGTTTACGCCGCATCCGACATCAACGCCTGCTGCGGCGCTTAACGCGTCTTATCAGCCCTACGCAAGCCCGAGAAACAGCCTGATATAGCGTGATTTCGTACGTCGCACAAg                                                                                                                                                                                                              >  2:333598/1‑281 (MQ=11)
                                                                        aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCTATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATCGCCTGAATTAGCGTGATTTTATAGGTCGGATAAGGc                                                                                                                                                                                                            >  1:471969/1‑280 (MQ=16)
                                                                         cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                                                    <  1:1015997/223‑1 (MQ=34)
                                                                         cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                                                    >  2:1015996/1‑223 (MQ=34)
                                                                              gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCTAGACAGCGCAATAGCCTGTTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATCTTGCAGGTCGGATAAGGCGTTTacg                                                                                                                                                                                                     >  1:160879/1‑279 (MQ=17)
                                                                                   cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATAAGGCGCTTATGCTGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCTAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTCGGTCGGATAACGCGTTTACGCAGc                                                                                                                                                                                                 >  2:287112/1‑280 (MQ=11)
                                                                                   cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCCTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATCAGGCGCTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGCCTTATCAGTCCTACGCCTGACAGCGCAATTGCATGATTGAGCGTGACTTTGAAGGTCGGATAAGTTGTTTACGACGc                                                                                                                                                                                                 >  1:762958/1‑280 (MQ=11)
                                                                                               aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCCTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAAAGCCGGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCgcaaccaa                                                                                                                                                                                     >  2:584240/1‑273 (MQ=255)
                                                                                                     cGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGACTTATCAGGCCTACCCCAGACAGCGCAATAGCCTGCTCTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCGCATCCGACCTCAACGCCTGATGCGACGCTTAACGCGTCTCATCCGGCCTACGCCAGACAGTCCAATAGCCTTATTTAGCGTGCTTTGCTAGGTCGGATAAGGCGTTTCGCCCGCTTCCGCTATCAACGcacga                                                                                                                                                                              >  2:746564/1‑277 (MQ=11)
                                                                                                                                ttcccccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                   <  1:1057891/276‑1 (MQ=9)
                                                                                                                                       gcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCTGCTAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTaactt                                                                                                                                            >  1:386196/1‑277 (MQ=255)
                                                                                                                                       gcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCGGATGCGACGCTTGCCGCGTCTTATCAGCCCttcctt                                                                                                                                            >  2:239930/1‑276 (MQ=9)
                                                                                                                                                                                          cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCt                                                                                                                                                 >  2:1001308/1‑225 (MQ=14)
                                                                                                                                                                                                             ctaaTTTGGCGTGATTTGTTAGGTCGGATAAGTCGTTTAGGCCGCATCCGACAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAGCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc                                                                     <  1:854315/278‑1 (MQ=11)
                                                                                                                                                                                                                                                  tccgCCGCGTCCGACATCAACGCCTGATGCGACGCTTGGCGCGTCTTATAAGGCCTACCCCAGACGACGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                <  2:839988/278‑1 (MQ=11)
                                                                                                                                                                                                                                                                      cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATATTGTAGGCCGGATAAGGCGTTTCAGCCGAACCCGCCTTCACTGCCGGATGCGACGTTTGCCGCGTCTTAGCGGGCATATCATAACCGTTGGtagat            >  2:577460/1‑277 (MQ=11)
                                                                                                                                                                                                                                                                           gATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt        >  1:1078957/1‑280 (MQ=16)
                                                                                                                                                                                                                                                                              tcGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa    <  1:287112/280‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa    <  1:894638/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                               caaCGCTTGCCGCGTATTATAAGGCCTACGCCAGACAGCGAAATGGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAGATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa     <  1:1084597/277‑1 (MQ=11)
                                                                                                                                                                                                                                                                                 aCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGACTGATGCGACGCTTAACGCGTCTTATCAGGCCTACCCAAGACAGCGCAATAGCCTGAATTAGCGTGATTTTGTAGGACGGACAAGGCGTTTACGCCGCATCCGACATCAATGACTGCCGCGACGCTTGCCCCGTCTTATAAGGCCTATCTTAACCGTTGGGTAATTTTTCCAAct  >  2:951259/1‑280 (MQ=11)
                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                 
GGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACT  >  NC_000913/4295911‑4296461

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: