Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 57,363 T→C 9.1% intergenic (‑254/‑1) imp ← / → djlA exported protein required for envelope biosynthesis and integrity/DnaJ‑like protein, membrane anchored

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb57,3630TC 9.1% 143.0 / 6.5 55intergenic (‑254/‑1)imp/djlAexported protein required for envelope biosynthesis and integrity/DnaJ‑like protein, membrane anchored
Reads supporting (aligned to +/- strand):  ref base T (16/34);  new base C (2/3);  total (18/37)
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.37e-01

AAACGGTAAAAGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCTTTTGGGGCGT  >  W3110S.gb/57297‑57431
                                                                  |                                                                    
aaaCGGTAAAAGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGc                                                                  <  1:3438750/71‑1 (MQ=255)
        aaaGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTgg                                                          >  1:1309475/1‑71 (MQ=255)
              ccTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGaaaa                                                    <  1:2780163/71‑1 (MQ=255)
                 ttATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAAt                                                  >  1:1726293/1‑70 (MQ=255)
                 ttATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAAt                                                  >  1:2910655/1‑70 (MQ=255)
                 ttATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAAt                                                  >  1:1338035/1‑70 (MQ=255)
                  tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATc                                                 <  1:1571185/70‑1 (MQ=255)
                  tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCa                                                >  1:2128308/1‑71 (MQ=255)
                  tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCa                                                >  1:2183788/1‑71 (MQ=255)
                      ttGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGt                                                               <  1:2340471/52‑1 (MQ=255)
                      ttGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTgg                                            <  1:1752307/71‑1 (MQ=255)
                      ttGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTgg                                            <  1:668312/71‑1 (MQ=255)
                        gtgGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATACATGCAGTATTGGGGAAAAATCAt                                               >  1:686473/1‑66 (MQ=255)
                         tgGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcg                                          <  1:594033/70‑1 (MQ=255)
                           ggAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGaaa                                                     <  1:73921/57‑1 (MQ=255)
                           ggAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGaaa                                                     <  1:3563137/57‑1 (MQ=255)
                           ggAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGaaa                                                     <  1:122036/57‑1 (MQ=255)
                           ggAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGaaa                                                     <  1:1010578/57‑1 (MQ=255)
                            gagaTAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtggc                                      <  1:2051480/71‑1 (MQ=255)
                                 aGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTa                                                              >  1:362249/1‑42 (MQ=255)
                                 aGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTa                                                              >  1:3245462/1‑42 (MQ=255)
                                 aGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTgg                                                          >  1:3552853/1‑46 (MQ=255)
                                  gCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtggccgtggc                                <  1:184431/71‑1 (MQ=255)
                                      tGATATCCGTGTGTCGATTTGGGGAATATATGCAg                                                                <  1:2670324/35‑1 (MQ=255)
                                      tGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTa                            >  1:804546/1‑71 (MQ=255)
                                        atatCCGTGTGTCGATTTGGGGAATATATGCAGTa                                                              >  1:2414714/1‑35 (MQ=255)
                                        atatCCGTGTGTCGATTTGGGGAATATATGCAGTa                                                              >  1:1319023/1‑35 (MQ=255)
                                        atatCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTAc                           >  1:420011/1‑70 (MQ=255)
                                         tatCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtg                                        <  1:160975/56‑1 (MQ=255)
                                         tatCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtg                                        <  1:831107/56‑1 (MQ=255)
                                         tatCCGTGTGTCGATTTGGGGAATACATGCAGTATTGGGGAAAAAt                                                  <  1:4150928/46‑1 (MQ=255)
                                           tCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGaaaa                                                    <  1:4569781/42‑1 (MQ=255)
                                           tCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCt                              <  1:2879381/64‑1 (MQ=255)
                                           tCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACt                          <  1:4398908/68‑1 (MQ=255)
                                               tgtgtCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggc                   >  1:1076771/1‑71 (MQ=255)
                                                 tgtCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcg                  <  1:1123518/70‑1 (MQ=255)
                                                 tgtCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcg                  <  1:262185/70‑1 (MQ=255)
                                                  gtCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGa                        <  1:409548/63‑1 (MQ=255)
                                                   tCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTgg                                            <  1:2988195/42‑1 (MQ=255)
                                                   tCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACt                          <  1:457636/60‑1 (MQ=255)
                                                   tCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGa                        <  1:2332731/62‑1 (MQ=255)
                                                   tCGATTTGGGGAATACATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACt                          <  1:2669417/60‑1 (MQ=255)
                                                   tCGATTTGGGGAATACATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcg               >  1:154341/1‑71 (MQ=255)
                                                     gATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcggc             <  1:437564/71‑1 (MQ=255)
                                                     gATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcgg              >  1:2833206/1‑70 (MQ=255)
                                                      aTTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtggc                                      <  1:3025540/45‑1 (MQ=255)
                                                       tttGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCtt           <  1:1754725/71‑1 (MQ=255)
                                                       tttGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCtt           <  1:978338/71‑1 (MQ=255)
                                                           gggAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCTTTTg        <  1:1462473/70‑1 (MQ=255)
                                                           gggAATACATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCTTTTg        <  1:3812498/70‑1 (MQ=255)
                                                            ggAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCTTTTggg      <  1:2584217/71‑1 (MQ=255)
                                                             taaTATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCt                              <  1:1926450/45‑1 (MQ=255)
                                                             gAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACt                          <  1:2503101/50‑1 (MQ=255)
                                                              aaTATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCt                              <  1:3413714/45‑1 (MQ=255)
                                                              aaTATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGGGCTTTTGGGGCg   >  1:1934138/1‑71 (MQ=255)
                                                                tatatGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCTTTTGGGGCGt  <  1:3944590/71‑1 (MQ=255)
                                                                  |                                                                    
AAACGGTAAAAGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCTTTTGGGGCGT  >  W3110S.gb/57297‑57431

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: