Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,927,471 T→A 9.3% P106P (CCT→CCA yobB → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,927,4710TA 9.3% 140.4 / 6.8 54P106P (CCT→CCAyobBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/21);  new base A (3/2);  total (31/23)
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.36e-01

TCCCGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAA  >  W3110S.gb/1927405‑1927540
                                                                  |                                                                     
tCCCGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGa                                                                   >  1:3734067/1‑71 (MQ=255)
 cccGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGa                                                                   >  1:1610516/1‑70 (MQ=255)
 cccGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGAt                                                                  >  1:1442778/1‑71 (MQ=255)
  ccGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATc                                                                 >  1:3767696/1‑71 (MQ=255)
  ccGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATc                                                                 >  1:2048796/1‑71 (MQ=255)
   cGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATc                                                                 >  1:1633867/1‑70 (MQ=255)
    gTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTa                                                               >  1:2141317/1‑71 (MQ=255)
       gAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGAt                                                                  >  1:1359812/1‑65 (MQ=255)
       gAATATAACGATCGCTTTATACGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCa                                                            >  1:2866111/1‑71 (MQ=255)
        aaTATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCAt                                                           <  1:3122141/71‑1 (MQ=255)
        aaTATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCAt                                                           <  1:2888408/71‑1 (MQ=255)
           ataACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCAt                                                           >  1:3825970/1‑68 (MQ=255)
           ataACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCAt                                                           >  1:3622702/1‑68 (MQ=255)
              aCGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCt                                                                       >  1:2384669/1‑53 (MQ=255)
                gATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCa                                                   >  1:922817/1‑71 (MQ=255)
                 aTCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATc                                                          <  1:258231/63‑1 (MQ=255)
                 aTCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCAt                                                  <  1:615304/71‑1 (MQ=255)
                 aTCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCAt                                                  <  1:528390/71‑1 (MQ=255)
                  tCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTggg                                                                    <  1:232377/52‑1 (MQ=255)
                  tCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTggg                                                                    <  1:1871475/52‑1 (MQ=255)
                    gCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATg                                                 >  1:2232979/1‑69 (MQ=255)
                       ttATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGc                                             <  1:3365099/70‑1 (MQ=255)
                        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGGTCTACCATCAGAGTCATGGTGc                                             <  1:3674819/69‑1 (MQ=255)
                        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGc                                             <  1:2086083/69‑1 (MQ=255)
                        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCAt                                           >  1:1136390/1‑71 (MQ=255)
                        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCAt                                           >  1:847286/1‑71 (MQ=255)
                        tATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCAt                                           >  1:1751288/1‑71 (MQ=255)
                                    tGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCcgtc                               >  1:2837209/1‑71 (MQ=255)
                                    tGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCcgtc                               >  1:2841839/1‑71 (MQ=255)
                                         tgtTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGt                                                     >  1:1810873/1‑44 (MQ=255)
                                             cGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGaa                      >  1:2127498/1‑71 (MQ=255)
                                             cGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGaa                      >  1:1602468/1‑71 (MQ=255)
                                             cGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGaa                      >  1:330852/1‑71 (MQ=255)
                                             cGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGa                       <  1:3077828/70‑1 (MQ=255)
                                                 ccGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTg                                      <  1:1397984/51‑1 (MQ=255)
                                                   gTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCa                   >  1:1302973/1‑68 (MQ=255)
                                                   gTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCa                <  1:501694/71‑1 (MQ=255)
                                                    tGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCa                <  1:2961525/70‑1 (MQ=255)
                                                    tGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCAc               <  1:1917853/71‑1 (MQ=255)
                                                      gCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACgg             <  1:2412206/71‑1 (MQ=255)
                                                      gCGAAAAACGCCTGGGATCTACAATCAGAGTCATGGTGCATGTCTGGGACGTCGTTCCAGAACCATCACgg             >  1:874922/1‑71 (MQ=255)
                                                        gAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACggtg           >  1:50938/1‑71 (MQ=255)
                                                        gAAAAACGCCTGGGATCTACCATCAGAGTCATGGAGCATGTCTGGGCCGTCGTTCCAGAACCATCACggtg           >  1:283765/1‑71 (MQ=255)
                                                        gAAAAACGCCAGGGATCTACCATCAGAGTTATGGt                                               <  1:851643/35‑1 (MQ=37)
                                                        gAAAAACGCCAGGGATCTACCATCAGAGTTATGGt                                               <  1:13170/35‑1 (MQ=37)
                                                          aaaaCGCCAGGGATCTACCATCAGAGTTATGGTGCATGTCTGGGCCGTCGTTcc                          >  1:613548/1‑54 (MQ=255)
                                                          aaaaCGCCAGGGATCTACCATCAGAGTTATGGTGCATGTCTGGGCCGTCGTTcc                          >  1:2270927/1‑54 (MQ=255)
                                                          aaaaCGCCAGGGATCTACCATCAGAGTTATGGTGCATGTCTGGGCCGTCGTTcc                          >  1:836011/1‑54 (MQ=255)
                                                           aaaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTg                                      >  1:999990/1‑41 (MQ=255)
                                                           aaaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTc        <  1:1872771/71‑1 (MQ=255)
                                                            aaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTc        >  1:3579386/1‑70 (MQ=255)
                                                            aaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCg       <  1:444366/71‑1 (MQ=255)
                                                            aaCGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCg       <  1:1085933/71‑1 (MQ=255)
                                                              cGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGAt     <  1:2363272/71‑1 (MQ=255)
                                                                 cTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGaa  <  1:3692656/71‑1 (MQ=255)
                                                                  tGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGaa  >  1:1990153/1‑70 (MQ=255)
                                                                  |                                                                     
TCCCGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCTGGGATCTACCATCAGAGTCATGGTGCATGTCTGGGCCGTCGTTCCAGAACCATCACGGTGGTCGATGAA  >  W3110S.gb/1927405‑1927540

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: