Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,079 C→T 62.3% 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,0790CT62.3% 41.1 / 58.5 53R231R (CGC→CGTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (4/16);  new base T (17/16);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 4.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGT  >  W3110S.gb/71015‑71148
                                                                |                                                                     
gTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAgc                                                                 <  1:622092/71‑2 (MQ=255)
gTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTAGCGCGAGGACCAACGTATCAgc                                                                 <  1:255237/71‑2 (MQ=255)
gTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCCTCTTAAGCTGGCGCGAGGACCAACGTATCAgc                                                                 <  1:221176/71‑2 (MQ=255)
  cTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGt                                                                 >  1:1147586/1‑69 (MQ=255)
  cTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGt                                                                 >  1:1723731/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:1552906/71‑1 (MQ=255)
                    cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                             <  1:1623118/71‑1 (MQ=255)
                    cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                             <  1:1505071/71‑1 (MQ=255)
                    cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                             <  1:1572234/71‑1 (MQ=255)
                    cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                             <  1:330608/71‑1 (MQ=255)
                    cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                             <  1:423051/71‑1 (MQ=255)
                    cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                             <  1:1603654/71‑1 (MQ=255)
                    cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                             <  1:1675946/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:1615796/71‑1 (MQ=255)
                         gTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTa                                        <  1:155394/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:409654/55‑1 (MQ=255)
                                 ttAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCtttt                                                <  1:1761698/55‑1 (MQ=255)
                                 ttAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCtttt                                                <  1:292580/55‑1 (MQ=255)
                                 ttAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCtttt                                                <  1:721403/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:1820069/1‑71 (MQ=255)
                                      gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                           >  1:816946/1‑71 (MQ=255)
                                      gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                           >  1:1606542/1‑71 (MQ=255)
                                      gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                           >  1:1095677/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:519084/1‑61 (MQ=255)
                                         ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGa                                  >  1:1280168/1‑61 (MQ=255)
                                         ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGc                        >  1:1580394/1‑71 (MQ=255)
                                         ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGc                        >  1:22384/1‑71 (MQ=255)
                                         ttAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGc                        >  1:952835/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:1244049/1‑70 (MQ=255)
                                                      aGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcggtcg            >  1:1797031/1‑70 (MQ=255)
                                                      aGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGc           <  1:1111290/71‑1 (MQ=255)
                                                               gCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGt  <  1:1155635/71‑1 (MQ=255)
                                                                |                                                                     
GTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGT  >  W3110S.gb/71015‑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: