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→AGC) araCDNA‑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

TGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGA  >  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)
              cGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑tttt                                                     <  1:1107146/71‑1 (MQ=255)
                ccAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGa                                                   <  1:1060487/71‑1 (MQ=255)
                ccAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGACGAAGCTGC‑TTTTGa                                                   <  1:15846/71‑1 (MQ=255)
                ccAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTg                                                    <  1:1166833/70‑1 (MQ=255)
                 cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAg                                                                    >  1:1434085/1‑54 (MQ=255)
                 cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAg                                                  <  1:234164/71‑1 (MQ=255)
                  agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTgcgc                                                 <  1:1132136/71‑1 (MQ=255)
                   gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAg                                                              >  1:2178080/1‑58 (MQ=255)
                   gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGAGAAg                                                              >  1:1622417/1‑58 (MQ=255)
                    cagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCa                                                >  1:100983/1‑70 (MQ=255)
                       ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑tt                                                       <  1:269037/60‑1 (MQ=255)
                       ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑tt                                                       <  1:1704009/60‑1 (MQ=255)
                       ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAg                                                  <  1:1798790/65‑1 (MQ=255)
                        tAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTAcc                                           >  1:38106/1‑71 (MQ=255)
                         aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCaccacc                                           >  1:902312/1‑70 (MQ=255)
                         aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCAccc                                          >  1:533024/1‑71 (MQ=255)
                           ggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgctgc‑t                                                        <  1:1656406/55‑1 (MQ=255)
                            gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCGGa                                       >  1:1878119/1‑71 (MQ=255)
                             aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAg                                                  <  1:1965947/59‑1 (MQ=255)
                              ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAg                                                                    <  1:636536/41‑1 (MQ=255)
                                aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCGGATGcc                                   >  1:909879/1‑71 (MQ=255)
                                aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCGGATGc                                    >  1:1670644/1‑70 (MQ=255)
                                aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCGGATGcc                                   >  1:1866459/1‑71 (MQ=255)
                                  cGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑ttt                                                      <  1:1210404/50‑1 (MQ=255)
                                   gTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCg                                         <  1:494681/62‑1 (MQ=255)
                                    tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgct                                                            <  1:2445914/43‑1 (MQ=255)
                                    tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAATCTGC‑tttt                                                     <  1:69328/49‑1 (MQ=255)
                                      ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTg                                                    <  1:1020360/48‑1 (MQ=255)
                                      ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCg                                         <  1:1021688/59‑1 (MQ=255)
                                      ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCGGATGCCTATCGc                             >  1:2241974/1‑71 (MQ=255)
                                      ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCGGATGCCTATCGc                             >  1:1773592/1‑71 (MQ=255)
                                       tAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTg                                                    <  1:2534438/47‑1 (MQ=255)
                                        aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCGGATGCCTATCGcc                            <  1:915396/70‑1 (MQ=255)
                                        aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCGGATGCCTATCGCCa                           >  1:228817/1‑71 (MQ=255)
                                        aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCGGATGCCTATCGCCa                           >  1:968964/1‑71 (MQ=255)
                                        aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTTGAGCACTACCCGGa                                       >  1:1584011/1‑60 (MQ=255)
                                          gCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCGGATGCCTATCGCCAcc                         <  1:325085/71‑1 (MQ=255)
                                          gCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCGGATGCCTATCGCCAcc                         <  1:1763541/71‑1 (MQ=255)
                                           cTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCGGATGCCTATCGCCACCg                        <  1:987292/71‑1 (MQ=255)
                                              gcgcGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCGGATGCCTATCGCCACCgtc                      <  1:2114991/70‑1 (MQ=255)
                                               cgcgAGGACCAACGTATCAGCCAGGCGAAgctgc‑t                                                        <  1:664068/35‑1 (MQ=25)
                                                     gACCAACGTATCAGCCAGGCGAAGCTGC‑TTTTGAGCACCACCCGGa                                       >  1:837221/1‑46 (MQ=255)
                                                            gCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGt       <  1:2119662/71‑1 (MQ=255)
                                                                 aGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa  >  1:1409125/1‑71 (MQ=255)
                                                                   |                                                                     
TGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGC‑TTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGA  >  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: ‑