Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,079 C→T 59.8% 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,0790CT59.8% 53.3 / 104.4 87R231R (CGC→CGTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (7/28);  new base T (26/26);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 6.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTT  >  W3110S.gb/71013‑71145
                                                                  |                                                                  
gcgTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGc                                                                    >  1:3694200/1‑67 (MQ=255)
         acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAgg                                                          >  1:3518296/1‑68 (MQ=255)
         acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAgg                                                          >  1:3213951/1‑68 (MQ=255)
         acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAgg                                                          >  1:2472398/1‑68 (MQ=255)
         acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAgg                                                          >  1:2413576/1‑68 (MQ=255)
         acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAgg                                                          >  1:1403421/1‑68 (MQ=255)
         acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAgg                                                          >  1:1880261/1‑68 (MQ=255)
           aTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAAGCGa                                                       <  1:1030032/69‑1 (MQ=255)
            tCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACTTATCAGCCAGGCGaa                                                      >  1:2255124/1‑69 (MQ=255)
            tCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGa                                                       <  1:1065913/68‑1 (MQ=255)
             cTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAg                                                     <  1:2673516/69‑1 (MQ=255)
              ttttCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAgc                                                    >  1:94743/1‑69 (MQ=255)
              ttttCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAgc                                                    >  1:355889/1‑69 (MQ=255)
              ttttCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgc                                                    >  1:2830984/1‑69 (MQ=255)
              ttttCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgc                                                    >  1:2761198/1‑69 (MQ=255)
                    gCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgctgct                                                <  1:310377/67‑1 (MQ=255)
                        gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                          <  1:517316/69‑1 (MQ=255)
                        gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                          <  1:2594540/69‑1 (MQ=255)
                        gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                          <  1:2329699/69‑1 (MQ=255)
                        gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                          <  1:3332765/69‑1 (MQ=255)
                        gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                          <  1:185228/69‑1 (MQ=255)
                         cagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                          <  1:1114633/68‑1 (MQ=255)
                          agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                         >  1:3572072/1‑68 (MQ=255)
                          agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                         >  1:734486/1‑68 (MQ=255)
                          agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                         >  1:355092/1‑68 (MQ=255)
                          agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                         >  1:1662973/1‑68 (MQ=255)
                             tAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAcca                                     >  1:3400531/1‑67 (MQ=255)
                             tAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAcca                                     >  1:1402530/1‑67 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGc                                                         >  1:1772280/1‑48 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGc                                                         <  1:1943593/48‑1 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGc                                                         <  1:3440765/48‑1 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGc                                                         <  1:3062897/48‑1 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGc                                                         <  1:2156438/48‑1 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                      <  1:2587024/51‑1 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                      <  1:1942895/51‑1 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                      <  1:817678/51‑1 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                      <  1:188001/51‑1 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                      <  1:1466269/51‑1 (MQ=255)
                              aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGcaccac                                    >  1:2204042/1‑69 (MQ=255)
                                 gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                 <  1:2851766/69‑1 (MQ=255)
                                 gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                 <  1:3539495/69‑1 (MQ=255)
                                 gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                 <  1:218584/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:846662/69‑1 (MQ=255)
                                  aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCgg                                <  1:1244059/69‑1 (MQ=255)
                                  aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCgg                                <  1:3365628/69‑1 (MQ=255)
                                  aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                 <  1:1194416/68‑1 (MQ=255)
                                  aTTAGCGGCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                      <  1:2418351/47‑1 (MQ=255)
                                    tAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGAt                              >  1:791437/1‑69 (MQ=255)
                                     aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                          >  1:267100/1‑56 (MQ=255)
                                     aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                          >  1:2083882/1‑56 (MQ=255)
                                     aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                 >  1:3427715/1‑65 (MQ=255)
                                     aGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                 >  1:1841666/1‑65 (MQ=255)
                                     aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATg                             >  1:2367736/1‑69 (MQ=255)
                                     aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATg                             >  1:620251/1‑69 (MQ=255)
                                     aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATg                             >  1:2431881/1‑69 (MQ=255)
                                        gTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                         <  1:3093146/54‑1 (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)
                                             aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAc                                       <  1:817229/51‑1 (MQ=255)
                                             aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAc                                       <  1:2819653/51‑1 (MQ=255)
                                              aGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCt                          <  1:1113384/63‑1 (MQ=255)
                                              aGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCt                          <  1:3194046/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:1468044/42‑1 (MQ=255)
                                                            ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                 <  1:1187062/42‑1 (MQ=255)
                                                            ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                 <  1:2601075/42‑1 (MQ=255)
                                                            ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                 <  1:3389160/42‑1 (MQ=255)
                                                            ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCaa      >  1:1860711/1‑69 (MQ=255)
                                                            ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCaa      >  1:1712234/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)
                                                                  |                                                                  
GCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTT  >  W3110S.gb/71013‑71145

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: