Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 108,684 G→A 9.1% A136T (GCG→ACG)  secA → preprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb108,6840GA 9.1% 171.4 / 6.2 55A136T (GCG→ACG) secApreprotein translocase subunit, ATPase that targets protein precursors to the SecYE core translocon
Reads supporting (aligned to +/- strand):  ref base G (20/30);  new base A (2/3);  total (22/33)
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.86e-01

ACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCC  >  W3110S.gb/108615‑108754
                                                                     |                                                                      
aCGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgc                                                                       >  1:3976339/1‑71 (MQ=255)
 cGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgcg                                                                      >  1:2740806/1‑71 (MQ=255)
  gCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgcgc                                                                     <  1:3631782/71‑1 (MQ=255)
  gCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgcgc                                                                     <  1:501953/71‑1 (MQ=255)
  gCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgcgc                                                                     <  1:4397629/71‑1 (MQ=255)
  gCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgcg                                                                      >  1:4598612/1‑70 (MQ=255)
   ctgcctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGgcgc                                                                     >  1:4380237/1‑70 (MQ=255)
   ctgcctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGACGCa                                                                    <  1:746333/71‑1 (MQ=255)
   ctgcctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGACGCa                                                                    <  1:628845/71‑1 (MQ=255)
   ctgcctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGACGCa                                                                    <  1:1481159/71‑1 (MQ=255)
      cctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACg                                                                 >  1:4302725/1‑71 (MQ=255)
       ctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGt                                                                >  1:3353293/1‑71 (MQ=255)
       ctgcTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGt                                                                <  1:1154035/71‑1 (MQ=255)
          cTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGAc                                                             >  1:3960825/1‑71 (MQ=255)
              ccTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGcc                                                          <  1:4746769/70‑1 (MQ=255)
                 gAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGaa                                                       >  1:4190489/1‑70 (MQ=255)
                        cTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCg                                               <  1:2829358/71‑1 (MQ=255)
                        cTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCg                                               >  1:3076552/1‑71 (MQ=255)
                          aaCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGACGCAACGTGACGCCGAAAACAACCGTc                                             >  1:3595296/1‑71 (MQ=255)
                           aCCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGACGCAACGTGACGCCGAAAACAACCGTc                                             >  1:1071098/1‑70 (MQ=255)
                               gTAAAGGCGTGCACGTAGTTACCGTCCACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCctct                                         >  1:2784699/1‑70 (MQ=255)
                                 aaaGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGt                                       >  1:1599535/1‑70 (MQ=255)
                                   aGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTg                                    <  1:4493451/71‑1 (MQ=255)
                                   aGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTg                                    <  1:4278638/71‑1 (MQ=255)
                                   aGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTg                                    <  1:2403843/71‑1 (MQ=255)
                                   aGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTg                                    <  1:1907617/71‑1 (MQ=255)
                                   aGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTg                                    <  1:4120265/71‑1 (MQ=255)
                                   aGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTg                                    <  1:4091940/71‑1 (MQ=255)
                                    ggCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGa                                   <  1:4205949/71‑1 (MQ=255)
                                          cACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTcc                              <  1:1960262/70‑1 (MQ=255)
                                          cACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCt                             >  1:3041272/1‑71 (MQ=255)
                                              tAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGc                         <  1:3835706/71‑1 (MQ=255)
                                                gTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCt                       <  1:4068430/71‑1 (MQ=255)
                                                gTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCt                       <  1:3616245/71‑1 (MQ=255)
                                                 ttACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCg                                               <  1:4531314/46‑1 (MQ=255)
                                                 ttACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCg                                               <  1:1475453/46‑1 (MQ=255)
                                                 ttACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCg                                               <  1:2358939/46‑1 (MQ=255)
                                                 ttACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGc                                          <  1:3311754/51‑1 (MQ=255)
                                                 ttACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAAt                                 <  1:766028/60‑1 (MQ=255)
                                                 ttACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAAt                                 <  1:5095571/60‑1 (MQ=255)
                                                 ttACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCt                       >  1:1036712/1‑70 (MQ=255)
                                                 ttACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCt                       >  1:1754682/1‑70 (MQ=255)
                                                  tACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGa                     >  1:2071655/1‑71 (MQ=255)
                                                  tACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGa                     >  1:738182/1‑71 (MQ=255)
                                                  aaCCGGCAACGACTACCTGGCGCCACGTGACGCCGAAAACAACCg                                               <  1:1864463/44‑1 (MQ=255)
                                                   aCCGTCAACGAGTACCTGGCGCAACGTGACGCCGAAAACAACCATCCGCTGTTTGaa                                  >  1:3591925/1‑57 (MQ=255)
                                                   aCCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGaa                                  >  1:2160982/1‑57 (MQ=255)
                                                   aCCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGaa                                  >  1:1431891/1‑57 (MQ=255)
                                                     cGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTg                                    <  1:1126285/53‑1 (MQ=255)
                                                     cGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTc                               <  1:3063738/58‑1 (MQ=255)
                                                       tCAACGACTACCTGGCGCAACGTGACGCCGAAaaca                                                   <  1:984590/36‑1 (MQ=255)
                                                       tCAACGACTACCTGGCGCAACGTGACGCCGAAaaca                                                   <  1:3239/36‑1 (MQ=255)
                                                       tCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTc                <  1:1828873/71‑1 (MQ=255)
                                                         aaCGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGttt                                     >  1:3213563/1‑48 (MQ=255)
                                                         aaCGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGaa                                  >  1:3926787/1‑51 (MQ=255)
                                                         aaCGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGaa                                  >  1:4066813/1‑51 (MQ=255)
                                                         aaCGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCg               >  1:334394/1‑70 (MQ=255)
                                                          aCGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTc                               >  1:3474379/1‑53 (MQ=255)
                                                                 ccTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCaa        <  1:3285073/69‑1 (MQ=255)
                                                                 ccTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCaa        <  1:12658/69‑1 (MQ=255)
                                                                  cTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTg                           <  1:334643/49‑1 (MQ=255)
                                                                   tGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAAt                                 <  1:311438/42‑1 (MQ=255)
                                                                    ggCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGaa                                  >  1:4585652/1‑40 (MQ=255)
                                                                    ggCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGaa                                  >  1:492683/1‑40 (MQ=255)
                                                                    ggCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGaa                                  >  1:1545074/1‑40 (MQ=255)
                                                                    ggCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGc   >  1:4655788/1‑71 (MQ=255)
                                                                     gcgcAACGTGACGCCGAAAACAACCGTCCGCTGTTTg                                    <  1:1041785/37‑1 (MQ=255)
                                                                     gcgcAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGcc  >  1:1828697/1‑71 (MQ=255)
                                                                     |                                                                      
ACGCTGCCTGCTTACCTGAACGCACTAACCGGTAAAGGCGTGCACGTAGTTACCGTCAACGACTACCTGGCGCAACGTGACGCCGAAAACAACCGTCCGCTGTTTGAATTCCTTGGCCTGACTGTCGGTATCAACCTGCC  >  W3110S.gb/108615‑108754

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: