Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,459,280 T→C 10.2% E15G (GAA→GGA)  secE ← preprotein translocase membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,459,2800TC10.2% 154.9 / 9.9 59E15G (GAA→GGA) secEpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base C (2/4);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TCGCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTCTAC  >  W3110S.gb/3459211‑3459349
                                                                     |                                                                     
tCGCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTc                                                                      >  1:1946/1‑71 (MQ=255)
     aTAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAgg                                                                  <  1:1717606/70‑1 (MQ=255)
      tAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGcc                                                                >  1:555207/1‑71 (MQ=255)
       aaaGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGcc                                                                >  1:2060834/1‑70 (MQ=255)
       aaaGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCcg                                                               >  1:2548613/1‑71 (MQ=255)
        aaGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTc                                                                      <  1:2796282/63‑1 (MQ=255)
          gATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCa                                                                    <  1:1802880/63‑1 (MQ=255)
            tAGTTGCCGACAATGGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGcc                                                           >  1:1122826/1‑70 (MQ=255)
                 gCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCtt                                                      <  1:2668787/70‑1 (MQ=255)
                   cGACAATCGCCACCAGGAGCAATGCCACAACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCt                                                   >  1:661323/1‑71 (MQ=255)
                    gACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCtt                                                  <  1:1679342/71‑1 (MQ=255)
                     aCAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCtt                                                  <  1:4198473/70‑1 (MQ=255)
                     aCAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCtt                                                  >  1:2082801/1‑70 (MQ=255)
                      cAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCtt                                                  >  1:3453249/1‑69 (MQ=255)
                      cAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCtt                                                  >  1:248540/1‑69 (MQ=255)
                      cAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGa                                                <  1:384720/71‑1 (MQ=255)
                         tCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGc                                                        <  1:3205012/60‑1 (MQ=255)
                         tCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCt                                             >  1:62221/1‑71 (MQ=255)
                          cGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGc                                                        <  1:881339/59‑1 (MQ=255)
                            ccaccaGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCg                                          >  1:3393180/1‑71 (MQ=255)
                             caccaGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCtt                                            >  1:4075479/1‑68 (MQ=255)
                             caccaGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCgg                                         <  1:2209972/71‑1 (MQ=255)
                             caccaGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCgg                                         <  1:3719060/71‑1 (MQ=255)
                               ccaGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTcc                                                    <  1:530360/58‑1 (MQ=255)
                                 aGGAGCAATGCCACCACAACGACCGACTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAt                                      >  1:3317972/1‑70 (MQ=255)
                                 aGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCa                                                                    <  1:3953255/40‑1 (MQ=255)
                                  ggAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAg                                               >  1:693501/1‑60 (MQ=255)
                                    aGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCg                                   >  1:2875511/1‑70 (MQ=255)
                                     gCAATGCCACCACAACGACCCACTTCATCGCTTCCAg                                                                   <  1:753265/37‑1 (MQ=255)
                                     gCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTc                                    >  1:1978194/1‑68 (MQ=255)
                                      cAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCa                                 >  1:1950929/1‑70 (MQ=255)
                                         tGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCa                             <  1:3601694/71‑1 (MQ=255)
                                           ccaccaCAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCTCa                                 >  1:4254777/1‑65 (MQ=255)
                                           ccaccaCAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATa                           <  1:715289/71‑1 (MQ=255)
                                           ccaccaCAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATa                           <  1:3080148/71‑1 (MQ=255)
                                           ccaccaCAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATa                           <  1:1802275/71‑1 (MQ=255)
                                             accacAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATaaa                         >  1:3990353/1‑71 (MQ=255)
                                                acaACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTc                                           <  1:4179393/50‑1 (MQ=255)
                                                acaACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaacc                       >  1:2975555/1‑70 (MQ=255)
                                                 caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAg                                               >  1:3073899/1‑45 (MQ=255)
                                                 caACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaacaac                     >  1:133142/1‑68 (MQ=255)
                                                   aCGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAaaccaacc                   >  1:266920/1‑71 (MQ=255)
                                                   aCGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTACTTGAGCTTCGGt                                        >  1:4244893/1‑50 (MQ=255)
                                                     gACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTg                 >  1:1880577/1‑71 (MQ=255)
                                                          aCTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTAtt                                     >  1:1809816/1‑46 (MQ=255)
                                                          aCTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGa            <  1:2443669/71‑1 (MQ=255)
                                                           cTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTa                                       <  1:2937870/43‑1 (MQ=255)
                                                           cTTCATCGCTCCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCAt                            <  1:2982277/54‑1 (MQ=255)
                                                             tCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGt         >  1:1838715/1‑71 (MQ=255)
                                                              cATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTc               >  1:104358/1‑64 (MQ=255)
                                                              cATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTc               >  1:81955/1‑64 (MQ=255)
                                                               aTCGCTTCCAGGCCGCTCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGtt        <  1:4028705/70‑1 (MQ=255)
                                                               aTCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGtt        <  1:4039423/70‑1 (MQ=255)
                                                               aTCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGtt        <  1:871599/70‑1 (MQ=255)
                                                                tCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGtttt      <  1:815316/71‑1 (MQ=255)
                                                                tCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGttt       >  1:2910282/1‑70 (MQ=255)
                                                                 cGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTc     >  1:1025652/1‑71 (MQ=255)
                                                                  gCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTc     <  1:224248/70‑1 (MQ=255)
                                                                  gCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTc     >  1:3528378/1‑70 (MQ=255)
                                                                     tCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTCTAc  >  1:407608/1‑70 (MQ=255)
                                                                     |                                                                     
TCGCGATAAAGATAGTTGCCGACAATCGCCACCAGGAGCAATGCCACCACAACGACCCACTTCATCGCTTCCAGGCCGCGCCCGCTTCCTTGAGCTTCGGTATTCGCACTCATAAACCAACCTGTCAGAAGTTTTCTAC  >  W3110S.gb/3459211‑3459349

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: