Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 70,434 C→T 42.7% 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
*AP00904870,4340CT42.7% 43.6 / 102.4 82N16N (AAC→AATaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (27/20);  new base T (14/21);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

AATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACC  >  AP009048/70366‑70502
                                                                    |                                                                    
aaTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGc                                                                    >  1:2283946/1‑71 (MQ=255)
 aTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGc                                                                    <  1:2180579/70‑1 (MQ=255)
 aTGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCACCTGCTGCCGGGATACTCGTTTAATGcc                                                                   >  1:522322/1‑71 (MQ=255)
    gtggGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCa                                                                 <  1:3249741/70‑1 (MQ=255)
    gtggGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCAt                                                                <  1:2101431/71‑1 (MQ=255)
     tggGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCa                                                                 <  1:3904768/69‑1 (MQ=255)
       ggAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtg                                                             <  1:493711/71‑1 (MQ=255)
          gTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtg                                                          <  1:2384438/71‑1 (MQ=255)
          gTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtggtg                                                          <  1:1158026/71‑1 (MQ=255)
          gTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtggtg                                                          >  1:4054739/1‑71 (MQ=255)
             tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                       >  1:1395365/1‑71 (MQ=255)
             tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCg                                                       >  1:4023209/1‑71 (MQ=255)
             tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGccc                                                                  <  1:1132016/60‑1 (MQ=255)
             tGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGc                                                        <  1:866319/70‑1 (MQ=255)
                 aaGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGtt                                                   >  1:2811259/1‑71 (MQ=255)
                  aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGttt                                                  <  1:1077351/71‑1 (MQ=255)
                  aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                  <  1:4775204/71‑1 (MQ=255)
                  aGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGttt                                                  <  1:2990345/71‑1 (MQ=255)
                   gTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTa                                                 <  1:3489732/71‑1 (MQ=255)
                    tATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTaa                                                >  1:236288/1‑71 (MQ=255)
                      tGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAc                                               <  1:1897149/70‑1 (MQ=255)
                       ggCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCtgg                                                            >  1:1396935/1‑56 (MQ=255)
                         cTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCg                                           <  1:5142933/71‑1 (MQ=255)
                           gAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc                                            <  1:4224755/68‑1 (MQ=255)
                            aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtg                                                          >  1:3155446/1‑53 (MQ=255)
                            aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                        <  1:3535216/71‑1 (MQ=255)
                            aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTa                                                 >  1:3583334/1‑62 (MQ=255)
                            aaGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                        <  1:4045807/71‑1 (MQ=255)
                             aGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtggtg                                                          <  1:290270/52‑1 (MQ=255)
                             aGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTg                                       <  1:4225072/71‑1 (MQ=255)
                             aGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGAGGGtt                                                   >  1:406984/1‑59 (MQ=255)
                             aGCGCAAAATGATCCCCAGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGAtt                                        >  1:1054951/1‑70 (MQ=255)
                              gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGtt                                                   >  1:1365859/1‑58 (MQ=255)
                              gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                      <  1:1251271/71‑1 (MQ=255)
                              gcgcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                      <  1:4876878/71‑1 (MQ=255)
                                gcAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                      <  1:309071/69‑1 (MQ=255)
                                 cAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTg                                       <  1:2462344/67‑1 (MQ=255)
                                  aaaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCggg                                                     >  1:4695114/1‑52 (MQ=255)
                                   aaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                  >  1:4890467/1‑70 (MQ=255)
                                   aaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                  <  1:1209196/70‑1 (MQ=255)
                                   aaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGcc                                  >  1:1923972/1‑70 (MQ=255)
                                   aaaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCa                                 <  1:545355/71‑1 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTaa                                                >  1:3653907/1‑55 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                                >  1:1518644/1‑71 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGAt                                         >  1:5088380/1‑62 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                                >  1:515131/1‑71 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                                >  1:3467896/1‑71 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCaa                                >  1:5135061/1‑71 (MQ=255)
                                    aaTGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCa                                 <  1:738988/70‑1 (MQ=255)
                                        aTCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc                                            <  1:3324292/55‑1 (MQ=255)
                                        aTCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACgg                             <  1:3575477/70‑1 (MQ=255)
                                        aTCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGt                            <  1:865804/71‑1 (MQ=255)
                                         tCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCa                                 <  1:3154783/65‑1 (MQ=255)
                                         tCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                      >  1:2617725/1‑60 (MQ=255)
                                         tCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGtt                           <  1:607536/71‑1 (MQ=255)
                                          ccccTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTa                          >  1:1783926/1‑71 (MQ=255)
                                            ccTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtc                        >  1:2462522/1‑71 (MQ=255)
                                             ctgctgCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtct                       >  1:88140/1‑71 (MQ=255)
                                             ctgctgCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtct                       >  1:1398741/1‑71 (MQ=255)
                                              tgctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                      >  1:1487345/1‑71 (MQ=255)
                                               gctgcCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                      >  1:321040/1‑70 (MQ=255)
                                                ctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCg                                           >  1:2941339/1‑48 (MQ=255)
                                                   ccGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttt                 <  1:1790302/71‑1 (MQ=255)
                                                     gggATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttttt               <  1:4503818/71‑1 (MQ=255)
                                                     gggATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttttt               <  1:1275607/71‑1 (MQ=255)
                                                        aTACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGa                                      >  1:1945039/1‑45 (MQ=255)
                                                         tACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTAt           >  1:5164111/1‑71 (MQ=255)
                                                         tACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTAt           >  1:5170344/1‑71 (MQ=255)
                                                           cTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGa                                          >  1:971502/1‑38 (MQ=255)
                                                           cTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATc          <  1:2707239/70‑1 (MQ=255)
                                                           cTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATc          <  1:4461275/70‑1 (MQ=255)
                                                           cTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATc          <  1:231098/70‑1 (MQ=255)
                                                              gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAtctc                      >  1:1443788/1‑55 (MQ=255)
                                                              gTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc      >  1:226810/1‑71 (MQ=255)
                                                               tttAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAtttt                <  1:2829779/60‑1 (MQ=255)
                                                               tttAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGAttt                 <  1:1820614/59‑1 (MQ=255)
                                                               tttAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgacc      >  1:1023301/1‑70 (MQ=255)
                                                                ttAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTa            >  1:3891218/1‑63 (MQ=255)
                                                                 tAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccgac   >  1:5315536/1‑71 (MQ=255)
                                                                 tAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccgac   >  1:2649946/1‑71 (MQ=255)
                                                                 tAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATcgaccga    <  1:1404334/70‑1 (MQ=255)
                                                                   aCGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCgaccgacc  >  1:346021/1‑70 (MQ=255)
                                                                    |                                                                    
AATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACC  >  AP009048/70366‑70502

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: