Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,109 T→C 60.8% T241T (ACT→ACC araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,1090TC60.8% 40.4 / 56.5 51T241T (ACT→ACCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (4/16);  new base C (22/9);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 5.31e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGC  >  W3110S.gb/71040‑71173
                                                                     |                                                                
gTTAGGGATTAGCGTCTTAAGCTGGCGCGATGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTa                                                                 <  1:1495630/71‑1 (MQ=255)
gTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTa                                                                 <  1:1615796/71‑1 (MQ=255)
gTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTa                                                                 <  1:155394/71‑1 (MQ=255)
    gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                             <  1:1471141/71‑1 (MQ=255)
    gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                             <  1:1382576/71‑1 (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:1606542/1‑71 (MQ=255)
             gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                                                    >  1:816946/1‑71 (MQ=255)
             gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                                                    >  1:1095677/1‑71 (MQ=255)
              tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:1251758/61‑1 (MQ=255)
              tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:401953/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:1580394/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: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)
                                           agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa                      >  1:154360/4‑71 (MQ=255)
                                           agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa                      >  1:278010/4‑71 (MQ=255)
                                           aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTg                       <  1:65356/70‑1 (MQ=255)
                                             ccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACg                    <  1:967638/70‑1 (MQ=255)
                                             ccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACg                    <  1:290148/70‑1 (MQ=255)
                                                 gCGAAGCTGCTTTTGAGCACCACCCGGATGCCTAt                                                    >  1:954661/1‑35 (MQ=255)
                                                      gctgctTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAActc            >  1:1223096/1‑70 (MQ=255)
                                                           ttttGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAt         <  1:345021/68‑1 (MQ=255)
                                                           ttttGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAt         <  1:1128498/68‑1 (MQ=255)
                                                           ttttGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAt         <  1:1289768/68‑1 (MQ=255)
                                                           ttttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGGTCAACTCTAtt        <  1:1255243/69‑1 (MQ=255)
                                                           ttttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAtt        <  1:1000551/69‑1 (MQ=255)
                                                           ttttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAtt        <  1:1067465/69‑1 (MQ=255)
                                                             ttGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtctc    >  1:222444/1‑71 (MQ=255)
                                                             ttGAGCACCACCCGGATGCCTATCGCCACCgtcggtc                                      >  1:769225/1‑37 (MQ=255)
                                                             ttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGttt                         >  1:526502/1‑50 (MQ=255)
                                                             ttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtc      >  1:1369061/1‑69 (MQ=255)
                                                             ttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtc      >  1:813131/1‑69 (MQ=255)
                                                             ttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtc      >  1:419078/1‑69 (MQ=255)
                                                             ttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtc      >  1:1770263/1‑69 (MQ=255)
                                                                 gCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTcgc  >  1:1673627/1‑69 (MQ=255)
                                                                  caccacCCGGATGCCTATCGCCACCGTCGGTCTCAATGTTGGTTTTGACGATCAACTCTAt         <  1:446539/61‑1 (MQ=255)
                                                                  caccacCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAt         <  1:57943/61‑1 (MQ=255)
                                                                     |                                                                
GTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGC  >  W3110S.gb/71040‑71173

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