Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 19.6% intergenic (+266/+376) 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,0600CT19.6% 70.7 / 16.2 46intergenic (+266/+376)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (24/13);  new base T (4/5);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

TTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTA  >  NC_000913/4295926‑4296200
                                                                                                                                      |                                                                                                                                            
tgctGCCGCGCCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATg                                                                                                                                      <  2:272647/140‑1 (MQ=12)
     ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                                                                                                 >  1:502855/1‑143 (MQ=17)
        cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                                                                                              >  1:501719/1‑143 (MQ=17)
          cccGGCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                            <  1:463281/142‑1 (MQ=21)
           ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCACCGCCTGACGCGACGCTTGcc                                                                                                                           >  2:396862/1‑143 (MQ=21)
               catgaccgcgTGCTGCGACGCTTAGCCGGTCTTATCAGGCGTACGCGAGACGGCGCAATAGCCTGATTTCGGGTGATTTTGTAGGTCGTATAAGGCGTGTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                       <  2:45771/137‑1 (MQ=11)
                aTCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                      <  2:256880/143‑1 (MQ=255)
                  cAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                    <  2:530805/143‑1 (MQ=255)
                   aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                   >  1:355894/1‑143 (MQ=25)
                       ccTGATGCGACGCTTAACGCGTCCTATCAGGCCCACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGTCGCGTGttttttg                                                                                                               <  1:150596/143‑6 (MQ=11)
                               gACGCTGTACTCGTCTTATCAGGTCTACGTGCGACCTCGCATTAGCATGATTTCGTGTCATTTTGTAGGTCGGACATGTCGTTTCCACCCGATTCGATATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                       <  2:51854/143‑1 (MQ=2)
                                aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctt                                                                                                      >  2:35096/1‑139 (MQ=9)
                                  gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCCGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAg                                                                                                    >  1:537210/1‑143 (MQ=38)
                                   ctgcaCGCGTCTTAGCGGGCCTACGCGAGACACCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAGGGGGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                   <  2:161773/139‑1 (MQ=11)
                                       aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                >  2:65312/1‑126 (MQ=40)
                                       aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                <  1:65312/126‑1 (MQ=40)
                                            tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCGGGCGgcgccaa                                                                                          >  2:413342/1‑140 (MQ=255)
                                             cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCCATa                                                                                         >  1:579875/1‑143 (MQ=38)
                                                      ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAGATCAACGCCTGATGCGAAGATTGGCGCGTCTTATTAGGGCTATGGACGACAGCCGAATAGCGTGAttt                                                                                >  1:294685/1‑143 (MQ=11)
                                                           gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCGTACGGCAGACAGCCCTGTAGCGTGATTTAgggg                                                                           >  2:150596/1‑140 (MQ=11)
                                                           gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                      <  1:23730/100‑1 (MQ=255)
                                                           gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                      >  2:23730/1‑100 (MQ=255)
                                                           gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                <  2:176074/106‑1 (MQ=255)
                                                           gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttcaccgg                                                                           >  2:165529/1‑137 (MQ=255)
                                                                 cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtttt                                                                     <  1:60213/143‑1 (MQ=38)
                                                                       aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCCCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGCCAGCGCATTAGCCTGCTATATCGTGATTGTGTAgtt                                                               >  2:181388/1‑141 (MQ=14)
                                                                              tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGCTac                                                        >  1:594649/1‑142 (MQ=255)
                                                                                                    tCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                       <  1:604117/58‑1 (MQ=255)
                                                                                                    tCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                       <  1:598560/58‑1 (MQ=255)
                                                                                                    tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                       >  2:598560/1‑58 (MQ=255)
                                                                                                    tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                       >  2:604117/1‑58 (MQ=255)
                                                                                                    tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATCGCCTGATTTAGCGGGAGTTTGTAGGTCGGAGAAGGCGTGTAGGCCGCAGCCGACa                                  >  2:24535/1‑143 (MQ=14)
                                                                                                    tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGCTTCAGCGTGAGTGTGCAGGTAGGATGACGCGTACAAGCCGCAGCGGACa                                  >  2:223267/1‑143 (MQ=11)
                                                                                                      ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATc                                >  2:149242/1‑143 (MQ=39)
                                                                                                      ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATGTAGCGTCATTTTGTGGGGCGGATCCGGCGTTTATGTCGCTTGCGACATc                                >  1:75562/1‑143 (MQ=14)
                                                                                                            ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGCTAAGGCGTTTGTGCCGCATCCGACATCAACGcc                          >  2:322019/1‑143 (MQ=255)
                                                                                                               gTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                       >  1:10072/1‑143 (MQ=37)
                                                                                                                  tACGCCGCATCCGACATCAATGCCTGATGCGACGc                                                                                                                                <  2:65807/35‑1 (MQ=1)
                                                                                                                  tACGCCGCATCCGACATCAACGCCTGATGCGACGc                                                                                                                                >  1:65807/1‑35 (MQ=40)
                                                                                                                  tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGc                    >  1:252084/1‑143 (MQ=38)
                                                                                                                              ggcgtCACCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc        <  2:411667/139‑1 (MQ=14)
                                                                                                                              gACATCAACGCCTGATGCGACGCTTGCCGCGTCTGATCAGGCCTACGCCAGACAGCGCATTGGGGTGGTTGAACAACATTTTGTCGGTCAGGTAATGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc        <  1:149242/143‑1 (MQ=11)
                                                                                                                                  tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                             <  1:526472/38‑1 (MQ=40)
                                                                                                                                  tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACa                                  >  2:276940/1‑113 (MQ=18)
                                                                                                                                   cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTGGCGCGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACGCCGGATGCGACGCTGATCGCGTCtt   >  1:11878/1‑143 (MQ=11)
                                                                                                                                    aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa  >  1:555222/1‑143 (MQ=17)
                                                                                                                                      |                                                                                                                                            
TTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTA  >  NC_000913/4295926‑4296200

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