Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,434 C→T 61.9% 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,4340CT61.9% 63.8 / 97.6 84N16N (AAC→AATaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (12/20);  new base T (28/24);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGAC  >  W3110S.gb/70373‑70501
                                                             |                                                                   
ggAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCt                                                             <  1:5158589/70‑1 (MQ=255)
ggAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtg                                                            <  1:485638/71‑1 (MQ=255)
   gTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggt                                                          <  1:1708498/70‑1 (MQ=255)
    tATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGc                                                                   <  1:4880477/60‑1 (MQ=255)
      tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggt                                                          <  1:2321280/67‑1 (MQ=255)
      tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCt                                                             <  1:1757416/64‑1 (MQ=255)
      tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                      >  1:694269/1‑71 (MQ=255)
      tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                      >  1:1520712/1‑71 (MQ=255)
      tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                      <  1:1861750/71‑1 (MQ=255)
      tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                      >  1:5196560/1‑71 (MQ=255)
      tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtggtgg                                                        >  1:333786/1‑69 (MQ=255)
      tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCg                                                      >  1:2921465/1‑71 (MQ=255)
       gAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCt                                                             <  1:1296403/63‑1 (MQ=255)
        aaaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCggg                                                    >  1:4939749/1‑71 (MQ=255)
        aaaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCgg                                                     <  1:2664745/70‑1 (MQ=255)
         aaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGt                                                   >  1:523281/1‑71 (MQ=255)
         aaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGt                                                   >  1:2074927/1‑71 (MQ=255)
            gTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                 >  1:5059455/1‑70 (MQ=255)
            gTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTa                                                <  1:4438243/71‑1 (MQ=255)
            gTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTa                                                <  1:4430488/71‑1 (MQ=255)
               tGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACg                                             >  1:1257576/1‑71 (MQ=255)
                ggCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGc                                            >  1:1000238/1‑71 (MQ=255)
                 gCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc                                           <  1:3627432/71‑1 (MQ=255)
                 gCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc                                           <  1:4059219/71‑1 (MQ=255)
                  cTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGcc                                           <  1:1804173/70‑1 (MQ=255)
                     aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                       >  1:612630/1‑71 (MQ=255)
                     aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                       >  1:4495172/1‑71 (MQ=255)
                      aGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCAt                                                               >  1:545384/1‑46 (MQ=255)
                      aGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTg                                      <  1:22471/71‑1 (MQ=255)
                       gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                       >  1:1919652/1‑69 (MQ=255)
                       gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                     <  1:3417829/71‑1 (MQ=255)
                          cAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                 <  1:3724883/56‑1 (MQ=255)
                          cAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGc                                  >  1:5061463/1‑71 (MQ=255)
                           aaaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGc                                  >  1:5087909/1‑70 (MQ=255)
                            aaaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                 >  1:1002112/1‑70 (MQ=255)
                             aaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGc                                                       <  1:3693330/47‑1 (MQ=255)
                              aTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAAc                              <  1:39740/71‑1 (MQ=255)
                              aTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAAc                              <  1:1898983/71‑1 (MQ=255)
                               tGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACg                             <  1:3919788/71‑1 (MQ=255)
                                 aTCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                                                    <  1:3679474/46‑1 (MQ=255)
                                 aTCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGt                           >  1:1603723/1‑71 (MQ=255)
                                 aTCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCg                                                      <  1:4908886/44‑1 (MQ=255)
                                     ccTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtc                       >  1:4627188/1‑71 (MQ=255)
                                     ccTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtc                       <  1:1848351/71‑1 (MQ=255)
                                      ctgctgCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGtt                          <  1:2551350/67‑1 (MQ=255)
                                      ctgctgCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtct                      >  1:1932510/1‑71 (MQ=255)
                                        gctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                     >  1:4594021/1‑70 (MQ=255)
                                        gctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                     >  1:4426136/1‑70 (MQ=255)
                                         ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                     >  1:3875568/1‑53 (MQ=255)
                                         ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                     >  1:4314434/1‑53 (MQ=255)
                                         ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                     <  1:340057/69‑1 (MQ=255)
                                         ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                     <  1:333848/69‑1 (MQ=255)
                                         ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                     <  1:3193917/69‑1 (MQ=255)
                                         ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGa                   >  1:180391/1‑71 (MQ=255)
                                            ccGGGATACTCGTTTAATGCCCATCTGGTGGCTGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttt                <  1:2812569/71‑1 (MQ=255)
                                             cGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttt               <  1:169437/71‑1 (MQ=255)
                                             cGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttt               <  1:3088058/71‑1 (MQ=255)
                                               ggATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                 <  1:581021/35‑1 (MQ=255)
                                               ggATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttttt             >  1:4297159/1‑71 (MQ=255)
                                               ggATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGcc                                           <  1:1370067/41‑1 (MQ=255)
                                                gATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttttttt            <  1:2360256/71‑1 (MQ=255)
                                                 aTACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTa           <  1:809761/71‑1 (MQ=255)
                                                  tACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttt                <  1:1787826/65‑1 (MQ=255)
                                                    cTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                       >  1:4772061/1‑40 (MQ=255)
                                                    cTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcg        <  1:4750430/71‑1 (MQ=255)
                                                    cTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcg        <  1:1890101/71‑1 (MQ=255)
                                                     tCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                     <  1:1746880/41‑1 (MQ=255)
                                                     tCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                 <  1:402727/45‑1 (MQ=255)
                                                     tCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGtt                          <  1:609290/52‑1 (MQ=255)
                                                     tCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGtt                          <  1:3609714/52‑1 (MQ=255)
                                                      cGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgac      >  1:4554429/1‑71 (MQ=255)
                                                      cGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgac      >  1:4857901/1‑71 (MQ=255)
                                                       gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                               >  1:4312791/1‑45 (MQ=255)
                                                       gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                               >  1:4911040/1‑45 (MQ=255)
                                                       gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc     >  1:1267152/1‑71 (MQ=255)
                                                       gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcga       <  1:4585152/69‑1 (MQ=255)
                                                       gTTTAACGCCCATCTGGAGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAt                        >  1:1658288/1‑52 (MQ=255)
                                                        tttAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccg    <  1:1017277/71‑1 (MQ=255)
                                                         ttAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccga   >  1:4421762/1‑71 (MQ=255)
                                                         ttAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccga   >  1:339392/1‑71 (MQ=255)
                                                          tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                     >  1:2340548/1‑36 (MQ=255)
                                                          tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                     >  1:899766/1‑36 (MQ=255)
                                                          tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccgac  >  1:2516136/1‑71 (MQ=255)
                                                          tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccgac  >  1:4623983/1‑71 (MQ=255)
                                                             |                                                                   
GGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGAC  >  W3110S.gb/70373‑70501

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: