Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 66,706 T→C 90.5% 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,7060TC90.5% 25938.2 / inf 8500intergenic (‑174/‑168)araD/araCpseudogene, truncated/ara regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (449/361);  new base C (4115/3564);  total (4572/3928)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

Only 100 of 8512 total aligned reads displayed.

ATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCG  >  CP009273/66639‑66772
                                                                   |                                                                  
   gCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                 >  1:17855969/1‑68 (MQ=255)
     aaaGTGTGGCGCCGTGCAAATAATCAATGTGGACTTTTCGGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:54437087/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACTCTTTTGCTAcgc                                                               >  1:56645667/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTAcgc                                                               >  1:48099358/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:29614096/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               <  1:22454808/68‑1 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:40105509/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:18885944/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:24960831/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:50255189/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                 <  1:42765262/66‑1 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                 <  1:11056031/66‑1 (MQ=255)
       aGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGt                                                             <  1:51289896/68‑1 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                           >  1:1995517/1‑68 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtttt                                                          >  1:5966173/1‑68 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtt                                                            <  1:14233215/66‑1 (MQ=255)
          gtgACGCCGTGCAAATAATCAAAGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtt                                                            <  1:21994974/66‑1 (MQ=255)
             aCGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTg                                                        >  1:51818777/1‑67 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGt                                                       <  1:10509726/67‑1 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTc                                                      >  1:13205233/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:27970264/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:36513069/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:35345124/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:45484558/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:47161327/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:655942/1‑68 (MQ=255)
                  gTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATgg                                                  <  1:46558186/68‑1 (MQ=255)
                  gTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATgg                                                  <  1:38604151/68‑1 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGc                                                 <  1:49140157/67‑1 (MQ=255)
                      aaaTAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGcct                                               <  1:24459096/67‑3 (MQ=255)
                       aataatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTg                                             <  1:8893173/68‑1 (MQ=255)
                        ataatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                            <  1:21000620/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTgg                                            >  1:9484405/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                            >  1:46050898/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                            >  1:41108001/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                           <  1:10029961/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                           <  1:17427492/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                           <  1:38016312/68‑1 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:40634939/1‑68 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:2545280/1‑68 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:35945089/1‑68 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:28486370/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:44959629/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:21494393/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:26921184/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:16353693/68‑1 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:52382277/1‑68 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:56144149/1‑68 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                    >  1:37500670/1‑68 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCttt                                   <  1:32873734/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:41595432/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:22882313/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:23450432/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:8507692/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:43256850/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:57194546/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:30724701/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:36973083/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:4864041/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:12611901/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:7534831/68‑1 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                              >  1:30157916/1‑68 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                              >  1:34813061/1‑68 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                         >  1:26453164/1‑66 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                         >  1:18392428/1‑66 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                         >  1:16879198/1‑66 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:52886138/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:14729816/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:20457778/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:49610934/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:47567311/68‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCtt                     >  1:53456226/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                     >  1:57669797/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                     >  1:31845430/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:50721810/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:44421473/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:19465942/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:33816280/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:39561406/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGGCCCGCTTTGTTACAGAATGCt                      <  1:43918644/67‑1 (MQ=255)
                                                gTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCttt                    <  1:29053955/68‑1 (MQ=255)
                                                  gATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTta                  <  1:39100318/68‑1 (MQ=255)
                                                    ttATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaat                >  1:7084065/1‑68 (MQ=255)
                                                    ttATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaat                >  1:42151381/1‑68 (MQ=255)
                                                     tataGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaata               <  1:54949835/68‑1 (MQ=255)
                                                      ataGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaataa              <  1:13732588/68‑1 (MQ=255)
                                                       taGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAg             <  1:15847030/68‑1 (MQ=255)
                                                         gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg           >  1:802384/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:3336894/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:12122997/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:15258671/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:34374431/1‑68 (MQ=255)
                                                           cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg           <  1:32351037/66‑1 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:31257746/68‑1 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:5392090/68‑1 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:23925878/68‑1 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:11550271/68‑1 (MQ=255)
                                                              ttttGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                      <  1:25912841/52‑1 (MQ=255)
                                                               tttGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGtt      >  1:55534933/1‑67 (MQ=255)
                                                               tttGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTa     <  1:27422375/68‑1 (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: