Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,085 T→C 56.9% 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,0850TC56.9% 25.8 / 62.6 51S233S (AGT→AGCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (5/17);  new base C (18/11);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 9.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGA  >  W3110S.gb/71018‑71153
                                                                   |                                                                     
tGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAg                                                                    >  1:862084/1‑71 (MQ=255)
 gTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAgg                                                                   >  1:184423/1‑71 (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)
                 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:269037/60‑1 (MQ=255)
                       ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtt                                                       <  1:1704009/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)
                                                                 aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa  >  1:1409125/1‑71 (MQ=255)
                                                                   |                                                                     
TGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGA  >  W3110S.gb/71018‑71153

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: