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.4% 172.9 / 0.8 56intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (26/27);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
*NC_0009134,296,3802.G96.4% 175.9 / 1.2 56intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (26/28);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296344‑4296470
                                       ||                                                                                            
gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                               >  2:3420746/1‑100 (MQ=255)
 aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                              >  1:3266955/1‑100 (MQ=255)
  tttCGCGTGATTTTGTATGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                               <  1:2507672/98‑1 (MQ=11)
  tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                               <  1:211394/98‑1 (MQ=255)
   ttAGCGTGATTTTGTTGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                            <  2:3502090/100‑1 (MQ=11)
   ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                            >  2:2024793/1‑100 (MQ=255)
    tAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                            <  1:3108898/99‑1 (MQ=255)
    tAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                           <  1:1712633/100‑1 (MQ=255)
    tAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                           <  2:437512/100‑1 (MQ=255)
     aGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                          <  1:1766382/100‑1 (MQ=255)
      gCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGTTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         >  2:2496053/1‑100 (MQ=14)
      gCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                          >  1:1080929/1‑99 (MQ=255)
      gCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         >  1:2923308/1‑100 (MQ=255)
       cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        >  2:2228154/1‑100 (MQ=255)
       cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        >  2:369052/1‑100 (MQ=255)
        gTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        >  1:16358/1‑99 (MQ=255)
        gTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  1:2600974/1‑100 (MQ=255)
        gTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CGCGGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  2:1567543/1‑100 (MQ=14)
         tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGCCATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                      >  2:1034019/1‑100 (MQ=14)
         tGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                      >  1:3168516/1‑100 (MQ=255)
          gATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                     >  1:950375/1‑100 (MQ=255)
          gATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                     >  1:950397/1‑100 (MQ=255)
           aTTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                    <  2:3275142/100‑1 (MQ=255)
             tttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  >  2:1914286/1‑100 (MQ=255)
              ttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                 >  2:745218/1‑100 (MQ=255)
              ttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                 <  2:128426/100‑1 (MQ=255)
               tGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                >  1:1978269/1‑100 (MQ=255)
               tGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  1:1215222/100‑1 (MQ=255)
                gTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa               <  1:2088137/100‑1 (MQ=255)
                gTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa               <  2:2381438/100‑1 (MQ=255)
                gTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa               <  2:911708/100‑1 (MQ=255)
                 tAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCan              <  2:1211410/100‑2 (MQ=17)
                 tAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa              <  2:465417/100‑1 (MQ=255)
                 tAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa              <  1:998041/100‑1 (MQ=255)
                   ggTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa              <  1:828607/98‑1 (MQ=255)
                   ggTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  2:990284/1‑100 (MQ=255)
                   ggTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  1:1567713/1‑100 (MQ=255)
                   ggTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  1:42134/1‑100 (MQ=255)
                   ggTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            <  1:3381544/100‑1 (MQ=255)
                    gTCGGATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct           >  1:440288/1‑100 (MQ=255)
                     tCGGCTAAGGCGTTTCCG‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct           <  1:2810432/99‑1 (MQ=17)
                       ggATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        <  1:3539507/100‑1 (MQ=255)
                       ggATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        >  1:211894/1‑100 (MQ=255)
                       ggATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        >  2:1091631/1‑100 (MQ=255)
                       ggATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCACGCCTATCTTAACCGTTCGTTAACTTTTCAAActcct        >  1:1713593/1‑98 (MQ=17)
                        gATAAGGGGTTT‑‑ATGC‑C‑GCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  1:3181921/100‑1 (MQ=17)
                        gATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        <  1:3250941/99‑1 (MQ=255)
                        gATAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  2:3052909/100‑1 (MQ=255)
                         aTAAGTCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCTCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        <  1:3375570/98‑1 (MQ=18)
                         aTAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca      <  1:2211384/100‑1 (MQ=255)
                         aTAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  1:336999/99‑1 (MQ=255)
                          tAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac     <  2:2334222/100‑1 (MQ=255)
                          tAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac     <  1:3266485/100‑1 (MQ=255)
                          tAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac     <  1:3050427/100‑1 (MQ=255)
                          tAAGGCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       >  1:3499199/1‑98 (MQ=18)
                             ggCGTTT‑‑ACGC‑C‑GCATCCGACATCAATGCCTGATGCGGCGCTTGCCGCGTCTTATCAGACCTATCTTAAACGTTGCTCAATTTTTCAAACTCTTCACAcc  >  1:3062345/1‑100 (MQ=17)
                                       ||                                                                                            
GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296344‑4296470

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