Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,413,273 G→T 23.1% R130R (CGC→CGA ydaO ← predicted C32 tRNA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,413,2730GT23.1% 145.0 / 22.1 56R130R (CGC→CGAydaOpredicted C32 tRNA thiolase
Reads supporting (aligned to +/- strand):  ref base G (29/14);  new base T (10/3);  total (39/17)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCAC  >  W3110S.gb/1413203‑1413332
                                                                      |                                                           
tatCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcg                                                             >  1:1729536/1‑71 (MQ=255)
tatCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcg                                                             >  1:1450464/1‑71 (MQ=255)
  tCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcgc                                                            <  1:3855314/70‑1 (MQ=255)
  tCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcgc                                                            <  1:2587062/70‑1 (MQ=255)
  tCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACcgc                                                            <  1:1549009/70‑1 (MQ=255)
   cGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCTcga                                                          >  1:3733934/1‑71 (MQ=255)
   cGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATAACTcga                                                          >  1:780204/1‑71 (MQ=255)
    gTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCAAGTTCCGTTGCGGTACGATAAAGGATACCGcgac                                                         >  1:2726938/1‑71 (MQ=255)
     tCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcgacg                                                        <  1:1187881/71‑1 (MQ=255)
       aCGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcgacg                                                        >  1:632765/1‑69 (MQ=255)
       aCGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcgacg                                                        >  1:3542533/1‑69 (MQ=255)
            ggTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGcga                                                          >  1:629754/1‑62 (MQ=255)
            ggTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGc                                                   >  1:1875447/1‑69 (MQ=255)
            ggTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGc                                                   >  1:3914789/1‑69 (MQ=255)
             gTGACCCAACGCGATCTTCGTCGTCCCCAGTTCCGTTGCGGTACGATAAAGGATACCTCGACGAAGGCgag                                                <  1:2666422/71‑1 (MQ=255)
             gTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGc                                                   >  1:854874/1‑68 (MQ=255)
              tGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCTCGACGAAGGCgag                                                <  1:380640/70‑1 (MQ=255)
                        cgATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCTCGACGAAGGCGAGAAcaca                                          >  1:2620773/1‑66 (MQ=255)
                        cgATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCTCGACGAAGGCGAGAAcaca                                          >  1:3893952/1‑66 (MQ=255)
                        cgATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCTCGACGAAGGCGAGAAcaca                                          >  1:658742/1‑66 (MQ=255)
                             ttCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGt                                >  1:3817437/1‑71 (MQ=255)
                               cgtcgCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGc                                    <  1:3108444/65‑1 (MQ=255)
                               cgtcgCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAg                                 <  1:3662855/68‑1 (MQ=255)
                                     ccccAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTg                         <  1:3821216/70‑1 (MQ=255)
                                     ccccAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTg                         <  1:411908/70‑1 (MQ=255)
                                     ccccAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTg                         <  1:1706323/70‑1 (MQ=255)
                                     ccccAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTg                         <  1:104023/70‑1 (MQ=255)
                                     ccccAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTg                         <  1:1263469/70‑1 (MQ=255)
                                          gTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAg                                                     >  1:2534936/1‑37 (MQ=255)
                                          gTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCaa                                  <  1:1457638/56‑1 (MQ=255)
                                          gTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCctct                   >  1:1919666/1‑71 (MQ=255)
                                          gTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCctct                   >  1:2604331/1‑71 (MQ=255)
                                          gTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGGGGTTTTGCCctct                   >  1:3690517/1‑71 (MQ=255)
                                           ttCCGTTGCGGTACGATAAAGGATACCTCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTg                  >  1:3101376/1‑71 (MQ=255)
                                           ttCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGc                                    <  1:204836/53‑1 (MQ=255)
                                              cGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGa                <  1:1514913/70‑1 (MQ=255)
                                                  gCGGTACGATAAAGGATACCTCGACGAAGGCGAGAACACAg                                         >  1:2734672/1‑41 (MQ=255)
                                                  gCGGTACGATAAAGGATACCTCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtct            <  1:1706524/70‑1 (MQ=255)
                                                  gCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtctt           >  1:1731217/1‑71 (MQ=255)
                                                  gCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtctt           >  1:1888908/1‑71 (MQ=255)
                                                  gCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtctt           >  1:3972704/1‑71 (MQ=255)
                                                    ggTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAAtcttct         >  1:3234633/1‑71 (MQ=255)
                                                      tACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGc                        >  1:2999739/1‑54 (MQ=255)
                                                      tACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGc                        >  1:2982866/1‑54 (MQ=255)
                                                      tACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGc                        >  1:2508137/1‑54 (MQ=255)
                                                      tACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGc                        >  1:1578277/1‑54 (MQ=255)
                                                       aCGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCtt      >  1:656976/1‑71 (MQ=255)
                                                       aCGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCtt      >  1:1917339/1‑71 (MQ=255)
                                                        cGATAAAGGATACCTCGACGAAGGCGAGAACACAGTGAGCAAGTGGttt                           >  1:307108/1‑49 (MQ=255)
                                                        cGATAAAGGATACCTCGACGAAGGCGAGAACACAGTGAGCAAGTGGttt                           >  1:2643647/1‑49 (MQ=255)
                                                        cGATAAAGGATACCGGGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCttt     >  1:290126/1‑71 (MQ=255)
                                                        cGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCttt     >  1:1416618/1‑71 (MQ=255)
                                                        cGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCttt     >  1:3556537/1‑71 (MQ=255)
                                                          aTAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCa   >  1:2229240/1‑71 (MQ=255)
                                                          aTAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCa   >  1:1929514/1‑71 (MQ=255)
                                                          aTAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCa   >  1:1555673/1‑71 (MQ=255)
                                                           tAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCAc  <  1:947189/71‑1 (MQ=255)
                                                           tAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCAc  <  1:3582913/71‑1 (MQ=255)
                                                           tAAAGGATACCGAGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCAc  <  1:3883197/71‑1 (MQ=255)
                                                            aaaGGATACCTCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTATCAc  >  1:2499392/1‑70 (MQ=255)
                                                            aaaGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCAc  >  1:3648369/1‑70 (MQ=255)
                                                                      |                                                           
TATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCAC  >  W3110S.gb/1413203‑1413332

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: