Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 66,726 T→C 86.8% intergenic (‑194/‑148) 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,7260TC86.8% 13250.6 / inf 4780intergenic (‑194/‑148)araD/araCpseudogene, truncated/ara regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (285/344);  new base C (1828/2315);  total (2117/2663)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 4803 total aligned reads displayed.

GCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCC  >  CP009273/66659‑66793
                                                                   |                                                                   
gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCt                                                                     <  1:29045656/68‑1 (MQ=255)
   aataatCAATGTGGACTTTTCTGCCGTGATTATAGCCACTTTTGTTACGCGTGTTTGTCATGGCTTTg                                                                  >  1:35303039/1‑68 (MQ=255)
   aataatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTg                                                                  <  1:37486064/68‑1 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                                                 >  1:23182640/1‑67 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                                                 >  1:21420246/1‑67 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                                                 >  1:13682812/1‑67 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                                                <  1:47806980/68‑1 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                                                <  1:14678365/68‑1 (MQ=255)
      aatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                                                >  1:18208132/1‑67 (MQ=255)
       atCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTc                                                               >  1:23614327/1‑67 (MQ=255)
       atCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTc                                                               >  1:39911946/1‑67 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTccc                                                             >  1:49090294/1‑68 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:17363401/1‑68 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:14169805/1‑68 (MQ=255)
         cAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCg                                                            <  1:46286249/68‑1 (MQ=255)
          aaTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTccc                                                             >  1:8745335/1‑66 (MQ=255)
          aaTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:25815558/1‑66 (MQ=255)
          aaTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:32473156/1‑66 (MQ=255)
           ttgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:45411942/67‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCt                                                          <  1:29498195/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:26282924/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:12737485/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:36944529/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:42836716/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:19223479/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:15234723/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:33451453/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:30877632/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCttt                                                        <  1:50427014/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:31489863/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        >  1:19629032/1‑68 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        >  1:38994846/1‑68 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:48654948/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:50871467/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:51314294/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:16133035/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:7816073/68‑1 (MQ=255)
              tgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTg                                                       >  1:39357207/1‑68 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:43375545/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      >  1:34464728/1‑68 (MQ=255)
                gACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGtt                                                     >  1:11765488/1‑68 (MQ=255)
                 aCTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGtt                                                     <  1:45824025/67‑1 (MQ=255)
                  cTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTAc                                                   >  1:25053547/1‑68 (MQ=255)
                  cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                                                   >  1:3036960/1‑68 (MQ=255)
                    tttCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                                 >  1:49519410/1‑68 (MQ=255)
                      tCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGaa                                               >  1:2987494/1‑68 (MQ=255)
                      tCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGaa                                               >  1:6087229/1‑68 (MQ=255)
                      tCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGaa                                               >  1:6453343/1‑68 (MQ=255)
                         gCCGTGATTATAGCCACTTTTGCTACGCGTTTATGTCATGGCCTTGGTCCCGCGTTGTAACACAAt                                              >  1:24614476/1‑66 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                                              >  1:28114459/1‑66 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                                              >  1:42347373/1‑66 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:41946698/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:15766035/68‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:10324932/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:9613095/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:22264326/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:10000518/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:12228376/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:20073720/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:46770532/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:38488511/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:21832942/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:26764156/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:27235240/67‑1 (MQ=255)
                              gATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTta                                       <  1:43893597/68‑1 (MQ=255)
                              gATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTta                                       <  1:10889118/68‑1 (MQ=255)
                                ttATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaa                                      <  1:31942408/67‑1 (MQ=255)
                                 tataGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaata                                    >  1:22741360/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:9229678/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:20496441/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:27662209/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:25415097/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:1648218/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:41431347/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:44368718/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:40427303/1‑68 (MQ=255)
                                       cacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              <  1:18690511/68‑1 (MQ=255)
                                       cacTTTTGCTTCGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              <  1:28543395/68‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              <  1:33978429/68‑1 (MQ=255)
                                        acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg                             <  1:44916694/68‑1 (MQ=255)
                                        acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg                             <  1:6978131/68‑1 (MQ=255)
                                        acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg                             <  1:40912376/68‑1 (MQ=255)
                                        acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg                             <  1:7356846/68‑1 (MQ=255)
                                         cTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGGGGGGt                            <  1:34851515/68‑1 (MQ=255)
                                         cTTTTGCTACGCGGTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGt                            >  1:49948680/1‑68 (MQ=255)
                                             tgctACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTAcc                        <  1:48243288/65‑1 (MQ=255)
                                              gctACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCg                       >  1:47310360/3‑68 (MQ=255)
                                                tACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGt                     >  1:8270122/1‑68 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCgg                      <  1:1772660/67‑1 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCgg                      <  1:330881/67‑1 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCgg                      <  1:37890097/67‑1 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGt                     >  1:11321516/1‑68 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGt                     >  1:36564628/1‑68 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGt                     >  1:36151775/1‑68 (MQ=255)
                                                     gTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTgg                 <  1:16821700/67‑1 (MQ=255)
                                                        tttGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAg             <  1:5644045/68‑1 (MQ=255)
                                                        tttGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAg             <  1:21040006/68‑1 (MQ=255)
                                                           gTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCga          >  1:24059939/1‑68 (MQ=255)
                                                                ggCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaaga     <  1:13241244/68‑1 (MQ=255)
                                                                ggCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaaga     <  1:35728334/68‑1 (MQ=255)
                                                                   |                                                                   
GCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCC  >  CP009273/66659‑66793

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