Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,085 T→C 57.7% 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,0850TC57.7% 28.5 / 63.5 52S233S (AGT→AGCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (2/20);  new base C (17/13);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 4.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTG  >  W3110S.gb/71017‑71152
                                                                    |                                                                   
cTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGt                                                                     >  1:1723731/1‑69 (MQ=255)
cTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGt                                                                     >  1:1147586/1‑69 (MQ=255)
      caTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCg                                                               <  1:59097/69‑1 (MQ=255)
        tCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAgc                                                           >  1:248033/1‑71 (MQ=255)
         cTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAgct                                                          >  1:639260/1‑71 (MQ=255)
                 ccAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCtttt                                                    >  1:1575570/1‑69 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                                 <  1:1623118/71‑1 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                                 <  1:1552906/71‑1 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                 <  1:423051/71‑1 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                 <  1:330608/71‑1 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                 <  1:1675946/71‑1 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                 <  1:1505071/71‑1 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                 <  1:1603654/71‑1 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                 <  1:1572234/71‑1 (MQ=255)
                  cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                 <  1:1640049/71‑1 (MQ=255)
                   agcagTTAGGGATTAGCGTCTTAAGCTTGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                                <  1:832324/71‑1 (MQ=255)
                       gTTAGGGATTAGCGTCTTAAGCTGGCGCGATGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTa                                            <  1:1495630/71‑1 (MQ=255)
                       gTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTa                                            <  1:155394/71‑1 (MQ=255)
                       gTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTa                                            <  1:1615796/71‑1 (MQ=255)
                        ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAcca                                            <  1:952942/70‑3 (MQ=255)
                          aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAgc                                                           >  1:1652774/1‑53 (MQ=255)
                           gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                        <  1:1471141/71‑1 (MQ=255)
                           gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                        <  1:1382576/71‑1 (MQ=255)
                               ttAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCtttt                                                    <  1:721403/55‑1 (MQ=255)
                               ttAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCtttt                                                    <  1:1761698/55‑1 (MQ=255)
                               ttAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCtttt                                                    <  1:292580/55‑1 (MQ=255)
                               ttAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCtttt                                                    <  1:409654/55‑1 (MQ=255)
                                tAGCGTCTTAAGCTGGCGCGAGGACCAACGTATAAGCCAGGCGaa                                                             >  1:424336/1‑45 (MQ=255)
                                  gCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGcc                                  <  1:612425/70‑1 (MQ=255)
                                    gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTa                                >  1:816799/1‑70 (MQ=255)
                                    gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                               >  1:1095677/1‑71 (MQ=255)
                                    gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                               >  1:1820069/1‑71 (MQ=255)
                                    gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                               >  1:816946/1‑71 (MQ=255)
                                    gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                               >  1:1606542/1‑71 (MQ=255)
                                     tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                        <  1:401953/61‑1 (MQ=255)
                                     tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                        <  1:1251758/61‑1 (MQ=255)
                                      cTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                        <  1:171921/60‑1 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGa                                      >  1:1280168/1‑61 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGa                                      >  1:519084/1‑61 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGc                            >  1:22384/1‑71 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGc                            >  1:952835/1‑71 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGc                            >  1:1580394/1‑71 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGa                                      >  1:931581/1‑61 (MQ=255)
                                       ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGc                            >  1:689862/1‑71 (MQ=255)
                                         aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGcc                           <  1:434297/70‑1 (MQ=255)
                                         aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGc                            <  1:1478267/69‑1 (MQ=255)
                                                    aGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcggtcg                >  1:1797031/1‑70 (MQ=255)
                                                    aGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcggtcg                >  1:1244049/1‑70 (MQ=255)
                                                    aGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGc               <  1:1111290/71‑1 (MQ=255)
                                                             gCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGt      <  1:1155635/71‑1 (MQ=255)
                                                                ttAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCt                                 <  1:787514/41‑1 (MQ=255)
                                                                ttAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCt                                 <  1:580991/41‑1 (MQ=255)
                                                                ttAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCt                                 <  1:49073/41‑1 (MQ=255)
                                                                  aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTg  <  1:65356/70‑1 (MQ=255)
                                                                    |                                                                   
CTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTG  >  W3110S.gb/71017‑71152

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: