Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,980,957 (G)6→7 coding (71/1386 nt) yifK → putative transporter YifK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9511.G96.0% 90.1 / ‑3.0 25coding (65/1386 nt)yifKputative transporter YifK
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (15/9);  total (16/9)
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

CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995
                                                  |                                             
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                <  1:1660843/50‑1 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgg                                             >  1:1807460/1‑50 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          >  1:1355648/1‑50 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                        >  1:2666430/1‑50 (MQ=255)
       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggggg                                       >  1:236024/1‑49 (MQ=255)
         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggc                                      <  1:1292670/49‑3 (MQ=255)
              tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  1:1135881/50‑1 (MQ=255)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                               <  1:1807474/50‑1 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             >  1:1994787/1‑50 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            <  1:798727/48‑1 (MQ=255)
                      aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                         >  1:159651/1‑50 (MQ=255)
                       aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                        >  1:154381/1‑50 (MQ=255)
                           cGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                    >  1:292145/1‑50 (MQ=255)
                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                  <  1:177565/50‑1 (MQ=255)
                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   >  1:19641/1‑49 (MQ=39)
                               atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                >  1:1891466/1‑50 (MQ=255)
                                tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg               >  1:362523/1‑50 (MQ=255)
                                tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg               >  1:2347780/1‑50 (MQ=255)
                                   cGAACTCATCGCCCTGGTGGGGCACCATTGGCGTCGGCCTGTTTATgggg            >  1:683354/1‑50 (MQ=39)
                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg           >  1:1222781/1‑50 (MQ=255)
                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg            <  1:859963/49‑1 (MQ=39)
                                      aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:2041844/50‑1 (MQ=255)
                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg        >  1:2075815/1‑50 (MQ=255)
                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg        <  1:2729306/50‑1 (MQ=255)
                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCAGCCTGTTTATGGGGgccg        >  1:1249424/1‑50 (MQ=255)
                                        tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:2052664/1‑50 (MQ=255)
                                         cATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:2572790/1‑49 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:2778904/1‑48 (MQ=39)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt   <  1:2337591/49‑1 (MQ=39)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt   <  1:175425/49‑1 (MQ=39)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  >  1:1239036/1‑50 (MQ=255)
                                                  |                                             
CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995

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