Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,978,702 C→T S68F (TCC→TTC)  yifK → putative amino acid/amine transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,978,7020CT99.3% 452.1 / NA 144S68F (TCC→TTC) yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (81/62);  total (81/63)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 147 total aligned reads displayed.

TCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTCCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCAC  >  NC_000913/3978629‑3978774
                                                                         |                                                                        
  cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCtg                                                                      <  2:5041507/76‑2 (MQ=255)
  cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTc                                                                      >  1:5493029/1‑76 (MQ=255)
     aTTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGaa                                                                   <  1:573232/76‑1 (MQ=255)
     aTTGTGGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGaa                                                                   <  2:1175739/76‑1 (MQ=255)
      ttgttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAac                                                                  >  2:4951757/1‑76 (MQ=255)
       tgttgtCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAacc                                                                 <  1:1918938/76‑1 (MQ=255)
       tgttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAacc                                                                 >  1:262409/1‑76 (MQ=255)
        gttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGCaccg                                                                >  2:310550/1‑76 (MQ=255)
        gttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccg                                                                >  1:4483754/1‑76 (MQ=255)
           ggCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggtt                                                             >  2:4494682/1‑76 (MQ=255)
           ggCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggtt                                                             <  1:3510138/76‑1 (MQ=255)
            gtcTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATCGGCGAAATGTTGTTCCTCGAaccggtta                                                            <  2:2547953/74‑1 (MQ=255)
            gCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttc                                                            >  2:4775524/1‑75 (MQ=255)
              cTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttacc                                                          <  2:4143168/76‑1 (MQ=255)
               tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccg                                                         <  1:2064976/76‑1 (MQ=255)
               tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccg                                                         <  2:2433838/76‑1 (MQ=255)
                atatCATCGCGGGGCTGTTCGTCTTTTTCATCATGCGTTCAAGGGGCGAAATGTTGTTCCTCGAaccggttaccgg                                                        <  2:5390483/76‑1 (MQ=255)
                 tatCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggt                                                       >  1:972166/1‑76 (MQ=255)
                  atcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                      >  1:962224/1‑76 (MQ=255)
                   tcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttc                                                     >  1:2949977/1‑76 (MQ=255)
                   gtgtcGCCGGGCTGTTCGTCTTTTTCATCGTGCGTTCAATGGGCGACATGTTGTTCCTGGAACCGGGTACCGgttc                                                     <  2:1837983/73‑1 (MQ=255)
                    catcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcg                                                    >  2:5467442/1‑76 (MQ=255)
                     atcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgt                                                   >  1:3676253/1‑76 (MQ=255)
                     atcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgt                                                   <  1:5252828/76‑1 (MQ=255)
                     atcGACGGGCTGTTCGTCTTTTTCATCATGCGTTCAATCACAGAAATGTTGTTACTCGAACAGGTTACCGgttagt                                                   >  2:2458897/1‑76 (MQ=255)
                      tcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgtt                                                  <  1:5701000/76‑1 (MQ=255)
                       cGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgttc                                                 >  1:3359381/1‑76 (MQ=255)
                        ggcGGGCTTTCCGTCTTTTTCGTCATGCGTTCATTGGGCGAGATGTTGTTCTTCGAACCGGTTTCCGGttcgttcg                                                <  2:2771123/74‑1 (MQ=255)
                        gCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGttcgttcg                                                <  1:4631931/76‑1 (MQ=255)
                        gCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGttcgttcg                                                >  2:4371083/1‑76 (MQ=255)
                        gCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGttcgttcg                                                >  2:389455/1‑76 (MQ=255)
                        gCCGGGCTGTTCGTCTTTTTCATCATGAGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGgttgttcg                                                >  1:78184/1‑76 (MQ=255)
                          cGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCTGTTACCGGTTCGTTCGcc                                              >  1:3539871/1‑76 (MQ=255)
                           gggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCg                                             >  1:4946728/1‑76 (MQ=255)
                           gggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCg                                             >  1:1793204/1‑76 (MQ=255)
                           gggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCg                                             <  1:5153512/76‑1 (MQ=255)
                           gggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCg                                             >  1:4005041/1‑76 (MQ=255)
                           gggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCg                                             >  1:2562665/1‑76 (MQ=255)
                             gCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGtt                                           <  2:5493027/76‑1 (MQ=255)
                               tGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTa                                         >  2:1678110/1‑76 (MQ=255)
                                gTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                        <  1:3612735/76‑1 (MQ=255)
                                gTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                        >  2:4208207/1‑76 (MQ=255)
                                 ttCGTCTTTTTCATCATGCGTACAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATg                                       >  2:3760807/1‑76 (MQ=255)
                                  tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgc                                      <  2:442340/76‑1 (MQ=255)
                                  tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgc                                      >  1:4733138/1‑76 (MQ=255)
                                  tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgc                                      <  1:3952739/76‑1 (MQ=255)
                                    gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgcgc                                    >  1:4179430/1‑76 (MQ=255)
                                    gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgcgc                                    >  1:5572441/1‑76 (MQ=255)
                                     tgtttttGAGCATGCGTTCAATGGGGGAGATGTTGTTCCTCGAACCGGTTACCGGTTTGTTCGCCGTTTATGCGCa                                   <  2:78184/74‑1 (MQ=255)
                                     tCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCa                                   >  1:4423730/1‑76 (MQ=255)
                                      cTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccgg                                                        >  2:5666760/1‑54 (MQ=255)
                                       tttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATc                                 >  2:3297239/1‑76 (MQ=255)
                                       tttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATc                                 <  2:5664519/76‑1 (MQ=255)
                                        ttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCg                                <  2:4958071/76‑1 (MQ=255)
                                        ttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCg                                >  2:5005109/1‑76 (MQ=255)
                                          ttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGtt                              >  1:4901291/1‑76 (MQ=255)
                                          ttCATCATGCGGTCAATGGGCGAAATGTTGGTCCTCGAACCGGTTACCGGATCGTTCGCACATTATGCGCAGCGtt                              >  2:2716245/1‑76 (MQ=255)
                                           tcatcaTGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGttt                             >  2:1019961/1‑75 (MQ=255)
                                            catcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtat                            >  1:2192464/1‑76 (MQ=255)
                                             atcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtata                           >  2:2195287/1‑76 (MQ=255)
                                              tcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTatat                          >  1:1160990/1‑76 (MQ=255)
                                                 tGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAg                       >  2:1871001/1‑76 (MQ=255)
                                                  gCGTTCAATGGTCGACATGTTGTTCCTCGAACCGGTTCCCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGc                      <  2:3639466/76‑1 (MQ=255)
                                                  gCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGc                      >  2:1909688/1‑76 (MQ=255)
                                                  gCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGc                      >  1:1867533/1‑76 (MQ=255)
                                                    gTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCTGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGccc                    >  1:4004660/1‑76 (MQ=255)
                                                    gTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGccc                    >  1:4249479/1‑76 (MQ=255)
                                                    gTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGccc                    >  2:103745/1‑76 (MQ=255)
                                                    gTTCAATGGGAGAAATGTTGTTCCTCGAACCGGTTACCGGGTCGTTCGCCGTTTATGAGCATCGCTATATGAGccc                    >  1:4780262/1‑76 (MQ=255)
                                                     ttCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCg                   >  2:5090500/1‑76 (MQ=255)
                                                     ttCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCg                   <  2:843473/76‑1 (MQ=255)
                                                     ttCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCg                   >  2:2535528/1‑76 (MQ=255)
                                                     ttCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCg                   <  2:2244778/76‑1 (MQ=255)
                                                     ttCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCg                   >  2:1484568/1‑76 (MQ=255)
                                                      tCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGt                  >  2:3412836/1‑76 (MQ=255)
                                                      tCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGt                  >  2:1441522/1‑76 (MQ=255)
                                                       cAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGtt                 <  2:463433/76‑1 (MQ=255)
                                                        aaTGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTc                >  2:2978813/1‑76 (MQ=255)
                                                        aaTGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTc                <  1:2923174/76‑1 (MQ=255)
                                                         aTGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCt               <  1:4668965/76‑1 (MQ=255)
                                                         aTGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCt               >  1:1024366/1‑76 (MQ=255)
                                                          tGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCtt              >  2:4628833/1‑76 (MQ=255)
                                                          tGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCtt              >  2:5290071/1‑76 (MQ=255)
                                                           gggCGAAATGTTGTTCCTCGAACCGGTTACCGgttcg                                                    <  1:2179702/37‑1 (MQ=39)
                                                           gggCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCttt             <  2:1949053/76‑1 (MQ=255)
                                                               gAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCt         >  1:2079233/1‑76 (MQ=255)
                                                               gAAATGTTGTCCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCt         <  1:1032948/76‑1 (MQ=255)
                                                                aaaTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTa        >  1:722197/1‑76 (MQ=255)
                                                                aaaTGCTGTTCCTCGAACAGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTGCTATCGCTa        >  2:3422372/1‑76 (MQ=255)
                                                                 aaTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAt       <  1:322904/76‑1 (MQ=255)
                                                                  aTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCtgtt      >  1:2642136/1‑73 (MQ=255)
                                                                  aTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAtc      <  2:4370335/76‑1 (MQ=255)
                                                                  aTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTCGCTAtc      >  2:5075179/1‑76 (MQ=255)
                                                                    gttgttCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAtctc    <  1:1900727/76‑1 (MQ=255)
                                                                    gttgttCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAtctc    >  2:1302961/1‑76 (MQ=255)
                                                                     ttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa   >  1:3006558/1‑76 (MQ=255)
                                                                     ttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa   >  1:2589580/1‑76 (MQ=255)
                                                                     ttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa   <  1:3519933/76‑1 (MQ=255)
                                                                     ttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa   <  1:4565558/76‑1 (MQ=255)
                                                                      tgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCAc  <  2:3955650/76‑1 (MQ=255)
                                                                         |                                                                        
TCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTCCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCAC  >  NC_000913/3978629‑3978774

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: