Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,079 C→T 55.8% R231R (CGC→CGT araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,0790CT55.8% 12.1 / 70.4 52R231R (CGC→CGTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (6/17);  new base T (18/11);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGT  >  W3110S.gb/71009‑71148
                                                                      |                                                                      
cGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGc                                                                        <  1:928212/71‑1 (MQ=255)
         tGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAg                                                               >  1:862084/1‑71 (MQ=255)
          gTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAgg                                                              >  1:184423/1‑71 (MQ=255)
          gTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCa                                                                       >  1:2497128/1‑62 (MQ=255)
           tCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGc                                                             >  1:706471/1‑71 (MQ=255)
             acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGa                                                           <  1:881256/71‑1 (MQ=255)
             acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGa                                                           <  1:50492/71‑1 (MQ=255)
                  ttttCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAgc                                                        >  1:1994297/1‑69 (MQ=255)
                     tCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCa                                                                >  1:1978761/1‑58 (MQ=255)
                       cGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtttt                                                <  1:1107146/71‑1 (MQ=255)
                         ccAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGa                                              <  1:1060487/71‑1 (MQ=255)
                         ccAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGACGAAGCTGCTTTTGa                                              <  1:15846/71‑1 (MQ=255)
                         ccAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTg                                               <  1:1166833/70‑1 (MQ=255)
                          cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAg                                                               >  1:1434085/1‑54 (MQ=255)
                          cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTa                                                                    >  1:468178/1‑49 (MQ=255)
                          cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                             <  1:234164/71‑1 (MQ=255)
                           agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTgcgc                                            <  1:1132136/71‑1 (MQ=255)
                            gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAg                                                         >  1:2178080/1‑58 (MQ=255)
                            gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGAGAAg                                                         >  1:1622417/1‑58 (MQ=255)
                             cagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCa                                           >  1:100983/1‑70 (MQ=255)
                                ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtt                                                  <  1:1704009/60‑1 (MQ=255)
                                ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtt                                                  <  1:269037/60‑1 (MQ=255)
                                ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                             <  1:1798790/65‑1 (MQ=255)
                                 tAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAcc                                      >  1:38106/1‑71 (MQ=255)
                                  aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCaccacc                                      >  1:902312/1‑70 (MQ=255)
                                  aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCAccc                                     >  1:533024/1‑71 (MQ=255)
                                    ggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgctgct                                                   <  1:1656406/55‑1 (MQ=255)
                                     gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGa                                  >  1:1878119/1‑71 (MQ=255)
                                      aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                             <  1:1965947/59‑1 (MQ=255)
                                       ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAg                                                               <  1:636536/41‑1 (MQ=255)
                                         aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGcc                              >  1:909879/1‑71 (MQ=255)
                                         aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGc                               >  1:1670644/1‑70 (MQ=255)
                                         aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGcc                              >  1:1866459/1‑71 (MQ=255)
                                           cGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCttt                                                 <  1:1210404/50‑1 (MQ=255)
                                            gTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                    <  1:494681/62‑1 (MQ=255)
                                             tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgct                                                       <  1:2445914/43‑1 (MQ=255)
                                             tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAATCTGCtttt                                                <  1:69328/49‑1 (MQ=255)
                                               ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTg                                               <  1:1020360/48‑1 (MQ=255)
                                               ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                    <  1:1021688/59‑1 (MQ=255)
                                               ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGc                        >  1:2241974/1‑71 (MQ=255)
                                               ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGc                        >  1:1773592/1‑71 (MQ=255)
                                                tAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTg                                               <  1:2534438/47‑1 (MQ=255)
                                                 aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGcc                       <  1:915396/70‑1 (MQ=255)
                                                 aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCa                      >  1:228817/1‑71 (MQ=255)
                                                 aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCa                      >  1:968964/1‑71 (MQ=255)
                                                 aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTTGAGCACTACCCGGa                                  >  1:1584011/1‑60 (MQ=255)
                                                   gCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCAcc                    <  1:325085/71‑1 (MQ=255)
                                                   gCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCAcc                    <  1:1763541/71‑1 (MQ=255)
                                                    cTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCg                   <  1:987292/71‑1 (MQ=255)
                                                       gcgcGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCgtc                 <  1:2114991/70‑1 (MQ=255)
                                                        cgcgAGGACCAACGTATCAGCCAGGCGAAgctgct                                                   <  1:664068/35‑1 (MQ=25)
                                                              gACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGa                                  >  1:837221/1‑46 (MQ=255)
                                                                     gCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGt  <  1:2119662/71‑1 (MQ=255)
                                                                      |                                                                      
CGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGT  >  W3110S.gb/71009‑71148

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: