Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,448,059 T→C 10.0% P1091P (CCA→CCG rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,448,0590TC10.0% 101.1 / 3.3 40P1091P (CCA→CCGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (14/22);  new base C (2/2);  total (16/24)
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.19e-01

CTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAGATCTTTACCA  >  W3110S.gb/3447990‑3448122
                                                                     |                                                               
cTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTg                                                                >  1:276477/1‑71 (MQ=255)
        cATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGa                                                        <  1:56073/71‑1 (MQ=255)
            ttCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAAcgt                                                    >  1:1983632/1‑71 (MQ=255)
                aGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGtt                                                <  1:1324581/71‑1 (MQ=255)
                     aaCAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTAcc                                             <  1:2351019/69‑1 (MQ=255)
                            gCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTg                                                                <  1:1096963/43‑1 (MQ=255)
                            gCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTg                                                                <  1:1238989/43‑1 (MQ=255)
                            gCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTg                                                                <  1:810949/43‑1 (MQ=255)
                               ttACCCGGCAGGAAGTACTGCGCAGGCATATCGGTACCCGGGATCAGAACGTCGTTACCCTGAGCATCAAc                                 >  1:378391/1‑71 (MQ=255)
                               ttACCCGGCAGGAAGTACTGCGCAGGCATATCGGTACCCGGGATCAGAACGTCGTTACCCTGAGCATCAAc                                 >  1:3346083/1‑71 (MQ=255)
                                  cccGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGa                               <  1:1524477/70‑1 (MQ=255)
                                     ggCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGAtttt                           <  1:3084531/71‑1 (MQ=255)
                                     ggCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGAtttt                           <  1:2233277/71‑1 (MQ=255)
                                      gCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTc                          <  1:539980/71‑1 (MQ=255)
                                       catgAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTc                          <  1:1989379/67‑1 (MQ=255)
                                        aGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAg                        <  1:194173/71‑1 (MQ=255)
                                          gAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTAc                                              >  1:3600450/1‑47 (MQ=255)
                                          gAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTg                      >  1:1026333/1‑71 (MQ=255)
                                           aaGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGc                     >  1:2751974/1‑71 (MQ=255)
                                           aaGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGc                     >  1:1918518/1‑71 (MQ=255)
                                            aGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGcc                    >  1:2315767/1‑71 (MQ=255)
                                                cTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGa                               <  1:1538269/56‑1 (MQ=255)
                                                cTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGa                               <  1:2929468/56‑1 (MQ=255)
                                                   cgcTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAg                        >  1:1719748/1‑60 (MQ=255)
                                                   cgcTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACg               <  1:2958797/69‑1 (MQ=255)
                                                    gcaggCATATCGGTACCCGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAg            <  1:2014990/68‑1 (MQ=255)
                                                    gcaggCATATCGGTACCCGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAg            <  1:3267181/68‑1 (MQ=255)
                                                        gCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCaa                                  >  1:3065006/1‑45 (MQ=255)
                                                        gCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAGATc         <  1:972907/70‑1 (MQ=255)
                                                         cATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAGATCtt       <  1:1870925/71‑1 (MQ=255)
                                                           tatCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAg            <  1:664155/64‑1 (MQ=255)
                                                            atCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGAtttt                           <  1:1743385/48‑1 (MQ=255)
                                                            atCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGAtttt                           <  1:1516023/48‑1 (MQ=255)
                                                             tCGGTACCTGTGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAg            <  1:1189649/62‑1 (MQ=255)
                                                             tCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAAc                                 >  1:1031921/1‑41 (MQ=255)
                                                             tCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAAc                                 >  1:1980670/1‑41 (MQ=255)
                                                             tCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAGATCttt      <  1:3797061/68‑1 (MQ=255)
                                                              cGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAGATCTTTacca  >  1:1152929/1‑71 (MQ=255)
                                                               ggTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAGATCTTTac    >  1:2076391/1‑68 (MQ=255)
                                                               ggTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAGATCTTTac    >  1:1383781/1‑68 (MQ=255)
                                                                     |                                                               
CTGTACGCCATCTTCCAGCTGAACAATCGCTTTACCCGGCAGGAAGTACTGCGCTGGCATATCGGTACCTGGGATCAGAACGTCGTTACCCTGAGCATCAACGATTTTCAGTGCCGGACGCAGATCTTTACCA  >  W3110S.gb/3447990‑3448122

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: