Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,434 C→T 38.1% N16N (AAC→AAT araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb70,4340CT38.1% 56.3 / 69.1 63N16N (AAC→AATaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (16/23);  new base T (11/13);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

AATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGC  >  W3110S.gb/70366‑70504
                                                                    |                                                                      
aaTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGc                                                                      >  1:252931/1‑71 (MQ=255)
 aTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGc                                                                      <  1:2298602/70‑1 (MQ=255)
    gtggGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCa                                                                   >  1:3150319/1‑70 (MQ=255)
       ggAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCt                                                                <  1:383796/70‑1 (MQ=255)
           tATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtggtgg                                                           >  1:1194988/1‑71 (MQ=255)
             tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGc                                                          <  1:1779261/70‑1 (MQ=255)
               aaaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtgg                                                              >  1:281975/1‑64 (MQ=255)
                aaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGt                                                      >  1:1854082/1‑71 (MQ=255)
                  aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                    <  1:1121674/71‑1 (MQ=255)
                  aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                    <  1:1040243/71‑1 (MQ=255)
                  aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                    <  1:2124218/71‑1 (MQ=255)
                    tATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGc                                                          >  1:1102912/1‑63 (MQ=255)
                    tATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGtt                                                     <  1:1200808/68‑1 (MQ=255)
                      tGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAAc                                                 <  1:2430297/70‑1 (MQ=255)
                        gCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGcc                                              <  1:2662713/71‑1 (MQ=255)
                         cTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCg                                             <  1:2124653/71‑1 (MQ=255)
                            aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                          <  1:1069058/71‑1 (MQ=255)
                              gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                        <  1:2787426/71‑1 (MQ=255)
                              gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                        <  1:391591/71‑1 (MQ=255)
                               cgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAg                                       <  1:286232/71‑1 (MQ=255)
                                 cAAAATGATCACCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGtt                                                     >  1:2266532/1‑55 (MQ=255)
                                  aaaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                        <  1:2239459/67‑1 (MQ=255)
                                  aaaaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                    <  1:1089878/55‑1 (MQ=255)
                                  aaaaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                    <  1:2055268/71‑1 (MQ=255)
                                   aaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtgg                                                              >  1:1834394/1‑44 (MQ=255)
                                   aaaTGATCCCCTGCTGCCGGGATACTCGTTTAACGccc                                                                    >  1:564493/1‑38 (MQ=255)
                                   aaaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCt                                                                >  1:385582/1‑42 (MQ=255)
                                   aaaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCt                                                                <  1:1074569/42‑1 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCt                                                                <  1:759930/41‑1 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                                  >  1:2961299/1‑71 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                                  >  1:1097036/1‑71 (MQ=255)
                                        aTCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGt                              <  1:1675408/71‑1 (MQ=255)
                                        aTCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACgg                               <  1:3075354/70‑1 (MQ=255)
                                            ccTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                    <  1:1143203/45‑1 (MQ=255)
                                             ctgctgCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGcc                                              <  1:2463401/50‑1 (MQ=255)
                                             ctgctgCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtct                         >  1:511091/1‑71 (MQ=255)
                                              tgctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAAc                                                 <  1:2372288/46‑1 (MQ=255)
                                              tgctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAAc                                                 <  1:151837/46‑1 (MQ=255)
                                                ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGa                      <  1:2462703/71‑1 (MQ=255)
                                                ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGa                      <  1:2458966/71‑1 (MQ=255)
                                                ctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                    >  1:1414535/1‑41 (MQ=255)
                                                     gggATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTTAGGCCAACGGTTATCTCGAttttt                 <  1:2657283/71‑1 (MQ=255)
                                                     gggATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttt                  >  1:1702319/1‑70 (MQ=255)
                                                      ggATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttttt                <  1:2347416/71‑1 (MQ=255)
                                                          aCTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                          <  1:54704/41‑1 (MQ=255)
                                                           cTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcg           <  1:1197633/71‑1 (MQ=255)
                                                              gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:2886230/1‑71 (MQ=255)
                                                              gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:36367/1‑71 (MQ=255)
                                                              gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:174714/1‑71 (MQ=255)
                                                              gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGAt                                           >  1:79791/1‑36 (MQ=255)
                                                               tttAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTa              <  1:730297/64‑1 (MQ=255)
                                                               tttAACGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                          <  1:54399/36‑1 (MQ=255)
                                                               tttAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttt                   <  1:2474390/59‑1 (MQ=255)
                                                                ttAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGc                                     <  1:671132/40‑1 (MQ=255)
                                                                 tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAAc                                 >  1:298946/1‑43 (MQ=255)
                                                                 tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttttt                 >  1:1781094/1‑59 (MQ=255)
                                                                 tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttttt                 >  1:1784806/1‑59 (MQ=255)
                                                                 tAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                                  <  1:411191/42‑1 (MQ=255)
                                                                 tAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttttt                 >  1:828615/1‑59 (MQ=255)
                                                                 tAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAta                    >  1:2550379/1‑55 (MQ=255)
                                                                 tAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATc            >  1:3005958/1‑64 (MQ=255)
                                                                  aaCGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttttt                 >  1:1753103/1‑58 (MQ=255)
                                                                   aCGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACg                                >  1:1962277/1‑42 (MQ=255)
                                                                    cGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGc  <  1:2315691/71‑1 (MQ=255)
                                                                    |                                                                      
AATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGC  >  W3110S.gb/70366‑70504

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: