Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,381 +GC 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.C93.0% 134.2 / ‑2.2 44intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (20/20);  minor base T (1/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009134,296,3802.G97.7% 146.3 / ‑3.0 44intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (21/22);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296341‑4296470
                                          ||                                                                                          
cggggttcaGCGTGATTTTGTCGGTCGGACACGCCGTTCCTG‑‑CCGCATCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                       <  1:339884/87‑1 (MQ=11)
   cttccaatctGATTTTGTAGGTCGGATAAGGCGTT‑‑TATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                               <  2:2616908/91‑1 (MQ=11)
       tAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                           <  2:1188746/100‑1 (MQ=255)
         gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACTTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  2:60661/100‑1 (MQ=14)
         gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         >  2:1704571/1‑100 (MQ=255)
         gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  2:140895/100‑1 (MQ=255)
         gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  2:60459/100‑1 (MQ=255)
         gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  2:60435/100‑1 (MQ=255)
         gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  1:2413380/100‑1 (MQ=255)
          cGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        >  2:1510531/1‑100 (MQ=255)
          cGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        <  2:2431101/100‑1 (MQ=255)
           gTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATGTTAACCGTTGGTTa                       >  1:1921004/1‑100 (MQ=11)
           gTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  1:269844/1‑100 (MQ=255)
           gTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  1:1550050/1‑100 (MQ=255)
            tGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                      >  2:516123/1‑100 (MQ=255)
             gATTTTGTAGGTCGGATAAGGCGTT‑‑TATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAccagatcggaag                     >  1:2616887/1‑90 (MQ=11)
             gATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                     <  2:220860/100‑1 (MQ=255)
             gATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                     >  2:936934/1‑100 (MQ=255)
               ttttGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                   >  2:1821823/1‑100 (MQ=255)
                tttGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  <  1:984774/100‑1 (MQ=255)
                tttGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  >  1:2460219/1‑100 (MQ=255)
                 ttGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                 <  2:2745012/100‑1 (MQ=255)
                  tGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  2:978827/100‑1 (MQ=255)
                  tGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  2:1315959/100‑1 (MQ=255)
                  tGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTCCCGCGTCTTAGCAGGCCTATCTTAACCGTTGGTTAAttttcc                <  2:1906435/100‑3 (MQ=17)
                    tAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa              <  2:2561033/100‑1 (MQ=255)
                     aGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa             >  1:184757/1‑100 (MQ=255)
                      ggTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  2:958680/1‑100 (MQ=255)
                      ggTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  2:2515542/1‑100 (MQ=255)
                       gTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct           >  2:1423936/1‑100 (MQ=255)
                       gTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct           >  1:459987/1‑100 (MQ=255)
                       gTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct           >  1:2332404/1‑100 (MQ=255)
                        tCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc          <  2:513938/100‑1 (MQ=255)
                        tCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTCATCCGGCCTATCTTCACCGTTGGTTAATTTTTCAAActc          >  2:1373858/1‑100 (MQ=17)
                         cGGATAAGGCGTT‑‑TATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct         <  1:1706480/100‑1 (MQ=255)
                         cGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct         >  2:1552326/1‑100 (MQ=255)
                           gATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  1:449815/100‑1 (MQ=255)
                           gATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  1:1340974/100‑1 (MQ=255)
                           gATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  2:1550062/100‑1 (MQ=255)
                            aTAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca      <  2:695413/100‑1 (MQ=255)
                              aaGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca    <  1:555354/100‑1 (MQ=255)
                              aaGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca    >  2:1357655/1‑100 (MQ=255)
                                ggCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc  >  2:2760774/1‑100 (MQ=255)
                                ggCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAACCTTTCAAACTCTTCACAcc  >  1:1069338/1‑100 (MQ=18)
                                          ||                                                                                          
CCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296341‑4296470

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