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.G98.7% 295.3 / ‑2.8 75coding (65/1386 nt)yifKputative transporter YifK
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (38/36);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01

Only 100 of 105 total aligned reads displayed.

CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995
                                                  |                                            
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                               <  2:4917541/50‑1 (MQ=255)
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                               >  1:1618763/1‑50 (MQ=255)
ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt                                               <  2:541720/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                             <  1:4472415/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                             <  1:3521063/50‑1 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                             >  1:2832191/1‑50 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                             >  2:4739901/1‑50 (MQ=255)
 cAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTg                                             >  1:1707059/1‑50 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgg                                            <  1:3682758/50‑1 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgg                                            >  2:2474441/1‑50 (MQ=255)
  agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgc                                            >  2:2992501/1‑49 (MQ=255)
   gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggg                                           >  2:4223040/1‑50 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          >  1:1781701/1‑50 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          >  2:4021590/1‑50 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          >  1:2182515/1‑50 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          >  2:1870197/1‑50 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          <  1:5394577/50‑1 (MQ=255)
    agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggg                                          >  1:4423099/1‑50 (MQ=255)
     gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggg                                         <  2:2826832/50‑1 (MQ=255)
     gCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggg                                         >  2:3439615/1‑50 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                        >  2:3770317/1‑50 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                        >  2:2301661/1‑50 (MQ=255)
      cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggg                                        <  1:1016290/50‑1 (MQ=255)
       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggggg                                       >  2:2736199/1‑49 (MQ=255)
       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggggg                                       >  1:5294522/1‑49 (MQ=255)
       tACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggggg                                       <  2:4954112/50‑2 (MQ=255)
       tACAGCGTGGGCTGCAAGCTCGACATATCGAACTCATCGCCCTggggggg                                       >  2:29735/1‑49 (MQ=255)
        aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTggggggg                                       >  2:2796552/1‑48 (MQ=255)
         cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTgggggggc                                      <  2:531370/49‑3 (MQ=255)
            cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCa                                   <  1:599752/50‑1 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  2:3423723/1‑50 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  <  2:4515226/50‑1 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:2146829/1‑50 (MQ=255)
             gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAt                                  >  1:1213576/1‑50 (MQ=255)
              tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  1:3579624/50‑1 (MQ=255)
              tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt                                 <  2:2536014/50‑1 (MQ=255)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg                               >  2:2714790/1‑50 (MQ=255)
                ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCCTTTGc                              >  2:1067454/1‑50 (MQ=255)
                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                              >  2:1454550/1‑50 (MQ=255)
                 gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc                              >  1:1096480/1‑50 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             >  1:5405922/1‑50 (MQ=255)
                  cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg                             <  1:1288966/50‑1 (MQ=255)
                    ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           >  2:4704925/1‑50 (MQ=255)
                    ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           >  1:2653102/1‑50 (MQ=255)
                    ggAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc                           >  2:2285776/1‑50 (MQ=255)
                     gAAGCTCGACATATCGATCTCAGCGCCCTGGGGGGGCACCATTGGCGTc                           <  1:3211920/49‑1 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                          <  2:2566457/50‑1 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                          <  1:2670008/50‑1 (MQ=255)
                     gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCg                          <  2:4829773/50‑1 (MQ=255)
                      aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                         >  2:913538/1‑50 (MQ=255)
                      aaGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCgg                         >  2:131730/1‑50 (MQ=255)
                       aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                        <  1:2093005/50‑1 (MQ=255)
                       aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                        >  1:3622315/1‑50 (MQ=255)
                       aGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGc                        >  1:5359848/1‑50 (MQ=255)
                        gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                       <  2:5518237/50‑1 (MQ=255)
                        gCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGcc                       >  1:2823288/1‑50 (MQ=255)
                         cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt                      <  2:2246480/50‑1 (MQ=255)
                          tCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTg                     >  2:2109993/1‑50 (MQ=255)
                           cGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                    >  1:3932621/1‑50 (MQ=255)
                           cGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                    >  2:291971/1‑50 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTg                     <  2:5592354/48‑1 (MQ=39)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   <  2:2724639/50‑1 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   >  1:5538219/1‑50 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   <  2:4188491/50‑1 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   >  2:3588663/1‑50 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   <  2:5168122/50‑1 (MQ=255)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                    <  1:1221247/49‑1 (MQ=39)
                            gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGt                    <  1:4281024/49‑1 (MQ=39)
                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                  <  2:1464808/50‑1 (MQ=255)
                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGttt                  <  2:2260726/50‑1 (MQ=255)
                             aCATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt                   >  1:4982578/1‑49 (MQ=39)
                              cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                 >  2:719886/1‑50 (MQ=255)
                              cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                 >  1:2832747/1‑50 (MQ=255)
                              cATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTa                 >  2:405416/1‑50 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                <  2:2523355/50‑1 (MQ=255)
                               atatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTAt                >  1:2444554/1‑50 (MQ=255)
                                tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg               <  2:1619826/50‑1 (MQ=255)
                                tatCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATg               <  2:5691569/50‑1 (MQ=255)
                                 atCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgg              >  1:1763134/1‑50 (MQ=255)
                                  tCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggg             <  2:5616764/50‑1 (MQ=255)
                                   cGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg            >  2:5517584/1‑50 (MQ=255)
                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATggggg           >  1:2288497/1‑50 (MQ=255)
                                    gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg            <  2:2863887/49‑1 (MQ=39)
                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc          <  1:5653826/50‑1 (MQ=255)
                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc          <  2:3404350/50‑1 (MQ=255)
                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc          >  2:3305437/1‑50 (MQ=255)
                                     aaCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgc          >  2:1321997/1‑50 (MQ=255)
                                      aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  1:2968384/50‑1 (MQ=255)
                                      aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc         <  2:2319054/50‑1 (MQ=255)
                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg        >  1:2418639/1‑50 (MQ=255)
                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg        <  2:3726293/50‑1 (MQ=255)
                                       cTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccg        >  2:5421285/1‑50 (MQ=255)
                                        tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  2:4737033/1‑50 (MQ=255)
                                        tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  2:4028239/1‑50 (MQ=255)
                                        tCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  2:5288757/1‑50 (MQ=255)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       >  1:4899912/1‑48 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc       <  1:5484362/48‑1 (MQ=39)
                                          aTCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCa     <  2:3624065/50‑1 (MQ=255)
                                             gCCCTGGGGGGGCACCATTTGCGTCGGCCTGTTTATGGGTGCCGCCAGTa  >  2:1404454/1‑50 (MQ=255)
                                             gCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTa  <  2:5229945/50‑1 (MQ=255)
                                                  |                                            
CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTA  >  NC_000913/3980902‑3980995

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