Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,381:1 +GC 100% intergenic (+587/+55) 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,3801.C96.2% 164.0 / ‑0.2 53intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (22/28);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,296,3802.G96.2% 168.0 / 0.3 53intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (23/28);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296343‑4296470
                                         ||                                                                                          
tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                <  2:208187/100‑1 (MQ=255)
tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                <  2:1338984/100‑1 (MQ=255)
tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                <  2:428340/100‑1 (MQ=255)
 gATTTAGCGTTATTTTGTAGGACGACTCAGGCGTGTTC‑T‑‑CCCCGTAAGGTATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                <  1:1101803/99‑1 (MQ=2)
 gATTTAGCGTGATTTTGTAGGTCGGATAAGGGGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                               <  1:2249712/100‑1 (MQ=11)
    ttAGCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                            <  2:2682393/100‑1 (MQ=255)
     tAGCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                           <  2:105810/100‑1 (MQ=255)
      aGCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                          <  1:1014212/100‑1 (MQ=255)
      aGCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                          <  2:39267/100‑1 (MQ=255)
       gCGTGATTTTGTAGGTCGGCTAAGACGT‑TTCCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  1:282880/100‑1 (MQ=11)
       gCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATGTTagacggcgca                         >  1:356456/1‑91 (MQ=11)
       gCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  1:491134/100‑1 (MQ=255)
        cttGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                              <  1:1037475/92‑1 (MQ=11)
         gTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  1:1639695/1‑100 (MQ=255)
         gTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       <  1:967168/100‑1 (MQ=255)
         gTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  2:2950395/1‑100 (MQ=255)
          tGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                      >  1:1393071/1‑100 (MQ=255)
          tGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  2:3098778/1‑99 (MQ=255)
            aTTTTGTAGGTCGGATAGGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                     >  1:3110712/1‑99 (MQ=14)
            aTTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                    >  1:2538717/1‑100 (MQ=255)
            aTTTTGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                    <  2:737953/100‑1 (MQ=255)
             ttttGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                   <  2:3072629/100‑1 (MQ=255)
               ttGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                 <  1:1755745/100‑1 (MQ=255)
                tGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  2:2593271/100‑1 (MQ=255)
                tGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  1:2594269/100‑1 (MQ=255)
                tGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  2:1377546/100‑1 (MQ=255)
                tGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                >  1:9756/1‑100 (MQ=255)
                tGTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  2:2807632/100‑1 (MQ=255)
                 gTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTATCTTAACCGTTGGTTAATTTTTCa               <  1:2249886/100‑1 (MQ=17)
                 gTAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa               <  1:1426856/100‑1 (MQ=255)
                  tAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa              >  2:134857/1‑100 (MQ=255)
                  tAGGTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa              <  1:2171371/100‑1 (MQ=255)
                    ggTCGGATAAGGGGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            <  2:1491858/100‑1 (MQ=17)
                    ggTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  1:2111477/1‑100 (MQ=255)
                    ggTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  1:2695861/1‑100 (MQ=255)
                    ggTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  1:2345977/1‑100 (MQ=255)
                    ggTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  2:2567796/1‑100 (MQ=255)
                    ggTCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  2:2438062/1‑100 (MQ=255)
                      tCGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc          <  1:3145297/100‑1 (MQ=255)
                       cGGATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct         <  1:298676/100‑1 (MQ=255)
                        ggATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        >  2:990919/1‑100 (MQ=255)
                        ggATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        >  2:2122451/1‑100 (MQ=255)
                        ggATAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        <  1:3300020/100‑1 (MQ=255)
                          aTAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGTTTAATTTTTCAAACTCTTca      <  1:2122434/100‑1 (MQ=18)
                           tAAGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac     >  1:2530660/1‑100 (MQ=255)
                            aaGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca    >  2:1839489/1‑100 (MQ=255)
                            aaGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca    >  2:2600511/1‑100 (MQ=255)
                             aGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac   >  1:3238624/1‑100 (MQ=255)
                             aGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac   >  1:1849339/1‑100 (MQ=255)
                             aGGCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCAACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac   >  1:1848461/1‑100 (MQ=21)
                              ggCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTTATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc  <  1:1553075/100‑1 (MQ=21)
                              ggCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc  >  2:1597494/1‑100 (MQ=255)
                              ggCGT‑TT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc  <  1:3249120/100‑1 (MQ=255)
                                         ||                                                                                          
TGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGT‑TT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296343‑4296470

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