Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,214 T→C 42.1% F276F (TTT→TTC araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,2140TC42.1% 24.0 / 70.9 57F276F (TTT→TTCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (18/15);  new base C (12/12);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGAC  >  W3110S.gb/71158‑71282
                                                        |                                                                    
cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGt                                                             >  1:1839835/1‑66 (MQ=255)
cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtg                                                           >  1:2862844/1‑68 (MQ=255)
cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtg                                                           >  1:2776719/1‑68 (MQ=255)
cAACTCTATTTCTAGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtg                                                           >  1:3132566/1‑68 (MQ=255)
 aaCTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                                                         <  1:2692626/69‑1 (MQ=255)
         ttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaa                                                  >  1:3138376/1‑68 (MQ=255)
           ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAg                                                >  1:444058/1‑68 (MQ=255)
           ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAg                                                >  1:3380708/1‑68 (MQ=255)
             cgcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg                                              <  1:2855071/68‑1 (MQ=255)
             cgcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg                                              <  1:2237035/68‑1 (MQ=255)
                gggTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat                                           <  1:2178218/66‑1 (MQ=255)
                gggTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat                                           <  1:1754136/66‑1 (MQ=255)
                gggTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat                                           <  1:3046040/66‑1 (MQ=255)
                     tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGt                                                                   >  1:2885642/1‑39 (MQ=255)
                       tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAg                                    >  1:2161802/1‑68 (MQ=255)
                       tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGc                                   >  1:3679111/1‑69 (MQ=255)
                       tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGc                                   >  1:1253918/1‑69 (MQ=255)
                       tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGc                                   >  1:2040144/1‑69 (MQ=255)
                       tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGc                                   >  1:3201/1‑69 (MQ=255)
                       tAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGc                                   >  1:2435100/1‑69 (MQ=255)
                        aaaaaaTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGcc                                  >  1:3634448/1‑69 (MQ=255)
                           aaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTa                                     >  1:189144/1‑63 (MQ=255)
                            aaTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCa                              <  1:2209996/69‑1 (MQ=255)
                            aaTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCa                              <  1:83502/69‑1 (MQ=255)
                               gCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGt                           >  1:2682729/1‑69 (MQ=255)
                                    ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg                                          >  1:2858865/1‑49 (MQ=255)
                                    ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg                                          >  1:2163015/1‑49 (MQ=255)
                                    ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg                                          >  1:874358/1‑49 (MQ=255)
                                     gggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                     <  1:2088848/69‑1 (MQ=255)
                                     gggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                     <  1:380218/69‑1 (MQ=255)
                                     gggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                     <  1:2738620/69‑1 (MQ=255)
                                     gggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                     <  1:1429089/69‑1 (MQ=255)
                                     gggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                     <  1:1158006/69‑1 (MQ=255)
                                      ggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTa                                     <  1:482414/52‑1 (MQ=255)
                                       gCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATaa                   <  1:3173596/69‑1 (MQ=255)
                                       gCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATaa                   <  1:2218133/69‑1 (MQ=255)
                                        ccAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAt                  <  1:2805301/69‑1 (MQ=255)
                                         cAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATa                    <  1:159161/66‑1 (MQ=255)
                                         cAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATa                    <  1:1312531/66‑1 (MQ=255)
                                          aGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAtt                 >  1:2711360/1‑68 (MQ=255)
                                          aGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAtt                 >  1:1502104/1‑68 (MQ=255)
                                            cccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGt              <  1:1083644/69‑1 (MQ=255)
                                            cccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGt              <  1:978746/69‑1 (MQ=255)
                                            cccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGt              <  1:1206241/69‑1 (MQ=255)
                                            cccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGt              <  1:1971359/69‑1 (MQ=255)
                                             ccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTa             >  1:2873113/1‑69 (MQ=255)
                                             ccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTa             >  1:1435164/1‑69 (MQ=255)
                                             ccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTa             >  1:771912/1‑69 (MQ=255)
                                             ccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTa             >  1:1980431/1‑69 (MQ=255)
                                             ccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTa             >  1:3708365/1‑69 (MQ=255)
                                             ccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTa             >  1:3704186/1‑69 (MQ=255)
                                             ccGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTa             >  1:3166257/1‑69 (MQ=255)
                                             ccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtc             <  1:228389‑M1/69‑7 (MQ=255)
                                             ccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtc             <  1:3458934‑M1/69‑7 (MQ=255)
                                              cgagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtc             >  1:585498‑M1/1‑62 (MQ=255)
                                                agcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtccct          <  1:1140353‑M1/69‑10 (MQ=255)
                                                     gTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAAagggtccct          <  1:2676897‑M1/64‑10 (MQ=255)
                                                       ttCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                     <  1:339755/51‑1 (MQ=255)
                                                       ttCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCAt                     <  1:787527/51‑1 (MQ=255)
                                                        tCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGt              <  1:3630256/57‑1 (MQ=255)
                                                        tCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGt              <  1:1525696/57‑1 (MQ=255)
                                                        tCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGAc  <  1:635650/69‑1 (MQ=255)
                                                        tCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGAc  <  1:2866315/69‑1 (MQ=255)
                                                        |                                                                    
CAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGAC  >  W3110S.gb/71158‑71282

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: