Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,514,320 (C)7→6 intergenic (+19/‑74) ftsN → / → hslV essential cell division protein/peptidase component of the HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,514,3140C.97.8% 176.7 / ‑2.2 46intergenic (+13/‑80)ftsN/hslVessential cell division protein/peptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (15/30);  minor base A (0/1);  total (15/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAATCGGTTGAAGATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAATCCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTGTA  >  W3110S.gb/3514244‑3514371
                                                                       |                                                         
cAATCGGTTGAAGATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAATc                                                           >  1:1734500/1‑71 (MQ=255)
         gAAGATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAATccccccatc                                                   <  1:4456011/70‑4 (MQ=255)
            gATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCtata                                              <  1:784447/71‑1 (MQ=255)
            gATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCtata                                              <  1:2179/71‑1 (MQ=255)
                gCGGGTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTg                                          <  1:3982880/71‑1 (MQ=255)
                   gggcATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:3696864/68‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:1032886/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:3112049/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:278236/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:2754947/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:4588372/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:1641777/71‑1 (MQ=255)
                    gTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:696635/70‑1 (MQ=255)
                     tCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTa                                     >  1:166611/1‑71 (MQ=255)
                     tCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTa                                     >  1:2786570/1‑71 (MQ=255)
                       aTACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATg                                   <  1:2845801/71‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCACAAT‑CCCCCCATCTATAATTGCAt                                       <  1:741561/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:2447453/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:2747770/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:746853/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:3493145/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:3388195/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAA‑ACCCCCAATCTATAATTGCAt                                       <  1:2316514/66‑1 (MQ=255)
                         aCAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTa                                     >  1:1272374/1‑67 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       >  1:2289944/1‑65 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAATc                                                           >  1:1665052/1‑45 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTa                                     >  1:1533990/1‑66 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTa                                     >  1:2067131/1‑66 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGccc                                >  1:3294296/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGccc                                >  1:870771/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGccc                                >  1:4049857/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGccc                                >  1:4061674/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGccc                                >  1:4287921/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGccc                                >  1:4371852/1‑71 (MQ=255)
                                     cGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGCCCCGTACTTTTg                      >  1:4438020/1‑70 (MQ=255)
                                      ggCTCGCCGCC‑GGGGGTTGAAACCCTCAAAAT‑CCCCCCATCtat                                               <  1:4601180/44‑1 (MQ=255)
                                               gc‑gGGGGGTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:3408220/42‑1 (MQ=255)
                                               cc‑gGGGGTTGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTg           <  1:3202348/71‑1 (MQ=255)
                                                       ttGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:422583/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:3853145/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:159851/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:703096/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCAt                                       <  1:3471757/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTg    <  1:25707/71‑1 (MQ=255)
                                                       ttGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTg    <  1:1580043/71‑1 (MQ=255)
                                                       ttGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTg    <  1:1291271/71‑1 (MQ=255)
                                                       ttGAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTg    <  1:947989/71‑1 (MQ=255)
                                                         gAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTGTa  >  1:450422/1‑71 (MQ=255)
                                                         gAAACCCTCAAAAT‑CCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTGTa  >  1:4541448/1‑71 (MQ=255)
                                                                       |                                                         
CAATCGGTTGAAGATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCC‑GGGGGTTGAAACCCTCAAAATCCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTGTA  >  W3110S.gb/3514244‑3514371

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: ‑