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.G97.9% 186.0 / ‑2.7 47coding (65/1386 nt)yifKputative transporter YifK
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (18/28);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995
                                                  |                                             
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                                <  1:3660391/50‑1 (MQ=255)
ccAGAGCTACAGCGTGGGCTGGAGCTCGACATATCGAACTCATCGCCCTg                                              >  1:356379/1‑50 (MQ=38)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                              <  1:1791375/50‑1 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgg                                             <  1:4328323/50‑1 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgg                                             <  1:4780741/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          <  1:5062004/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          <  1:2109598/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          <  1:4438405/50‑1 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                        <  1:2911877/50‑1 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                        <  1:1170539/50‑1 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                        >  1:4792671/1‑50 (MQ=255)
       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggggg                                       >  1:470551/1‑49 (MQ=255)
        aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggggg                                       >  1:1636366/1‑48 (MQ=255)
         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggc                                      <  1:2898910/49‑3 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  <  1:4405100/50‑1 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:2938950/1‑50 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:1147411/1‑50 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:3816467/1‑50 (MQ=255)
              tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  1:476944/50‑1 (MQ=255)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTg                                <  1:5897965/49‑1 (MQ=39)
                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                              <  1:2783144/50‑1 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             <  1:5434231/50‑1 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             <  1:4566127/50‑1 (MQ=255)
                   tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            >  1:4926963/1‑50 (MQ=255)
                   tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            >  1:1905742/1‑50 (MQ=255)
                   tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            >  1:4529898/1‑50 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt                            <  1:933205/48‑1 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           <  1:1047769/49‑1 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           <  1:1099047/49‑1 (MQ=255)
                      aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                         >  1:5222819/1‑50 (MQ=255)
                      aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                         <  1:917986/50‑1 (MQ=255)
                      aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                         <  1:1083833/50‑1 (MQ=255)
                      aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                         >  1:4116370/1‑50 (MQ=255)
                        gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                       >  1:3818477/1‑50 (MQ=255)
                        gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                       >  1:1579312/1‑50 (MQ=255)
                         cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                      >  1:52066/1‑50 (MQ=255)
                           cGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                    >  1:3627928/1‑50 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   <  1:4171502/50‑1 (MQ=255)
                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                  <  1:3724877/50‑1 (MQ=255)
                              cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                 <  1:5176434/50‑1 (MQ=255)
                              cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                 >  1:1242815/1‑50 (MQ=255)
                              cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                 <  1:1155908/50‑1 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                <  1:1095453/50‑1 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                >  1:5395832/1‑50 (MQ=255)
                                tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg               <  1:5819192/50‑1 (MQ=255)
                                 atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg              >  1:380194/1‑50 (MQ=255)
                                  tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg             >  1:2695830/1‑50 (MQ=255)
                                      aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:3248566/50‑1 (MQ=255)
                                      aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         >  1:4395527/1‑50 (MQ=255)
                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg        <  1:4107410/50‑1 (MQ=255)
                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg        <  1:4347142/50‑1 (MQ=255)
                                         cATCGCCCTGGTGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:3527411/1‑49 (MQ=38)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgcc      <  1:935019/49‑1 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:570152/50‑1 (MQ=255)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:2594350/50‑1 (MQ=255)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:2805380/50‑1 (MQ=255)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  1:4120967/50‑1 (MQ=255)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt   <  1:2055919/49‑1 (MQ=39)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt   <  1:5780600/49‑1 (MQ=39)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt   <  1:3680464/49‑1 (MQ=39)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  >  1:3209573/1‑50 (MQ=255)
                                                  |                                             
CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995

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