Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,141,886 A→G 7.8% intergenic (+264/+10) yegH → / ← asmA fused predicted membrane proteins/predicted assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,141,8860AG 7.8% 131.8 / 3.9 51intergenic (+264/+10)yegH/asmAfused predicted membrane proteins/predicted assembly protein
Reads supporting (aligned to +/- strand):  ref base A (28/19);  new base G (2/2);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

ATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGC  >  W3110S.gb/2141822‑2141953
                                                                |                                                                   
atgatgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCGCCTgg                                                                <  1:895592/70‑1 (MQ=255)
atgatgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGt                                                               <  1:2184519/71‑1 (MQ=255)
atgatgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGGCACCTGGt                                                               <  1:1352019/71‑1 (MQ=255)
           aCGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtct                                                    >  1:3979057/1‑71 (MQ=255)
             gAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtcttc                                                  <  1:3473157/71‑1 (MQ=255)
             gAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtcttc                                                  >  1:2049052/1‑71 (MQ=255)
              aaCGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCAcc                                                                   >  1:499611/1‑53 (MQ=255)
                     gATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTa                                                              <  1:1088279/51‑1 (MQ=255)
                     gATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtc                                                     <  1:1418013/60‑1 (MQ=255)
                      aTGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                       >  1:789958/1‑57 (MQ=255)
                       tGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAAct                                        >  1:3825824/1‑71 (MQ=255)
                       tGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAAct                                        >  1:4647730/1‑71 (MQ=255)
                           aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGTTCACCTGGTAACCTACATCTTCTCCAGCAActtctt                                    <  1:1430745/71‑1 (MQ=255)
                           aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCGcc                                                                   <  1:3750870/40‑1 (MQ=255)
                           aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtctt                                    <  1:263374/71‑1 (MQ=255)
                           aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAAct                                        <  1:371273/67‑1 (MQ=255)
                           aCGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAAct                                        <  1:1715349/67‑1 (MQ=255)
                               tttACGCGTCTTATCAGGCCTACAGGACTGGTCAcc                                                                   <  1:921404/36‑1 (MQ=255)
                                  aCGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCaa                                          >  1:1435347/1‑58 (MQ=255)
                                     cgTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTg                          <  1:3571097/71‑1 (MQ=255)
                                      gTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTcac                                 <  1:1219086/63‑1 (MQ=255)
                                       tCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtctt                                    >  1:2774683/1‑59 (MQ=255)
                                       tCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGc                         >  1:4476337/1‑70 (MQ=255)
                                        cTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtctt                                                   >  1:2347709/1‑43 (MQ=255)
                                        cTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtctt                                                   >  1:4246312/1‑43 (MQ=255)
                                        cTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCGACTTCTTCACAt                               >  1:4409166/1‑63 (MQ=255)
                                         ttATCAGGCCTACAGGACTGGTCGCCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCAt                      >  1:2286517/1‑71 (MQ=255)
                                         ttATCAGGCCTACAGGACTGGTCACCTGGTAACCTACa                                                       >  1:3917975/1‑38 (MQ=255)
                                         ttATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCAt                      >  1:1152815/1‑71 (MQ=255)
                                           aTCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTa                    >  1:1929748/1‑71 (MQ=255)
                                              aGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTa                    >  1:2622733/1‑68 (MQ=255)
                                              aGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACgg                 >  1:3709071/1‑71 (MQ=255)
                                                gCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTa                    >  1:2470992/1‑66 (MQ=255)
                                                   tACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCAt                      <  1:3417647/61‑1 (MQ=255)
                                                   tACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGa               <  1:455161/68‑1 (MQ=255)
                                                   tACAGGACTGGACGCCTGGTAACCTACATCTTCTCCAGCAACTTCTTcac                                 >  1:4627116/1‑50 (MQ=39)
                                                    aCAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATc            <  1:3846821/70‑1 (MQ=255)
                                                     cAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtctt                                    >  1:3642676/1‑45 (MQ=255)
                                                     cAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActtctt                                    >  1:892980/1‑45 (MQ=255)
                                                     cAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActt                                       >  1:3639657/1‑42 (MQ=255)
                                                     cAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt          >  1:965830/1‑71 (MQ=255)
                                                     cAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt          >  1:4635117/1‑71 (MQ=255)
                                                      aGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAActt                                       >  1:3271499/1‑41 (MQ=255)
                                                      aGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCa                       >  1:2762017/1‑57 (MQ=255)
                                                       ggACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACttcttc                                   <  1:2947905/44‑1 (MQ=255)
                                                       ggACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACttcttc                                   <  1:2204866/44‑1 (MQ=255)
                                                        gACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt          <  1:3608287/68‑1 (MQ=255)
                                                        gACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTAt       >  1:2819146/1‑71 (MQ=255)
                                                          cTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTc     <  1:1555211/71‑1 (MQ=255)
                                                             gTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt          >  1:3777843/1‑63 (MQ=255)
                                                             gTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGc  >  1:4638605/1‑71 (MQ=255)
                                                                aCCTGGTAACCTACATCTTCTCCAGCAACTTCTTAACATCTTTGCCATTACGGGAATCTTTATTc     >  1:2831722/1‑65 (MQ=255)
                                                                |                                                                   
ATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGC  >  W3110S.gb/2141822‑2141953

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: