Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 66,706 T→C 88.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,7060TC88.1% 16750.8 / inf 5822intergenic (‑174/‑168)araD/araCpseudogene, truncated/ara regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (380/312);  new base C (2736/2375);  total (3128/2694)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 5831 total aligned reads displayed.

ATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCG  >  CP009273/66639‑66772
                                                                   |                                                                  
    cAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                 >  1:10392401/1‑67 (MQ=255)
    cAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                 >  1:6119694/1‑67 (MQ=255)
     aaaGTGTGACGCCGTGCAAATCATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               <  1:47277960/68‑1 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:27135561/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:42563517/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:14854666/1‑68 (MQ=255)
       aGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGt                                                             <  1:45044807/68‑1 (MQ=255)
       aGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGt                                                             >  1:19745307/1‑68 (MQ=255)
       aGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGt                                                             <  1:16170591/68‑1 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGCCTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                           >  1:40184712/1‑68 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                           >  1:6484401/1‑68 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                           >  1:43950465/1‑68 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                           >  1:26165722/1‑68 (MQ=255)
         tgtgACGCCGTGAAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                           >  1:16517257/1‑68 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGtttt                                                          >  1:157950/1‑68 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtttt                                                          >  1:17910208/1‑68 (MQ=255)
             aCGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTg                                                        >  1:22968893/1‑67 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGt                                                       <  1:38191143/67‑1 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGt                                                       <  1:24432427/67‑1 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTc                                                      >  1:46311459/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCa                                                     >  1:11911677/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:39179972/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:28942327/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:35387383/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:25738127/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:42119253/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:30644307/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:50162901/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:34096489/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:15316839/1‑68 (MQ=255)
                  gTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATgg                                                  <  1:16949616/68‑1 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGc                                                 <  1:128700/67‑1 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGcc                                                <  1:20138635/68‑2 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGcc                                                <  1:345583/68‑2 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGc                                                 <  1:20625866/67‑1 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGc                                                 <  1:49157643/67‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                            >  1:22520667/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                            >  1:51036012/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                           <  1:44539582/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                           <  1:50554215/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                           <  1:31141111/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                           <  1:36303719/68‑1 (MQ=255)
                           atCAATGTGGACTTTTCTGCGGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTc                                          >  1:8795051/1‑67 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:48762843/1‑68 (MQ=255)
                              aaTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:40687670/1‑66 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:37755190/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:41202805/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:33126860/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:11373557/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:25329988/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGACGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:21565768/68‑1 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:7484208/1‑68 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:12320747/1‑68 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:37162510/1‑68 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:27532215/1‑68 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:32554908/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:9258825/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:10902806/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:4778640/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:18804308/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:23887010/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGt                                 <  1:7961306/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:14438185/68‑1 (MQ=255)
                                    gACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:28497168/67‑1 (MQ=255)
                                     aCTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGtt                                >  1:13382809/1‑67 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                              >  1:28028837/1‑68 (MQ=255)
                                        tttCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAg                            <  1:38746071/68‑1 (MQ=255)
                                        tttCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAg                            >  1:31616170/1‑68 (MQ=255)
                                          tCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGaa                          >  1:1748960/1‑68 (MQ=255)
                                            tGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAAt                         >  1:35816255/1‑67 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGc                       <  1:22046145/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       >  1:48282559/1‑68 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:4299717/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:5717669/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:24856198/68‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:21088839/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                     >  1:10000518/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                     >  1:9613095/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                     >  1:29400722/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:8318771/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:33514633/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:7074982/67‑1 (MQ=255)
                                                     tataGACACTTTTGATACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTtaata               <  1:39541519/68‑1 (MQ=255)
                                                      ataGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaataa              <  1:19422540/68‑1 (MQ=255)
                                                         gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg           >  1:26711940/1‑68 (MQ=255)
                                                         gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg           >  1:41678177/1‑68 (MQ=255)
                                                         gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg           >  1:23437174/1‑68 (MQ=255)
                                                          acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:5236616/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:4352738/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:34928431/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:32102817/1‑68 (MQ=255)
                                                           cacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg         >  1:1383756/1‑68 (MQ=255)
                                                           cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTATTAAGCggg         >  1:36741204/1‑68 (MQ=255)
                                                           cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg         <  1:45879465/68‑1 (MQ=255)
                                                            acTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:46778767/68‑1 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        >  1:49741649/1‑68 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:45495514/68‑1 (MQ=255)
                                                             cTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGt       <  1:18274299/68‑1 (MQ=255)
                                                             cTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGt       <  1:29793848/68‑1 (MQ=255)
                                                               tttGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGtt      >  1:302122/1‑67 (MQ=255)
                                                                   |                                                                  
ATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCG  >  CP009273/66639‑66772

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: