Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,434 C→T 71.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,4340CT71.9% 45.6 / 22.2 32N16N (AAC→AATaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (17/6);  total (21/11)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCG  >  W3110S.gb/70366‑70499
                                                                    |                                                                 
aaTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAAc                                                                   <  1:551852/69‑1 (MQ=255)
aaTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGc                                                                 >  1:1230304/1‑71 (MQ=255)
           tATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtggtgg                                                      <  1:1100684/71‑1 (MQ=255)
            aTGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGc                                                     >  1:2180845/1‑71 (MQ=255)
            aTGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGc                                                     >  1:1390863/1‑71 (MQ=255)
            aTGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTCGTGGc                                                     >  1:2179527/1‑71 (MQ=255)
            aTGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtggtgg                                                      <  1:1491693/70‑1 (MQ=255)
            aTGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGc                                                     >  1:1341310/1‑71 (MQ=255)
             tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGTGATACTCGTTTAACGCCCATCTGGTGGCg                                                    >  1:595618/1‑71 (MQ=255)
             tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                    >  1:83790/1‑71 (MQ=255)
             tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                    >  1:247296/1‑71 (MQ=255)
               aaaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATc                                                            >  1:1915028/1‑61 (MQ=255)
               aaaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                                                  >  1:1513699/1‑71 (MQ=255)
                aaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGt                                                 >  1:1165171/1‑71 (MQ=255)
                aaaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGt                                                 >  1:1201171/1‑71 (MQ=255)
                  aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                               <  1:1916460/71‑1 (MQ=255)
                  aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                               <  1:211318/71‑1 (MQ=255)
                    tATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                               <  1:1417847/69‑1 (MQ=255)
                             aGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTa                                              <  1:1865662/61‑1 (MQ=255)
                              gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                     >  1:488676/1‑69 (MQ=255)
                                 cAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCa                                                              <  1:1221921/41‑1 (MQ=255)
                                   aaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCAGATTGAGGcc                               >  1:499204/1‑70 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCAt                                                             >  1:68726/1‑39 (MQ=255)
                                     aTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAAc                            >  1:1833781/1‑71 (MQ=255)
                                      tGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACg                           >  1:840951/1‑71 (MQ=255)
                                           cccTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAt                      <  1:1193293/71‑1 (MQ=255)
                                           cccTGATGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                               <  1:493441/46‑1 (MQ=255)
                                                  gcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtt               >  1:1875707/1‑71 (MQ=255)
                                                       gATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttttt           >  1:831035/1‑70 (MQ=255)
                                                              gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc   >  1:1029882/1‑71 (MQ=255)
                                                              gTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc   <  1:103319/71‑1 (MQ=255)
                                                               tttAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccg  <  1:1344657/71‑1 (MQ=255)
                                                                 tAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcga     >  1:351536/1‑66 (MQ=255)
                                                                    |                                                                 
AATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCG  >  W3110S.gb/70366‑70499

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: