Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,109 T→C 47.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,1090TC47.8% 2.1 / 101.8 69T241T (ACT→ACCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (7/29);  new base C (15/18);  total (22/47)
Fisher's exact test for biased strand distribution p-value = 3.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

AGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGC  >  W3110S.gb/71043‑71173
                                                                  |                                                                
aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGcaccac                                                                >  1:2204042/1‑69 (MQ=255)
   gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                             <  1:3539495/69‑1 (MQ=255)
   gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                             <  1:218584/69‑1 (MQ=255)
   gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                             <  1:2851766/69‑1 (MQ=255)
   gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                             <  1:2131557/69‑1 (MQ=255)
   gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:2032709/69‑1 (MQ=255)
   gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:1737523/69‑1 (MQ=255)
    aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCgg                                                            <  1:3365628/69‑1 (MQ=255)
    aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCgg                                                            <  1:1244059/69‑1 (MQ=255)
    aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCgg                                                            <  1:846662/69‑1 (MQ=255)
    aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:1194416/68‑1 (MQ=255)
      tAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGAt                                                          >  1:791437/1‑69 (MQ=255)
       aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                             >  1:3427715/1‑65 (MQ=255)
       aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                             >  1:1841666/1‑65 (MQ=255)
       aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATg                                                         >  1:620251/1‑69 (MQ=255)
       aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATg                                                         >  1:2431881/1‑69 (MQ=255)
       aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATg                                                         >  1:2367736/1‑69 (MQ=255)
          gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCt                                                      <  1:878697/69‑1 (MQ=255)
          gTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:2102013/62‑1 (MQ=255)
          gTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:1841284/62‑1 (MQ=255)
          gTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGcc                                                       <  1:2535249/68‑1 (MQ=255)
            cTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTAt                                                    <  1:1342725/69‑1 (MQ=255)
                aGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCt                                                      <  1:3194046/63‑1 (MQ=255)
                aGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCt                                                      <  1:1113384/63‑1 (MQ=255)
                 gCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGcc                                                <  1:3422963/68‑1 (MQ=255)
                 gCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGcc                                                <  1:3626327/68‑1 (MQ=255)
                 gCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGcc                                                <  1:2102458/68‑1 (MQ=255)
                   tGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCAcc                                             <  1:1376136/69‑1 (MQ=255)
                   tGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCAcc                                             <  1:2871150/69‑1 (MQ=255)
                            gACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                             <  1:2737538/44‑1 (MQ=25)
                            gACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATc                                                   >  1:529083/1‑54 (MQ=255)
                            gACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATc                                                   >  1:533498/1‑54 (MQ=255)
                            gACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATc                                                   >  1:1462158/1‑54 (MQ=255)
                            gACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGc                                    >  1:114474/1‑69 (MQ=255)
                            gACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGc                                    <  1:928174/69‑1 (MQ=255)
                              ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:2601075/42‑1 (MQ=255)
                              ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:1187062/42‑1 (MQ=255)
                              ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:1468044/42‑1 (MQ=255)
                              ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:3389160/42‑1 (MQ=255)
                              ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCaa                                  >  1:1712234/1‑69 (MQ=255)
                              ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCaa                                  >  1:1860711/1‑69 (MQ=255)
                               cAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                             <  1:644174/41‑1 (MQ=255)
                                  cGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGtt                              <  1:522450/69‑1 (MQ=255)
                                  cGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGt                               <  1:542056/68‑1 (MQ=255)
                                        cgccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                        <  1:914327/66‑1 (MQ=255)
                                        agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                        <  1:229607/66‑1 (MQ=255)
                                        agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                        <  1:1128633/66‑1 (MQ=255)
                                        agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                        <  1:185540/66‑1 (MQ=255)
                                        aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                        >  1:2163775/1‑69 (MQ=255)
                                        aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                        <  1:214551/69‑1 (MQ=255)
                                        aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                        <  1:1821343/69‑1 (MQ=255)
                                               cGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                        >  1:461466/1‑62 (MQ=255)
                                               cGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                        >  1:448730/1‑62 (MQ=255)
                                               cGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                 <  1:3517811/69‑1 (MQ=255)
                                                gAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCAc                                              >  1:2711602/1‑39 (MQ=255)
                                                 aaGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGc                                    >  1:3398597/1‑48 (MQ=255)
                                                 aaGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGc                                    >  1:1799592/1‑48 (MQ=255)
                                                 aaGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa                      >  1:153736/1‑62 (MQ=255)
                                                 aaGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCaa               >  1:3580042/1‑69 (MQ=255)
                                                     tgctTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAActct           <  1:1908677/69‑1 (MQ=255)
                                                      gctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACg                    <  1:3409293/59‑1 (MQ=255)
                                                      gctTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTa          >  1:3284932/1‑69 (MQ=255)
                                                         tttGAGCACTACCCTGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                 <  1:3596259/59‑1 (MQ=255)
                                                         tttGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                 <  1:2092994/59‑1 (MQ=255)
                                                         tttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAtt        <  1:1412177/68‑1 (MQ=255)
                                                         tttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAt         <  1:2248508/67‑1 (MQ=255)
                                                         tttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAt         <  1:2336487/67‑1 (MQ=255)
                                                          ttGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAActct           <  1:3312734/64‑1 (MQ=255)
                                                           tGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtct     >  1:1586512/1‑69 (MQ=255)
                                                               cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTcgc  <  1:3447188/68‑1 (MQ=255)
                                                                  |                                                                
AGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGC  >  W3110S.gb/71043‑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: