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

CAATCGGTTGAAGATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTGTA  >  W3110S.gb/3514244‑3514371
                                                                       |                                                         
cAATCGGTTGAAGATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATc                                                           >  1:1734500/1‑71 (MQ=255)
         gAAGATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATccccccatc                                                   <  1:4456011/70‑4 (MQ=255)
            gATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCtata                                              <  1:784447/71‑1 (MQ=255)
            gATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCtata                                              <  1:2179/71‑1 (MQ=255)
                gCGGGTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTg                                          <  1:3982880/71‑1 (MQ=255)
                   gggcATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:3696864/68‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:1032886/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:3112049/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:278236/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:2754947/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:4588372/71‑1 (MQ=255)
                   ggTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:1641777/71‑1 (MQ=255)
                    gTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:696635/70‑1 (MQ=255)
                     tCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTa                                     >  1:166611/1‑71 (MQ=255)
                     tCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTa                                     >  1:2786570/1‑71 (MQ=255)
                       aTACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATg                                   <  1:2845801/71‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCACAATCCCCCCATCTATAATTGCAt                                       <  1:741561/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:2447453/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:2747770/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:746853/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:3493145/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:3388195/66‑1 (MQ=255)
                        tACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAAACCCCCAATCTATAATTGCAt                                       <  1:2316514/66‑1 (MQ=255)
                         aCAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTa                                     >  1:1272374/1‑67 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       >  1:2289944/1‑65 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATc                                                           >  1:1665052/1‑45 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTa                                     >  1:1533990/1‑66 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTa                                     >  1:2067131/1‑66 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGccc                                >  1:3294296/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGccc                                >  1:870771/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGccc                                >  1:4049857/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGccc                                >  1:4061674/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGccc                                >  1:4287921/1‑71 (MQ=255)
                          cAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGccc                                >  1:4371852/1‑71 (MQ=255)
                                     cGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTg                      >  1:4438020/1‑70 (MQ=255)
                                      ggCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCATCtat                                               <  1:4601180/44‑1 (MQ=255)
                                               gcgGGGGGTGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:3408220/42‑1 (MQ=255)
                                               ccgGGGGTTGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTg           <  1:3202348/71‑1 (MQ=255)
                                                       ttGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:422583/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:3853145/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:159851/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:703096/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAATCCCCCCATCTATAATTGCAt                                       <  1:3471757/36‑1 (MQ=38)
                                                       ttGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTg    <  1:25707/71‑1 (MQ=255)
                                                       ttGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTg    <  1:1580043/71‑1 (MQ=255)
                                                       ttGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTg    <  1:1291271/71‑1 (MQ=255)
                                                       ttGAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTg    <  1:947989/71‑1 (MQ=255)
                                                         gAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTGTa  >  1:450422/1‑71 (MQ=255)
                                                         gAAACCCTCAAAATCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTGTa  >  1:4541448/1‑71 (MQ=255)
                                                                       |                                                         
CAATCGGTTGAAGATGGCGGGTCATACAAACTGCATTCGGCTCGCCGCCGGGGGTTGAAACCCTCAAAATCCCCCCCATCTATAATTGCATTATGCCCCGTACTTTTGTACGGGGTTTGTACTCTGTA  >  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: