Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,434 C→T 55.4% 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,4340CT55.4% 17.7 / 104.8 74N16N (AAC→AATaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (19/14);  new base T (23/18);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGA  >  W3110S.gb/70367‑70500
                                                                   |                                                                  
aTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGcc                                                                 >  1:412119/1‑71 (MQ=255)
    gtggAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATc                                                             <  1:4471665/69‑1 (MQ=255)
      ggAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCt                                                            <  1:4180171/70‑1 (MQ=255)
           aTGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGc                                                      >  1:36160/1‑71 (MQ=255)
            tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                     >  1:4707868/1‑71 (MQ=255)
            tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                     >  1:3885540/1‑71 (MQ=255)
            tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtggtgg                                                       >  1:219469/1‑69 (MQ=255)
            tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCg                                                     >  1:1305601/1‑71 (MQ=255)
            tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCg                                                     <  1:691473/71‑1 (MQ=255)
             gAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCgg                                                    >  1:2892915/1‑71 (MQ=255)
              aaaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCgg                                                    <  1:2210358/70‑1 (MQ=255)
               aaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGt                                                  >  1:3637196/1‑71 (MQ=255)
                 aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCgg                                                    >  1:2788686/1‑67 (MQ=255)
                 aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCgg                                                    >  1:3093666/1‑67 (MQ=255)
                 aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                <  1:3358440/71‑1 (MQ=255)
                 aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                <  1:2595054/71‑1 (MQ=255)
                   tATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCt                                                            <  1:4468589/57‑1 (MQ=255)
                   tATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTaa                                              <  1:3114492/71‑1 (MQ=255)
                    aTGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa                                              >  1:1923782/1‑70 (MQ=255)
                     tGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtg                                                           <  1:208326/56‑1 (MQ=255)
                     tGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtg                                                           <  1:2263199/56‑1 (MQ=255)
                     tGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtg                                                           <  1:74214/56‑1 (MQ=255)
                     tGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCggg                                                   >  1:3276521/1‑64 (MQ=255)
                     tGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCggg                                                   >  1:3691129/1‑64 (MQ=255)
                     tGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                >  1:3533925/1‑67 (MQ=255)
                       gCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtgg                                                          >  1:472749/1‑55 (MQ=255)
                         tGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGa                                        >  1:1918482/1‑71 (MQ=255)
                         tGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCt                                                            >  1:610052/1‑51 (MQ=255)
                          gAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAt                                       >  1:4131361/1‑71 (MQ=255)
                          gAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAt                                       >  1:2213646/1‑71 (MQ=255)
                          gAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGccc                                                                <  1:4025547/46‑1 (MQ=255)
                           aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGa                                        <  1:2186294/69‑1 (MQ=255)
                             gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCt                                                            <  1:1019478/47‑1 (MQ=255)
                             gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                <  1:4703120/59‑1 (MQ=255)
                               gcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtgg                                                          >  1:3215360/1‑47 (MQ=255)
                               gcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtgg                                                          >  1:4224615/1‑47 (MQ=255)
                                cAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                <  1:2973251/56‑1 (MQ=255)
                                 aaaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                <  1:1661005/71‑1 (MQ=255)
                                     tGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACg                            >  1:710166/1‑71 (MQ=255)
                                          cccTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                              >  1:2180888/1‑64 (MQ=255)
                                          cccTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAAAGGTTAt                       <  1:622975/71‑1 (MQ=255)
                                             tgctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACgg                           <  1:4623985/64‑1 (MQ=255)
                                             tgctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACgg                           <  1:1136172/64‑1 (MQ=255)
                                             tgctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                    >  1:1914061/1‑71 (MQ=255)
                                             tgctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                    >  1:2786727/1‑71 (MQ=255)
                                              gctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAgg                                  >  1:4260863/1‑56 (MQ=255)
                                              gctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                    >  1:870977/1‑70 (MQ=255)
                                              gctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGa                                        >  1:3793968/1‑50 (MQ=255)
                                               ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGa                                        >  1:4264316/1‑49 (MQ=255)
                                               ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                    <  1:33975/69‑1 (MQ=255)
                                               ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGa                  >  1:783316/1‑71 (MQ=255)
                                               ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGa                  <  1:1357669/71‑1 (MQ=255)
                                               ctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTa                                               >  1:2991794/1‑42 (MQ=255)
                                               ctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTa                                               >  1:1576269/1‑42 (MQ=255)
                                               ctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTa                                               >  1:1587572/1‑42 (MQ=255)
                                               ctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTa                                               >  1:1969240/1‑42 (MQ=255)
                                                     ggATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttttt            <  1:3435535/71‑1 (MQ=255)
                                                     ggATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttttt            <  1:1324605/71‑1 (MQ=255)
                                                      gATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttttttt           <  1:4026226/71‑1 (MQ=255)
                                                        tACTCGTTTAATGCCCACCTGGTGGCGGGTTTAACGCCGAtt                                      <  1:4614660/42‑1 (MQ=38)
                                                         aCTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGtt                         >  1:1199026/1‑54 (MQ=255)
                                                          cTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcg       <  1:2253535/71‑1 (MQ=255)
                                                          cTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcg       <  1:1860543/71‑1 (MQ=255)
                                                           tCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                    >  1:3445541/1‑41 (MQ=255)
                                                           tCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                    >  1:1379421/1‑41 (MQ=255)
                                                            cGTTTAACGCCCATCTGGTGGCGGGTTTAACGCAGATTGAGGCCAACGGTTATCTCGATTTTTTTATcga      >  1:4270729/1‑70 (MQ=255)
                                                             gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc    <  1:4584007/71‑1 (MQ=255)
                                                             gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc    <  1:2700277/71‑1 (MQ=255)
                                                             gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc    >  1:2909591/1‑71 (MQ=255)
                                                             gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgac     >  1:1314874/1‑70 (MQ=255)
                                                              tttAATGCCCATCTGGTGGCGGGTTTAACGCCGATTg                                     <  1:4239391/37‑1 (MQ=255)
                                                              tttAACGCCCATCTGGTGGCGGGTTTAACGCCGATTg                                     <  1:2957278/37‑1 (MQ=255)
                                                                tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAg                                   >  1:1312686/1‑37 (MQ=255)
                                                                tAATGCACATCTGGTGGCGGGTTTAACGCCGATTGAg                                   >  1:10769/1‑37 (MQ=38)
                                                                   cGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttt              <  1:4373058/55‑1 (MQ=255)
                                                                   cGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccga  >  1:2791515/1‑67 (MQ=255)
                                                                   |                                                                  
ATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGA  >  W3110S.gb/70367‑70500

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: