Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 70.9% R263R (CGA→CGG araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,1750AG70.9% 167.5 / 84.6 107R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (8/23);  new base G (41/34);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 9.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 107 total aligned reads displayed.

CACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAAT  >  W3110S.gb/71106‑71241
                                                                     |                                                                  
cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                   <  1:231596/71‑1 (MQ=255)
cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                   <  1:2173535/71‑1 (MQ=255)
cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                   <  1:2169245/71‑1 (MQ=255)
 accaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGt                                                                  <  1:2621913/68‑1 (MQ=255)
 accaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGt                                                                  <  1:1134673/68‑1 (MQ=255)
  ccaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                 <  1:1793384/69‑1 (MQ=255)
  ccaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                 <  1:631933/69‑1 (MQ=255)
    aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                               >  1:1865573/1‑71 (MQ=255)
       cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                            <  1:2811968/71‑1 (MQ=255)
       cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                            <  1:1570671/71‑1 (MQ=255)
       cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                             >  1:2466522/1‑70 (MQ=255)
       cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                             >  1:2707422/1‑70 (MQ=255)
         gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                          <  1:2010058/71‑1 (MQ=255)
         gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                          <  1:3456096/71‑1 (MQ=255)
           ttccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:1216545/3‑71 (MQ=255)
           tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:2972060/1‑71 (MQ=255)
           tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:1581377/1‑71 (MQ=255)
           tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                        >  1:3575533/1‑71 (MQ=255)
           tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGGGTATTTAAAAAATg                                                      <  1:252349/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                      <  1:577140/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                      <  1:548966/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                      <  1:1682595/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                      <  1:393741/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                      <  1:759396/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                      <  1:3239291/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                      <  1:3099312/71‑1 (MQ=255)
              cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                     <  1:873487/71‑1 (MQ=255)
              cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                     <  1:1920642/71‑1 (MQ=255)
              cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                     <  1:2259910/71‑1 (MQ=255)
              cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                     <  1:210115/71‑1 (MQ=255)
              cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                     <  1:3155031/71‑1 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                           <  1:2550744/63‑1 (MQ=255)
                 tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                         >  1:1602041/1‑64 (MQ=255)
                 tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                         >  1:2220288/1‑64 (MQ=255)
                 tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGATATTTaaaaa                                                         >  1:530663/1‑64 (MQ=255)
                 tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                  <  1:1344143/71‑1 (MQ=255)
                   gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:3538199/1‑63 (MQ=255)
                   gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:1539473/1‑63 (MQ=255)
                   gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:716014/1‑63 (MQ=255)
                    ccACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCggg                                               <  1:3555250/71‑1 (MQ=255)
                          tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGc                                         >  1:1106527/1‑71 (MQ=255)
                            ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                       >  1:929873/1‑71 (MQ=255)
                            ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCAAGccc                                       >  1:322005/1‑71 (MQ=255)
                               cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:1841393/1‑51 (MQ=255)
                               cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:3563751/1‑51 (MQ=255)
                               cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:420643/1‑51 (MQ=255)
                               cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:513520/1‑51 (MQ=255)
                               cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:1070189/1‑51 (MQ=255)
                                gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:537627/1‑55 (MQ=255)
                                gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                   >  1:2195973/1‑71 (MQ=255)
                                 cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGTGGCCAGCCcgagcg                                  <  1:2121460/71‑1 (MQ=255)
                                 cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                  <  1:1069406/71‑1 (MQ=255)
                                   aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                <  1:2835959/71‑1 (MQ=255)
                                   aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                <  1:3229440/71‑1 (MQ=255)
                                    tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                     <  1:1935529/49‑1 (MQ=255)
                                    tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                     <  1:2333792/49‑1 (MQ=255)
                                    tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                  <  1:3549507/68‑1 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                >  1:13605/1‑68 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                >  1:1805188/1‑68 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                             >  1:1031662/1‑70 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                             >  1:2810099/1‑70 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                             >  1:417619/1‑70 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                             >  1:2448072/1‑70 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                             <  1:1828861/71‑2 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                              >  1:963192/1‑70 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                             >  1:1755813/1‑71 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                             <  1:1693906/70‑1 (MQ=255)
                                          ttttGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                          <  1:168206/38‑1 (MQ=255)
                                          ttttGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                          <  1:3291062/38‑1 (MQ=255)
                                          ttttGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                          <  1:1402970/38‑1 (MQ=255)
                                          ttttGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                          <  1:3542446/38‑1 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGcc                                            >  1:2901102/1‑50 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                       >  1:2721808/1‑71 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                       >  1:531086/1‑71 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                       >  1:2726282/1‑71 (MQ=255)
                                              gACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:3028831/1‑41 (MQ=255)
                                              gACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:1807786/1‑41 (MQ=255)
                                              gACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:3517153/1‑41 (MQ=255)
                                              gACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:2018602/1‑41 (MQ=255)
                                                  aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                  <  1:3441485/54‑1 (MQ=255)
                                                    cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                >  1:3522458/1‑38 (MQ=255)
                                                    cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                >  1:2322833/1‑70 (MQ=255)
                                                    cAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                >  1:2430464/1‑70 (MQ=255)
                                                         ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa           >  1:740830/1‑70 (MQ=255)
                                                         ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa           >  1:1538422/1‑70 (MQ=255)
                                                         ctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaagaa           >  1:3629736/1‑70 (MQ=255)
                                                          tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa           <  1:2381263/69‑1 (MQ=255)
                                                          tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa           <  1:3034855/69‑1 (MQ=255)
                                                            tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa        <  1:683207/70‑1 (MQ=255)
                                                            tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaa         <  1:1150980/69‑1 (MQ=255)
                                                            tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAg       <  1:1129159/71‑1 (MQ=255)
                                                              tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt      >  1:2714032/1‑70 (MQ=255)
                                                               ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     <  1:2980772/70‑1 (MQ=255)
                                                               ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     <  1:1135943/70‑1 (MQ=255)
                                                               ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     <  1:1131590/70‑1 (MQ=255)
                                                               ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     <  1:3310042/70‑1 (MQ=255)
                                                               ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     <  1:944168/70‑1 (MQ=255)
                                                                  gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat  <  1:3630443/70‑1 (MQ=255)
                                                                  gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat  <  1:3161513/70‑1 (MQ=255)
                                                                  gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTCAat  <  1:1790530/70‑1 (MQ=255)
                                                                     |                                                                  
CACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAAT  >  W3110S.gb/71106‑71241

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: