Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,515,054 A→G 9.1% N41D (AAC→GAC)  hslU → molecular chaperone and ATPase component of HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,515,0540AG 9.1% 113.8 / 4.2 44N41D (AAC→GAC) hslUmolecular chaperone and ATPase component of HslUV protease
Reads supporting (aligned to +/- strand):  ref base A (23/17);  new base G (2/2);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

ATCGGCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTAAAAC  >  W3110S.gb/3514985‑3515124
                                                                     |                                                                      
atcGGCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCaa                                                                       <  1:3360862/71‑1 (MQ=255)
atcGGCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCaa                                                                       >  1:953748/1‑71 (MQ=255)
   ggCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACg                                                                     >  1:3473388/1‑70 (MQ=255)
   ggCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACg                                                                     <  1:2302520/70‑1 (MQ=255)
    gCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACg                                                                     <  1:2617679/69‑1 (MQ=255)
     ccAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAag                                                                  <  1:3716935/71‑1 (MQ=255)
            aaCGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAaga                                                                 <  1:48926/65‑1 (MQ=255)
              cGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGc                                                               >  1:1473086/1‑65 (MQ=255)
               gCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCa                                                        >  1:2615924/1‑71 (MQ=255)
                 cAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCGACGAAGAGCTGCGCCATg                                                      <  1:3148397/71‑1 (MQ=255)
                 cAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCGACGAAGAGCTGCGCCATg                                                      <  1:4492960/71‑1 (MQ=255)
                 cAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCAt                                                       >  1:1291851/1‑70 (MQ=255)
                   aGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTgcgc                                                          <  1:4439478/65‑1 (MQ=255)
                         cTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAagag                                                                >  1:1601963/1‑53 (MQ=255)
                         cTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTgcg                                                           >  1:1441375/1‑58 (MQ=255)
                         cTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGAcc                                              >  1:3467744/1‑71 (MQ=255)
                         cTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGAcc                                              >  1:1984207/1‑71 (MQ=255)
                           gtgGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGaa                                                                   >  1:2401495/1‑48 (MQ=255)
                                   tGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGa                                                     >  1:769613/1‑54 (MQ=255)
                                   tGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGa                                                     >  1:3319018/1‑54 (MQ=255)
                                    gCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGaa                                                                   >  1:4184079/1‑39 (MQ=255)
                                    gCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGaa                                                                   >  1:1030120/1‑39 (MQ=255)
                                    gCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAatat                                   >  1:2140393/1‑71 (MQ=255)
                                    gCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAatat                                   >  1:3508451/1‑71 (MQ=255)
                                    gCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAatat                                   >  1:3813454/1‑71 (MQ=255)
                                      tctGCGTAACCGCTGGCGTCGCATGCAGCTCGACGAAGAGCTGCGCCATGaa                                                    >  1:3880025/1‑52 (MQ=255)
                                      tctGCGTAACCGCTGGCGTCGCATGCAGCTCGACGAAGAGCTGCGCCATGaa                                                    >  1:4172248/1‑52 (MQ=255)
                                              aCCGCTGGCGTCGCATGCAGCGCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTgatgat                          <  1:3771959/70‑1 (MQ=255)
                                                cGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGaaaaa                                      <  1:981026/56‑1 (MQ=255)
                                                  cTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGc                      <  1:2944740/70‑1 (MQ=255)
                                                     gCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTgatgat                          >  1:2967554/1‑63 (MQ=255)
                                                         cGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCg               >  1:3905993/1‑70 (MQ=255)
                                                          gCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCgg              <  1:2766729/70‑1 (MQ=255)
                                                             tGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTc          <  1:1869873/71‑1 (MQ=255)
                                                               cAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCgg        >  1:2177749/1‑71 (MQ=255)
                                                                aGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCgg        >  1:2876849/1‑70 (MQ=255)
                                                                aGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCg         >  1:4072141/1‑69 (MQ=255)
                                                                 gCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGAcc                <  1:1859729/61‑1 (MQ=255)
                                                                  cTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCtgatga                           <  1:2500801/49‑1 (MQ=255)
                                                                  cTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGgt             <  1:3723008/63‑1 (MQ=255)
                                                                  cTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTaa     <  1:3825903/71‑1 (MQ=255)
                                                                  cTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTaa     <  1:2210445/71‑1 (MQ=255)
                                                                   tCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCtgatg                            <  1:2983132/47‑1 (MQ=255)
                                                                   tCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTaa     <  1:2147638/70‑1 (MQ=255)
                                                                    cAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCtga                              >  1:2035704/1‑44 (MQ=255)
                                                                    cAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCgt                       >  1:3558950/1‑50 (MQ=255)
                                                                     aaCGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTAAAAc  >  1:2590807/1‑71 (MQ=255)
                                                                     |                                                                      
ATCGGCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTAAAAC  >  W3110S.gb/3514985‑3515124

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