Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,085 T→C 57.8% 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.8% 54.0 / 110.8 90S233S (AGT→AGCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (8/30);  new base C (26/26);  total (34/56)
Fisher's exact test for biased strand distribution p-value = 7.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTT  >  W3110S.gb/71022‑71151
                                                               |                                                                  
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:1880261/1‑68 (MQ=255)
acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAgg                                                                >  1:1403421/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:3332765/69‑1 (MQ=255)
               gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                <  1:2594540/69‑1 (MQ=255)
               gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                <  1:185228/69‑1 (MQ=255)
               gcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                <  1:2329699/69‑1 (MQ=255)
                cagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                                <  1:1114633/68‑1 (MQ=255)
                 agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                               >  1:1662973/1‑68 (MQ=255)
                 agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                               >  1:355092/1‑68 (MQ=255)
                 agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                               >  1:3572072/1‑68 (MQ=255)
                 agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                               >  1:734486/1‑68 (MQ=255)
                    tAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAcca                                           >  1:3400531/1‑67 (MQ=255)
                    tAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAcca                                           >  1:1402530/1‑67 (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:1772280/1‑48 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGc                                                               <  1:2156438/48‑1 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                            <  1:2587024/51‑1 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                            <  1:1466269/51‑1 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                            <  1:817678/51‑1 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                            <  1:1942895/51‑1 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGaa                                                            <  1:188001/51‑1 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGcaccac                                          >  1:2204042/1‑69 (MQ=255)
                        gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                       <  1:218584/69‑1 (MQ=255)
                        gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                       <  1:3539495/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:1244059/69‑1 (MQ=255)
                         aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCgg                                      <  1:3365628/69‑1 (MQ=255)
                         aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCgg                                      <  1:846662/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:2431881/1‑69 (MQ=255)
                            aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATg                                   >  1:2367736/1‑69 (MQ=255)
                            aGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATg                                   >  1:620251/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:2819653/51‑1 (MQ=255)
                                    aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAc                                             <  1:817229/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:533498/1‑54 (MQ=255)
                                                 gACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATc                             >  1:1462158/1‑54 (MQ=255)
                                                 gACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATc                             >  1:529083/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: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)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt  >  1:2163775/1‑69 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt  <  1:214551/69‑1 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt  <  1:1821343/69‑1 (MQ=255)
                                                               |                                                                  
ACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTT  >  W3110S.gb/71022‑71151

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: