Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,750,169 C→T 38.4% T228M (ACG→ATG)  yfjB → NAD kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,750,1690CT38.4% 61.0 / 79.2 73T228M (ACG→ATG) yfjBNAD kinase
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base T (18/10);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

GCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAG  >  W3110S.gb/2750099‑2750231
                                                                      |                                                              
gCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCAc                                                                >  1:869219/1‑71 (MQ=255)
  gATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGa                                                              <  1:1895053/71‑1 (MQ=255)
      aCCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATcc                                                           <  1:863094/70‑1 (MQ=255)
      aCCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCg                                                          >  1:2200770/1‑71 (MQ=255)
       cccTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCg                                                          <  1:1384161/70‑1 (MQ=255)
        ccTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTc                                                        <  1:3319751/71‑1 (MQ=255)
           ggTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGc                                                     >  1:3476571/1‑71 (MQ=255)
            gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATc                                                            >  1:916336/1‑63 (MQ=255)
            gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATc                                                            >  1:1194529/1‑63 (MQ=255)
            gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCg                                                    <  1:20550/71‑1 (MQ=255)
            gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCg                                                    <  1:1679304/71‑1 (MQ=255)
               cccATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGtt                                                  <  1:3525437/70‑1 (MQ=255)
               cccATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGtt                                                  <  1:1184463/70‑1 (MQ=255)
                   tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                             <  1:3498258/71‑1 (MQ=255)
                   tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                             <  1:831150/71‑1 (MQ=255)
                   tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                             >  1:1591403/1‑71 (MQ=255)
                   tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                             >  1:1175133/1‑71 (MQ=255)
                   tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                             >  1:1003463/1‑71 (MQ=255)
                   tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                             <  1:2542206/71‑1 (MQ=255)
                   tGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCg                                             >  1:2128975/1‑71 (MQ=255)
                     ttCCCGTATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGc                                            <  1:1012380/70‑1 (MQ=255)
                     ttCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGc                                            <  1:1146991/70‑1 (MQ=255)
                     ttCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCa                                           >  1:1389293/1‑71 (MQ=255)
                     ttCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCa                                           >  1:3320056/1‑71 (MQ=255)
                     ttCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCa                                           >  1:1077134/1‑71 (MQ=255)
                       cccGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACg                                                               >  1:3762710/1‑49 (MQ=255)
                       cccGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGa                                                              >  1:1219636/1‑50 (MQ=255)
                       cccGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGa                                                              >  1:1643058/1‑50 (MQ=255)
                       cccGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGa                                                              >  1:132392/1‑50 (MQ=255)
                           cATACGTTGTCAGCACGACCACTGGTGATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCg                                     >  1:4081166/1‑71 (MQ=255)
                           cATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCg                                     >  1:3732304/1‑71 (MQ=255)
                           cATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCg                                     >  1:203141/1‑71 (MQ=255)
                           cATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGAAGCATGATCCGTCTGCGTTTTTCGCATCGCCg                                     >  1:585068/1‑71 (MQ=255)
                               cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGa                                                              >  1:3097872/1‑42 (MQ=255)
                               cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGa                                                              >  1:1013588/1‑42 (MQ=255)
                               cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATc                                         >  1:3664687/1‑63 (MQ=255)
                               cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATc                                         >  1:2327123/1‑63 (MQ=255)
                               cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATc                                         >  1:2227409/1‑63 (MQ=255)
                                 ttGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGtttt                                                <  1:371208/54‑1 (MQ=255)
                                 ttGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGtttt                                                <  1:1018951/54‑1 (MQ=255)
                                 ttGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGttttt                                               <  1:2725927/55‑1 (MQ=255)
                                 ttGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGttttt                                               <  1:3506455/55‑1 (MQ=255)
                                 ttGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGttttt                                               <  1:3031744/55‑1 (MQ=255)
                                   gTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAAcaga                              >  1:1863375/1‑67 (MQ=255)
                                    tCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGc                                            <  1:3023131/55‑1 (MQ=255)
                                    tCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGc                                            <  1:3086648/55‑1 (MQ=255)
                                    tCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGc                                            <  1:268414/55‑1 (MQ=255)
                                      aGCACGCCCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTg                           >  1:3168970/1‑70 (MQ=255)
                                         aCGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGaaa                       >  1:69972/1‑71 (MQ=255)
                                                  ggTCATAAAGAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                  >  1:3377458/1‑67 (MQ=255)
                                                  ggTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                  >  1:2583949/1‑67 (MQ=255)
                                                  ggTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                  >  1:3502151/1‑67 (MQ=255)
                                                  ggTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                  >  1:400700/1‑67 (MQ=255)
                                                  ggTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGc               >  1:3691222/1‑70 (MQ=255)
                                                  ggTCATAAACAGCAGCAGCACGATCCGTCTGCGTATTTCGCAACGCCGTAACGa                               >  1:1543483/1‑54 (MQ=255)
                                                   gTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGc               <  1:2693442/69‑1 (MQ=255)
                                                   gTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGc               <  1:1052394/69‑1 (MQ=255)
                                                   gTCATAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGc               <  1:1287933/69‑1 (MQ=255)
                                                   gTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTgg                          >  1:2461374/1‑58 (MQ=255)
                                                   gTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGc               <  1:982186/69‑1 (MQ=255)
                                                   gTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGa             >  1:1510174/1‑71 (MQ=255)
                                                   gTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGa             >  1:365454/1‑71 (MQ=255)
                                                       tAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                  <  1:3793181/62‑1 (MQ=255)
                                                       tAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                  <  1:3797791/62‑1 (MQ=255)
                                                       tAAACAGCAGCAGCACGATCCATCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                  <  1:2999488/62‑1 (MQ=255)
                                                        aaaCAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGc         >  1:3175929/1‑70 (MQ=255)
                                                           cagcagcagCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCa                    >  1:1551685/1‑56 (MQ=255)
                                                           cagcagcagCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCa                    >  1:3626998/1‑56 (MQ=255)
                                                              cagcagcaTGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAg  <  1:2296978/71‑1 (MQ=255)
                                                                gcagcaTGATCCGTCTTCGTTTTTCGCATCGCCGTAACGACCTGGAAATc                     <  1:2507576/50‑1 (MQ=255)
                                                                gcagcaTGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATc                     <  1:1699872/50‑1 (MQ=255)
                                                                  agcaTGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCa       >  1:2902476/1‑62 (MQ=255)
                                                                  agcaTGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCa       >  1:1991622/1‑62 (MQ=255)
                                                                  agcaTGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCa       >  1:3445454/1‑62 (MQ=255)
                                                                      |                                                              
GCGATTACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAG  >  W3110S.gb/2750099‑2750231

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