Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,434 C→T 66.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,4340CT66.1% 60.1 / 55.2 56N16N (AAC→AATaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (14/5);  new base T (20/17);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGC  >  W3110S.gb/70368‑70504
                                                                  |                                                                      
tggtggGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGccc                                                                    <  1:1128307/71‑1 (MQ=255)
      gAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtg                                                               >  1:965908/1‑70 (MQ=255)
        gTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtggtg                                                            >  1:863144/1‑71 (MQ=255)
           tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                         >  1:886201/1‑71 (MQ=255)
           tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCg                                                         <  1:1254773/71‑1 (MQ=255)
             aaaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTgggggg                                                       >  1:931080/1‑71 (MQ=255)
             aaaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCggg                                                       >  1:832891/1‑71 (MQ=255)
             aaaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCggg                                                       >  1:577480/1‑71 (MQ=255)
                aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                    <  1:273790/71‑1 (MQ=255)
                aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                    <  1:1083410/71‑1 (MQ=255)
                aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                    <  1:395626/71‑1 (MQ=255)
                aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                    <  1:430941/71‑1 (MQ=255)
                aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGCGGCGGGttt                                                    <  1:574591/71‑1 (MQ=255)
                aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                    <  1:1248616/71‑1 (MQ=255)
                aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                    <  1:14634/71‑1 (MQ=255)
                aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                    <  1:1339050/71‑1 (MQ=255)
                      gCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc                                              >  1:754034/1‑71 (MQ=255)
                          aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                          <  1:718072/71‑1 (MQ=255)
                          aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                          <  1:1028492/71‑1 (MQ=255)
                          aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                          <  1:997640/71‑1 (MQ=255)
                          aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                          <  1:1558531/71‑1 (MQ=255)
                                aaaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                    <  1:713035/71‑1 (MQ=255)
                                aaaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                    <  1:1680614/71‑1 (MQ=255)
                                aaaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                    <  1:119004/71‑1 (MQ=255)
                                 aaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                    <  1:1620000/70‑1 (MQ=255)
                                   aTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtgg                                                           >  1:1661533/1‑45 (MQ=255)
                                   aTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtgg                                                           >  1:1442977/1‑45 (MQ=255)
                                   aTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtgg                                                           >  1:190600/1‑45 (MQ=255)
                                   aTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAAc                                 >  1:182434/1‑71 (MQ=255)
                                    tGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACg                                >  1:1121521/1‑71 (MQ=255)
                                       tCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGt                              >  1:1329963/1‑70 (MQ=255)
                                            tgctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                    <  1:118599/59‑1 (MQ=255)
                                            tgctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                    <  1:611374/59‑1 (MQ=255)
                                             gctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                        >  1:868106/1‑70 (MQ=255)
                                             gctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                        >  1:1257528/1‑70 (MQ=255)
                                             gctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                        >  1:67462/1‑70 (MQ=255)
                                              ctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCg                       <  1:584353/70‑1 (MQ=255)
                                                     gATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttttttt               <  1:1441319/71‑1 (MQ=255)
                                                         cTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                        >  1:1088161/1‑42 (MQ=255)
                                                         cTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                        >  1:752082/1‑42 (MQ=255)
                                                         cTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                        >  1:968992/1‑42 (MQ=255)
                                                         cTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                        >  1:654796/1‑42 (MQ=255)
                                                         cTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATc            >  1:588405/1‑70 (MQ=255)
                                                          tCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcga          >  1:793021/1‑71 (MQ=255)
                                                          tCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcga          >  1:1087270/1‑71 (MQ=255)
                                                            gTTTAACGCCCTTCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:832042/1‑71 (MQ=255)
                                                            gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:341638/1‑71 (MQ=255)
                                                            gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:1008116/1‑71 (MQ=255)
                                                            gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:950416/1‑71 (MQ=255)
                                                            gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:939448/1‑71 (MQ=255)
                                                            gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:1769961/1‑71 (MQ=255)
                                                            gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:458494/1‑71 (MQ=255)
                                                            gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc        >  1:1315524/1‑71 (MQ=255)
                                                               tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccgac     >  1:408999/1‑71 (MQ=255)
                                                               tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccgac     >  1:807207/1‑71 (MQ=255)
                                                               tAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaacgac     >  1:1797170/1‑71 (MQ=255)
                                                                  cGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGc  <  1:851164/71‑1 (MQ=255)
                                                                  cGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGc  <  1:104067/71‑1 (MQ=255)
                                                                  cGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGc  <  1:1757812/71‑1 (MQ=255)
                                                                  |                                                                      
TGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGC  >  W3110S.gb/70368‑70504

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: