Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,515,054 A→G 15.0% 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,0540AG15.0% 96.6 / 10.6 40N41D (AAC→GAC) hslUmolecular chaperone and ATPase component of HslUV protease
Reads supporting (aligned to +/- strand):  ref base A (13/21);  new base G (4/2);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GGCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTAA  >  W3110S.gb/3514988‑3515121
                                                                  |                                                                   
ggCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCGACg                                                                  >  1:2651264/1‑70 (MQ=255)
ggCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACg                                                                  <  1:1317870/70‑1 (MQ=255)
 gCCATGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACg                                                                  <  1:1091449/69‑1 (MQ=255)
          aCGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTgcgc                                                       <  1:770299/71‑1 (MQ=255)
           cGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTgcgc                                                       >  1:3380716/1‑70 (MQ=255)
              cAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATg                                                   <  1:2071697/71‑1 (MQ=255)
              cAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATg                                                   <  1:207710/71‑1 (MQ=255)
              cAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATg                                                   <  1:3139657/71‑1 (MQ=255)
                  cGTTCTGTGGCGTTTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGt                                               <  1:3353823/71‑1 (MQ=255)
                         tgGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGaa                                                                <  1:1844971/47‑1 (MQ=255)
                         tgGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGaa                                                                <  1:2536225/47‑1 (MQ=255)
                               ttGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAat                                  <  1:3571320/71‑1 (MQ=255)
                               ttGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAat                                  <  1:1290226/71‑1 (MQ=255)
                                tGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGaa                                                 >  1:907330/1‑55 (MQ=255)
                                tGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGaa                                                 >  1:402369/1‑55 (MQ=255)
                                tGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGaa                                                 >  1:2507224/1‑55 (MQ=255)
                                tGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAata                                 >  1:3190206/1‑71 (MQ=255)
                                 gCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCGACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAatat                                >  1:3520738/1‑71 (MQ=255)
                                 gCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCGACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAata                                 >  1:2202517/1‑70 (MQ=255)
                                 gCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAata                                 <  1:596841/70‑1 (MQ=255)
                                 gCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAata                                 <  1:689514/70‑1 (MQ=255)
                                  ctctGCGTAACCGCTGGCGTCGCATGCAGCTCGACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAatat                                <  1:3500702/70‑1 (MQ=255)
                                       cGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCtgat                          >  1:1817525/1‑71 (MQ=255)
                                       cGTAACCGCTGGCGTCGCATACAGCTCGACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCt                             <  1:4009714/68‑1 (MQ=255)
                                         tAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGaaa                                     <  1:3255942/58‑1 (MQ=255)
                                              gCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGc                   >  1:933704/1‑71 (MQ=255)
                                                  gCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGccc                 <  1:308422/69‑1 (MQ=255)
                                                  gCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGccc                 <  1:1943939/69‑1 (MQ=255)
                                                    gTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTgatgat                       >  1:3275411/1‑61 (MQ=255)
                                                    gTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTgatgat                       >  1:2129737/1‑61 (MQ=255)
                                                       gCATGCAGCTCGACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTgatgat                       >  1:2520604/1‑58 (MQ=255)
                                                       gCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGgt          >  1:321011/1‑71 (MQ=255)
                                                       gCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGgt          >  1:1195596/1‑71 (MQ=255)
                                                       gCATGCAGCTCAACGAAGAGCTGAGCCATGAAGTGACCCCGAAAAATATCCTgatgat                       >  1:3973270/1‑58 (MQ=255)
                                                            cAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGcc                  >  1:3051892/1‑58 (MQ=255)
                                                               cTCAACGAAGAGCTGCGCCATGAAGTGACCCCGaaaaa                                   <  1:3789991/38‑1 (MQ=255)
                                                               cTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTaa  <  1:3129800/71‑1 (MQ=255)
                                                               cTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTaa  <  1:3863521/71‑1 (MQ=255)
                                                               cTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGACCGACCGGTGTCGGTaa  <  1:825513/71‑1 (MQ=255)
                                                                tCAACGAAGAGCTGCGCCATGAAGTGACCCCGaaaaa                                   <  1:2289602/37‑1 (MQ=255)
                                                                  |                                                                   
GGCCAGGACAACGCCAAGCGTTCTGTGGCGATTGCTCTGCGTAACCGCTGGCGTCGCATGCAGCTCAACGAAGAGCTGCGCCATGAAGTGACCCCGAAAAATATCCTGATGATCGGCCCGACCGGTGTCGGTAA  >  W3110S.gb/3514988‑3515121

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: