Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,234,526 T→C 19.1% T28A (ACG→GCG)  dsbB ← oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,234,5260TC19.1% 100.2 / 20.4 47T28A (ACG→GCG) dsbBoxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base C (6/3);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

ACGCGCAGCGTTCATAAATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACAT  >  W3110S.gb/1234457‑1234587
                                                                     |                                                             
aCGCGCAGCGTTCATAAATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTc                                                              <  1:3245820/71‑1 (MQ=255)
   cgcAGCGTTCATAAATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTc                                                              <  1:270278/68‑1 (MQ=255)
   cgcAGCGTTCATAAATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTc                                                              <  1:264196/68‑1 (MQ=255)
     cAGCGTTCATAAATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTc                                                         <  1:3593372/71‑1 (MQ=255)
         gTTCATAAATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGt                                                     >  1:858988/1‑71 (MQ=255)
             aTAAATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCa                                                             >  1:1380600/1‑59 (MQ=255)
               aaaTACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGt                                                     <  1:929932/65‑1 (MQ=255)
                    cAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGt                                          <  1:3229437/71‑1 (MQ=255)
                    cAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGt                                          <  1:1566813/71‑1 (MQ=255)
                     agagCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTa                                         <  1:3088359/71‑1 (MQ=255)
                         cACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACg                                     >  1:1091464/1‑71 (MQ=255)
                           cGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGc                                    >  1:4677647/1‑70 (MQ=255)
                              aaGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATc                                <  1:4738394/71‑1 (MQ=255)
                              aaGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATc                                <  1:4737814/71‑1 (MQ=255)
                               aGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATc                                >  1:4818023/1‑70 (MQ=255)
                                ggTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCaa                              >  1:1617288/1‑71 (MQ=255)
                                ggTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCaa                              >  1:3827136/1‑71 (MQ=255)
                                 gTTTCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTaa                                        >  1:1430598/1‑60 (MQ=255)
                                  tttCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACa                            <  1:3996088/71‑1 (MQ=255)
                                   ttCAGTAACATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTa                                         <  1:100281/57‑1 (MQ=255)
                                         aaCATCACATGCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTaaa                                       >  1:30600/1‑53 (MQ=255)
                                             tCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGccc                 >  1:3186209/1‑71 (MQ=255)
                                             tCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGccc                 >  1:4561736/1‑71 (MQ=255)
                                             tCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATAACAGCCACGCGcccc                >  1:4648613/1‑71 (MQ=255)
                                               acaTGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCa                        >  1:1169389/1‑62 (MQ=255)
                                               acaTGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCg               >  1:4306991/1‑71 (MQ=255)
                                                  tGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTa                                         <  1:752101/42‑1 (MQ=255)
                                                  tGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc            <  1:3198423/71‑1 (MQ=255)
                                                  tGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGcc            <  1:4896957/71‑1 (MQ=255)
                                                   gCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCt           >  1:5033711/1‑71 (MQ=255)
                                                   gCTGGAACCACAGCGCCGCCTGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCt           >  1:2758001/1‑71 (MQ=255)
                                                   gCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCt           >  1:1187675/1‑71 (MQ=255)
                                                   gCTGGAACCACAGCGCCGCCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCt           >  1:5005557/1‑71 (MQ=255)
                                                    cTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTaaa                                       >  1:1087280/1‑42 (MQ=255)
                                                    cTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTaaa                                       >  1:4605975/1‑42 (MQ=255)
                                                     tGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACTCCATCAACAGCCACGCGCCCCGGCCt           <  1:2508341/69‑1 (MQ=255)
                                                     tGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCt           <  1:3736456/69‑1 (MQ=255)
                                                     tGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCt           <  1:3636711/69‑1 (MQ=255)
                                                      ggAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATc                                <  1:2946194/47‑1 (MQ=255)
                                                       gAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTtgtg       <  1:686516/71‑1 (MQ=255)
                                                       gAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTtg         >  1:2343013/1‑69 (MQ=255)
                                                       gAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTtg         >  1:2020500/1‑69 (MQ=255)
                                                        aaCCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCaa                              >  1:606930/1‑47 (MQ=255)
                                                            acaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGaa     >  1:753903/1‑68 (MQ=255)
                                                            acaGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACAt  >  1:2042176/1‑71 (MQ=255)
                                                             caGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGcc                  >  1:122744/1‑54 (MQ=255)
                                                             caGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACAt  >  1:3536106/1‑70 (MQ=255)
                                                                     tCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCAcg                      <  1:4354522/42‑1 (MQ=255)
                                                                     tCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCAcg                      <  1:1912415/42‑1 (MQ=255)
                                                                     tCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCAcg                      <  1:1889001/42‑1 (MQ=255)
                                                                     tCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCAcg                      <  1:1878546/42‑1 (MQ=255)
                                                                     |                                                             
ACGCGCAGCGTTCATAAATACAGAGCACGCAAGGTTTCAGTAACATCACATGCTGGAACCACAGCGCCGTCAGTTCCAGTGCCAGAGCAGTAAACGCCATCAACAGCCACGCGCCCCGGCCTTGTGAACAT  >  W3110S.gb/1234457‑1234587

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: