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% 426.3 / NA 139S68F (TCC→TTC) yifKputative amino acid/amine transport protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (76/62);  minor base G (0/1);  total (76/63)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 140 total aligned reads displayed.

CCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTCCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCA  >  NC_000913/3978630‑3978773
                                                                        |                                                                       
ccGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCt                                                                      >  2:1148561/1‑76 (MQ=255)
 cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTc                                                                     <  1:3805721/76‑1 (MQ=255)
 cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTc                                                                     >  1:535834/1‑76 (MQ=255)
 cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTc                                                                     <  2:1663799/76‑1 (MQ=255)
 cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTc                                                                     >  1:287205/1‑76 (MQ=255)
 cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTc                                                                     >  2:3661988/1‑76 (MQ=255)
 cGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTc                                                                     <  2:4550859/76‑1 (MQ=255)
  gtgttgttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCg                                                                    <  2:1972590/73‑1 (MQ=255)
  gTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCg                                                                    >  2:2797825/1‑76 (MQ=255)
  gTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCg                                                                    <  2:2049146/76‑1 (MQ=255)
    aTTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGaa                                                                  <  2:4294674/76‑1 (MQ=255)
    aTTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGaa                                                                  <  1:4315121/76‑1 (MQ=255)
    aTTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGaa                                                                  <  1:126951/76‑1 (MQ=255)
    aTTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGaa                                                                  <  1:468398/76‑1 (MQ=255)
     ttgttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGCAATGTTGTTCCTCGAac                                                                 >  2:4082635/1‑76 (MQ=255)
      tgttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAacc                                                                >  1:1870606/1‑76 (MQ=255)
      tgttgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAacc                                                                <  1:4085678/76‑1 (MQ=255)
         tgGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggt                                                             <  2:3019869/76‑1 (MQ=255)
              tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccg                                                        <  1:46219/76‑1 (MQ=255)
              tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccg                                                        <  1:3379990/76‑1 (MQ=255)
               atatCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccgg                                                       >  2:4193222/1‑76 (MQ=255)
               atatCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGCTGTTCCTCGAaccggttaccgg                                                       >  1:4189548/1‑76 (MQ=255)
                tatCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggt                                                      <  1:2684167/76‑1 (MQ=255)
                 atcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                     >  1:931082/1‑76 (MQ=255)
                 atcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                     <  2:4757346/76‑1 (MQ=255)
                 atcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                     >  1:1378014/1‑76 (MQ=255)
                  tcatcGCCGGGCTGTTCGTCTTTTTCCTCATGCGTTCAATGGGCGAAATGTTGTTCCCCGAACCGGTTACCGgttc                                                    >  2:4837431/1‑76 (MQ=255)
                  tcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttc                                                    <  1:4093709/76‑1 (MQ=255)
                  tcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttc                                                    <  1:4344137/76‑1 (MQ=255)
                     tcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgtt                                                 <  2:3455996/76‑1 (MQ=255)
                      cGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgttc                                                >  1:2991877/1‑76 (MQ=255)
                      cGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgttc                                                <  1:3508669/76‑1 (MQ=255)
                      cGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgttc                                                <  1:2964637/76‑1 (MQ=255)
                      cGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgttc                                                <  1:2322542/76‑1 (MQ=255)
                       gCCTGGCTGTTCGTCTTTTTCATCATTCTTTCAATGCGCGAAATGTCTTTCCTCGAACCGGTTACCGgctccttcc                                               >  2:1457423/1‑75 (MQ=255)
                       gCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGttcgttcg                                               <  2:897634/76‑1 (MQ=255)
                        ccGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGc                                              >  2:655691/1‑76 (MQ=255)
                        ccGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGc                                              >  1:2889941/1‑76 (MQ=255)
                         cGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGcc                                             <  1:394435/76‑1 (MQ=255)
                         cGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGcc                                             >  2:1369034/1‑76 (MQ=255)
                         cGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGcc                                             >  1:480883/1‑76 (MQ=255)
                         cGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGcc                                             >  2:573818/1‑76 (MQ=255)
                          gggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCg                                            >  2:858835/1‑76 (MQ=255)
                            gCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGACCCGGTTACCGGTTCGTTCGCCGtt                                          <  1:3856954/76‑1 (MQ=255)
                            gCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGtt                                          >  1:1847398/1‑76 (MQ=255)
                              tGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTa                                        >  1:17566/1‑76 (MQ=255)
                              tGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTa                                        >  2:3365959/1‑76 (MQ=255)
                               gTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                       >  2:562309/1‑76 (MQ=255)
                               gTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAt                                       <  1:1004950/76‑1 (MQ=255)
                                ttCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGGACCGGTTACCGGTTCGTTCGCCGTTTATg                                      >  1:2486845/1‑76 (MQ=255)
                                ttCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATg                                      <  2:1184850/76‑1 (MQ=255)
                                 tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAtcc                                     >  2:343399/1‑74 (MQ=255)
                                 tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATATTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgc                                     >  1:667117/1‑76 (MQ=255)
                                  cGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTCATgcc                                    >  2:568169/1‑75 (MQ=255)
                                   gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGATCCGGTTACCGGTTCGTTCGCCGTTTATgcgc                                   <  1:959492/76‑1 (MQ=255)
                                   gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgcgc                                   <  1:52881/76‑1 (MQ=255)
                                   gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgcgc                                   >  1:1315859/1‑76 (MQ=255)
                                   gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgcgc                                   >  2:1548209/1‑76 (MQ=255)
                                      tttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCCTc                                >  2:63934/1‑76 (MQ=255)
                                      tttttCATCATGCGTTCAAGGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTCCGCCGTTTATGCGCATc                                >  1:4554883/1‑76 (MQ=255)
                                       ttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCg                               >  1:2105202/1‑76 (MQ=255)
                                          tcatcaTGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTta                            >  2:937323/1‑76 (MQ=255)
                                           catcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtat                           >  2:2409382/1‑76 (MQ=255)
                                           catcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtat                           >  2:2387207/1‑76 (MQ=255)
                                            atcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtata                          >  2:4827884/1‑76 (MQ=255)
                                            atcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtata                          <  1:1194313/76‑1 (MQ=255)
                                            atcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtata                          >  2:1856430/1‑76 (MQ=255)
                                             tcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTatat                         >  2:1862060/1‑76 (MQ=255)
                                             tcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTatat                         >  2:3082615/1‑76 (MQ=255)
                                              catGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATg                        >  2:790746/1‑76 (MQ=255)
                                               atGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGa                       >  1:2484446/1‑76 (MQ=255)
                                               atGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGa                       <  2:3798995/76‑1 (MQ=255)
                                                tGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAg                      <  2:4128354/76‑1 (MQ=255)
                                                tGCGCTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTatatcac                      >  2:2696727/1‑73 (MQ=255)
                                                  cGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGcc                    >  1:3627105/1‑76 (MQ=255)
                                                  cGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGcc                    <  2:2284488/76‑1 (MQ=255)
                                                  cGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGcc                    <  2:4420369/76‑1 (MQ=255)
                                                   gTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGccc                   <  1:4826310/76‑1 (MQ=255)
                                                   gTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGccc                   >  2:391209/1‑76 (MQ=255)
                                                   gTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGccc                   <  1:100192/76‑1 (MQ=255)
                                                     tCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGt                 <  1:3166734/76‑1 (MQ=255)
                                                      cAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATTGTTATATGAGCCCGtt                <  1:1470350/76‑1 (MQ=255)
                                                      cAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGtt                <  1:4543107/76‑1 (MQ=255)
                                                       aaTGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGTGCCCGTTc               >  2:84527/1‑76 (MQ=255)
                                                       aaTGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTc               >  1:1944460/1‑76 (MQ=255)
                                                         tGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCtt             >  1:522409/1‑76 (MQ=255)
                                                           ggCGAAAGGTTGTGCCTCGAACCGGTGACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTg           <  2:3935111/76‑1 (MQ=255)
                                                            gCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg          <  1:2849370/76‑1 (MQ=255)
                                                            gCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg          <  1:4771758/76‑1 (MQ=255)
                                                            gCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg          <  2:1066334/76‑1 (MQ=255)
                                                              gAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCt        >  1:103889/1‑76 (MQ=255)
                                                                aaTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAt      <  2:522239/76‑1 (MQ=255)
                                                                 aTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAtc     >  2:1377243/1‑76 (MQ=255)
                                                                  tgttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAtct    >  1:3282435/1‑76 (MQ=255)
                                                                   gttgttCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTAtctc   <  1:2472870/76‑1 (MQ=255)
                                                                    ttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa  >  1:1334592/1‑76 (MQ=255)
                                                                    ttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa  >  2:3640558/1‑76 (MQ=255)
                                                                    ttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa  <  2:1496545/76‑1 (MQ=255)
                                                                    ttgtTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCa  <  1:663689/76‑1 (MQ=255)
                                                                    ttgtTCCACGAACCTGTTACCGGTCCGTTCGCCGTTTATGCGCCTCTTTATATGAGCCCATTCTTCGTCCATCTCa  <  1:2536734/76‑1 (MQ=255)
                                                                        |                                                                       
CCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTCCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCA  >  NC_000913/3978630‑3978773

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