Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,085 T→C 37.5% S233S (AGT→AGC araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,0850TC37.5% 36.7 / 51.9 48S233S (AGT→AGCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (13/17);  new base C (6/12);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGAC  >  W3110S.gb/71016‑71154
                                                                     |                                                                     
tCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTc                                                                      <  1:1292103/71‑1 (MQ=255)
      acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGATGACCAACGTATCAGCCAGGCGa                                                                <  1:1567134/71‑1 (MQ=255)
       caTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCg                                                                 <  1:1702549/69‑1 (MQ=255)
         tCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGt                                                                       <  1:1583579/61‑1 (MQ=255)
             ttCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAgctgc                                                          <  1:193084/70‑1 (MQ=255)
              tCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgctgct                                                         <  1:1088847/70‑1 (MQ=255)
              tCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGAATTAGTCAGGCGAAgctgct                                                         <  1:1325499/70‑1 (MQ=255)
                cGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtttt                                                      <  1:1220072/71‑1 (MQ=255)
                   cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                   >  1:936518/1‑71 (MQ=255)
                   cagcagTTAGGGATTAGCGTCTTAAGCTGGAGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                                   <  1:1198936/71‑1 (MQ=255)
                     gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                               >  1:1153614/1‑57 (MQ=255)
                         ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCa                                                                     <  1:190994/47‑1 (MQ=255)
                         ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCg                                                                 <  1:362460/51‑1 (MQ=255)
                         ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGt                                                                       <  1:1553282/45‑1 (MQ=255)
                         ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGa                                                                <  1:1106285/52‑1 (MQ=255)
                         ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTg                                                     <  1:1772567/63‑1 (MQ=255)
                         ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                   <  1:1587896/65‑1 (MQ=255)
                           aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCttt                                                       >  1:378156/1‑59 (MQ=255)
                           aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAccc                                           >  1:723223/1‑71 (MQ=255)
                           aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAcc                                            >  1:1811410/1‑70 (MQ=255)
                            gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                          <  1:1586818/71‑1 (MQ=255)
                            gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                          <  1:504736/71‑1 (MQ=255)
                            gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                          <  1:480224/71‑1 (MQ=255)
                            gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                          <  1:1212758/71‑1 (MQ=255)
                                ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAg                                                              <  1:2183866/47‑1 (MQ=255)
                                 tAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCg                                                                 <  1:916111/43‑1 (MQ=37)
                                 tAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                   >  1:1595009/1‑57 (MQ=255)
                                   gCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                          <  1:2195743/64‑1 (MQ=255)
                                      tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                          <  1:588106/61‑1 (MQ=255)
                                        ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGa                                        >  1:1373399/1‑61 (MQ=255)
                                        ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGa                                                    >  1:1844775/1‑49 (MQ=255)
                                        ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTa                                              >  1:752692/1‑55 (MQ=255)
                                        ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                          <  1:536925/59‑1 (MQ=255)
                                          aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGa                                        >  1:195739/1‑59 (MQ=255)
                                          aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGc                              <  1:1084591/69‑1 (MQ=255)
                                          aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                            >  1:504204/1‑71 (MQ=255)
                                                gcgcGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtc                       <  1:1362532/70‑1 (MQ=255)
                                                gcgcGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgt                        <  1:1075670/69‑1 (MQ=255)
                                                 cgcgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCgtcgg                     <  1:119649/71‑1 (MQ=255)
                                                   cgAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcggtc                   >  1:934848/1‑71 (MQ=255)
                                                    gAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcggtc                   >  1:1689097/1‑70 (MQ=255)
                                                      ggcccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGAACTACCCGGATGCCTATCGCCACCGTCGGTCGc                 >  1:2128127/4‑70 (MQ=255)
                                                       gACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCaa               >  1:926012/1‑71 (MQ=255)
                                                         ccAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTa                                  <  1:1206202/50‑1 (MQ=255)
                                                            aCGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                                 >  1:1860587/1‑48 (MQ=255)
                                                            aCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTg          >  1:706743/1‑71 (MQ=255)
                                                            aCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTg          >  1:1104435/1‑71 (MQ=255)
                                                              gCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGt        >  1:1249649/1‑71 (MQ=255)
                                                                   aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGt            <  1:462567/62‑1 (MQ=255)
                                                                    gTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGAc  <  1:1719171/71‑1 (MQ=255)
                                                                     tCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTg          >  1:1622363/1‑62 (MQ=255)
                                                                     |                                                                     
TCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGAC  >  W3110S.gb/71016‑71154

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: