Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 66,706 T→C 63.1% intergenic (‑174/‑168) araD ← / → araC pseudogene, truncated/ara regulon transcriptional activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927366,7060TC63.1% 1006.7 / inf 1099intergenic (‑174/‑168)araD/araCpseudogene, truncated/ara regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (230/174);  new base C (383/307);  total (617/482)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

Only 100 of 1100 total aligned reads displayed.

ATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCG  >  CP009273/66639‑66772
                                                                   |                                                                  
    cAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTAc                                                                 >  1:645748/1‑67 (MQ=255)
    cAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                 >  1:81481/1‑67 (MQ=255)
    cAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                 >  1:18024369/1‑67 (MQ=255)
    cAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                 >  1:25819201/1‑67 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTAcgc                                                               >  1:20325648/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTAcgc                                                               >  1:12667815/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTAcgc                                                               >  1:22839434/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:16358199/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:20145759/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:10106855/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               <  1:29229702/68‑1 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               <  1:21038202/68‑1 (MQ=255)
       aGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGt                                                             >  1:8814617/1‑68 (MQ=255)
        gtgtGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtt                                                            >  1:31471155/1‑68 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGttt                                                           >  1:11750553/1‑68 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                           >  1:5822072/1‑68 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGtttt                                                          >  1:13936023/1‑68 (MQ=255)
            gACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTg                                                        <  1:9390321/68‑1 (MQ=255)
            gACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTg                                                        >  1:32106194/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCa                                                     >  1:17125750/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCa                                                     >  1:8505364/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCa                                                     >  1:30725039/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCa                                                     >  1:10448895/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTATCa                                                     >  1:21458607/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGAGTTTTTGTCa                                                     >  1:1421918/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:2739707/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:18248064/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:10676501/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:25326472/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:15668767/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:9115724/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:14780792/1‑68 (MQ=255)
               gCCGTGCAAATAATCAAGGTGGACTTTTCTGCCGTGAATATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:22183502/1‑68 (MQ=255)
                 cGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATg                                                   >  1:26534836/1‑68 (MQ=255)
                 cGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATg                                                   <  1:3967907/68‑1 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCTTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCt                                                <  1:14554506/68‑1 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGcc                                                <  1:27183025/68‑2 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTgg                                            >  1:28275700/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTgg                                            >  1:6866618/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGt                                           <  1:6324613/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                            >  1:13644147/1‑67 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:21873771/1‑68 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:7918900/1‑68 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:7554848/1‑68 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCt                                     <  1:4174453/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCt                                     <  1:22430407/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGGGTTTTTGTCATGGCCTTGGTCGCGCt                                     <  1:26012119/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:20774037/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:16607411/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:20534306/68‑1 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:31768981/1‑68 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:2353593/1‑68 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:5968441/1‑68 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCttt                                   <  1:26841474/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:16814318/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:11023512/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:11355743/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:27730994/68‑1 (MQ=255)
                                   tgACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:30249770/67‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGt                                 <  1:19727619/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:13336046/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:25036048/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:4480915/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:24131025/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:17795552/68‑1 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTAc                              >  1:27921080/1‑68 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                              >  1:28703925/1‑68 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                              >  1:9503083/1‑68 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                              >  1:19991849/1‑68 (MQ=255)
                                        tttCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                            >  1:15022475/1‑68 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGc                       <  1:24477520/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGc                       <  1:4982757/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:548040/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:4784194/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:18835820/68‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:15994103/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                     >  1:19084873/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:9631776/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:25566759/67‑1 (MQ=255)
                                                gTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCttt                    <  1:29608143/68‑1 (MQ=255)
                                                 tGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtttt                   >  1:3480072/1‑68 (MQ=255)
                                                  gATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTta                  <  1:23894740/68‑1 (MQ=255)
                                                    ttATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTtaat                >  1:15454905/1‑68 (MQ=255)
                                                         gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg           >  1:19492478/1‑68 (MQ=255)
                                                          acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:5215801/1‑68 (MQ=255)
                                                          acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:31216658/1‑68 (MQ=255)
                                                          acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:24696568/1‑68 (MQ=255)
                                                          acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:12216610/1‑68 (MQ=255)
                                                          acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:13017066/1‑68 (MQ=255)
                                                          acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:17524580/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:30523253/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:29016014/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:26377051/1‑68 (MQ=255)
                                                           cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg         >  1:7204101/1‑68 (MQ=255)
                                                           cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg         >  1:23254268/1‑68 (MQ=255)
                                                           cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTATGTTACAGAATGCTTTTAATAAGCggg         <  1:18614368/68‑1 (MQ=255)
                                                            acTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:23809358/68‑1 (MQ=255)
                                                              ttttGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGtt      >  1:29952392/1‑68 (MQ=255)
                                                                  gTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCg  >  1:366809/1‑68 (MQ=255)
                                                                  gTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCg  >  1:11948566/1‑68 (MQ=255)
                                                                   |                                                                  
ATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCG  >  CP009273/66639‑66772

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: