Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 57,363 T→C 10.2% 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,3630TC10.2% 156.5 / 9.5 59intergenic (‑254/‑1)imp/djlAexported protein required for envelope biosynthesis and integrity/DnaJ‑like protein, membrane anchored
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base C (3/3);  total (33/26)
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.85e-01

AAAACGGTAAAAGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCTTT  >  W3110S.gb/57296‑57423
                                                                   |                                                            
aaaaCGGTAAAAGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATg                                                           <  1:1873238/71‑1 (MQ=255)
   aCGGTAAAAGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAg                                                        <  1:399178/71‑1 (MQ=255)
   aCGGTAAAAGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAg                                                        <  1:3636459/71‑1 (MQ=255)
    cGGTAAAAGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAg                                                        <  1:1281775/70‑1 (MQ=255)
         aaaGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAatat                                                              <  1:347897/59‑1 (MQ=255)
         aaaGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCa                                                         >  1:4022194/1‑64 (MQ=255)
         aaaGCACCTTTATATTGTGCGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAatat                                                              <  1:629872/59‑1 (MQ=255)
                 tttATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGa                                               <  1:1712908/66‑1 (MQ=255)
                   tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTggg                                                 >  1:3331658/1‑62 (MQ=255)
                   tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAAt                                          <  1:4020952/69‑1 (MQ=255)
                   tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAAt                                          <  1:1439848/69‑1 (MQ=255)
                   tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCa                                        >  1:2456532/1‑71 (MQ=255)
                   tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCa                                        >  1:3092848/1‑71 (MQ=255)
                   tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATACATGCAGTATTGGGGAAAAAt                                          <  1:1143983/69‑1 (MQ=255)
                   tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATACATGCAGTATTGGGGAAAAAt                                          <  1:1918685/69‑1 (MQ=255)
                   tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATACATGCAGTATTGGGGAAAAATCa                                        >  1:3700373/1‑71 (MQ=255)
                   tataTTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATACATGCAGTATTGGGGAAAAATCa                                        >  1:2377356/1‑71 (MQ=255)
                    atatTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCa                                        <  1:1311380/70‑1 (MQ=255)
                     tatTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCa                                        <  1:2557594/69‑1 (MQ=255)
                       ttGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTgg                                    >  1:1905105/1‑71 (MQ=255)
                              agaTAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGaaa                                             <  1:728995/55‑1 (MQ=255)
                               gaTAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtggccg                            >  1:3326584/1‑71 (MQ=255)
                                aTAGCCCTGATATCCGTGTGTCGATTTGGGGAATACATGCAGTATTGGGGAAAAATCa                                        <  1:3346507/58‑1 (MQ=255)
                                  aGCCCTGATATCCGTGTGTCGATTTGGGGAATACATGCAGTATTGGGGAAAAATCATTGGcgtggccgtg                          >  1:2395712/1‑70 (MQ=255)
                                   gCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtggccgtggc                        >  1:3000640/1‑71 (MQ=255)
                                   gCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtggccgtggc                        >  1:1277951/1‑71 (MQ=255)
                                     gccgATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGGGGCCGTGGCCt                      <  1:1924017/68‑1 (MQ=255)
                                     ccTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCt                      <  1:2209194/71‑1 (MQ=255)
                                      cTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCtt                     <  1:2365613/71‑1 (MQ=255)
                                       tGATATCCGTGTGTCGATTTGGGGAATATATGCAg                                                        <  1:1483049/35‑1 (MQ=255)
                                        gATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTAc                   >  1:1118297/1‑71 (MQ=255)
                                        gATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTAc                   >  1:3476254/1‑71 (MQ=255)
                                         atatCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCa                                        >  1:2561839/1‑49 (MQ=255)
                                         atatCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCa                                        >  1:511078/1‑49 (MQ=255)
                                         atatCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCt                      <  1:3282980/67‑1 (MQ=255)
                                           atCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTa                    >  1:3829890/1‑67 (MQ=255)
                                           atCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTa                    >  1:3578105/1‑67 (MQ=255)
                                           atCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGa                >  1:547217/1‑71 (MQ=255)
                                            tCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAAt                                          <  1:1294284/44‑1 (MQ=255)
                                            tCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAAt                                          <  1:3293522/44‑1 (MQ=255)
                                            tCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtggccgt                           <  1:1986774/59‑1 (MQ=255)
                                            tCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTa                    >  1:1208191/1‑66 (MQ=255)
                                            tCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGAt               >  1:2687686/1‑71 (MQ=255)
                                             ccGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCAtt                                      >  1:3045318/1‑47 (MQ=255)
                                             ccGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCAtt                                      >  1:1709695/1‑47 (MQ=255)
                                             ccGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCAtt                                      >  1:3027325/1‑47 (MQ=255)
                                             ccGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCAtt                                      >  1:138823/1‑47 (MQ=255)
                                             ccGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCAtt                                      >  1:3839526/1‑47 (MQ=255)
                                              cGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATg              <  1:3568291/70‑1 (MQ=255)
                                               gtgtgtCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGa                >  1:3685416/1‑67 (MQ=255)
                                               gtgtgtCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGa                >  1:2229330/1‑67 (MQ=255)
                                                    tCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcg       >  1:1283200/1‑71 (MQ=255)
                                                      gATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcggc     >  1:949004/1‑71 (MQ=255)
                                                      gATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcggc     >  1:3920504/1‑71 (MQ=255)
                                                      gATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcggc     >  1:3252109/1‑71 (MQ=255)
                                                      gATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcgg      >  1:2554778/1‑70 (MQ=255)
                                                      gATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcg       >  1:3078860/1‑69 (MQ=255)
                                                        tttGGGGAATATATGCAGTATTGGGGAAAAATCATTGGcgtggcc                             <  1:1834545/45‑1 (MQ=255)
                                                        tttGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCtt   <  1:3776083/71‑1 (MQ=255)
                                                          tGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCttt  <  1:3412707/70‑1 (MQ=255)
                                                          tGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCttt  <  1:2994816/70‑1 (MQ=255)
                                                                atatatGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGggcggcg       >  1:1218646/1‑59 (MQ=255)
                                                                 tatatGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACt                  <  1:690002/47‑1 (MQ=255)
                                                                   |                                                            
AAAACGGTAAAAGCACCTTTATATTGTGGGAGATAGCCCTGATATCCGTGTGTCGATTTGGGGAATATATGCAGTATTGGGGAAAAATCATTGGCGTGGCCGTGGCCTTACTGATGGGCGGCGGCTTT  >  W3110S.gb/57296‑57423

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: