Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 66,726 T→C 70.4% 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,7260TC70.4% 1784.4 / inf 1186intergenic (‑194/‑148)araD/araCpseudogene, truncated/ara regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (159/191);  new base C (370/463);  total (530/656)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1204 total aligned reads displayed.

GCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCC  >  CP009273/66659‑66793
                                                                   |                                                                   
    ataatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTgg                                                                 <  1:34069064/68‑1 (MQ=255)
     taatCAATGTGGTCTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                                                <  1:32537554/68‑1 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                                                 >  1:25086927/1‑67 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                                                 >  1:26981342/1‑67 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                                                 >  1:35290956/1‑67 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                                                <  1:19079635/68‑1 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGAGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                                                <  1:6340787/68‑1 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTccc                                                             >  1:11285113/1‑68 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTccc                                                             >  1:19709751/1‑68 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTccc                                                             >  1:13435413/1‑68 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:16558146/1‑68 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:2440526/1‑68 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:18810493/1‑68 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:10735250/1‑68 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCt                                                          <  1:14297958/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCt                                                          <  1:12101622/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:19424027/68‑1 (MQ=255)
            tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCtt                                                         >  1:34945661/1‑68 (MQ=255)
            tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                                         >  1:34712774/1‑68 (MQ=255)
            tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                                         >  1:32120190/1‑68 (MQ=255)
             gttgACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:1171534/65‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCttt                                                        <  1:21106037/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCttt                                                        <  1:4544765/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCttt                                                        <  1:8884900/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:29874872/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:10484347/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:26669987/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:29602592/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:33298175/68‑1 (MQ=255)
              tgGACTTTTCTGACGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTg                                                       >  1:28577298/1‑68 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGt                                                      <  1:36019135/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGt                                                      <  1:10973047/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACTCGTTTTTGTAATGGCCTTGGTCCCGCTTTGt                                                      >  1:303916/1‑68 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:4008319/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:2636191/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:10002049/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:34423707/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:35705073/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:16351998/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:14828780/68‑1 (MQ=255)
                gACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGtt                                                     >  1:27742393/1‑68 (MQ=255)
                 aCTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGtt                                                     >  1:31811405/1‑67 (MQ=255)
                  cTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTAc                                                   >  1:13856793/1‑68 (MQ=255)
                      tCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGaa                                               <  1:22294304/68‑1 (MQ=255)
                      tCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGaa                                               >  1:7780726/1‑68 (MQ=255)
                      tCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGaa                                               >  1:4799146/1‑68 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:23920968/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:12517284/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:14593876/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:12909122/68‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:17627967/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:36698964/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:29284588/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:9268627/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:3772856/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:28186128/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:20831296/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:7491986/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:9519020/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:17934519/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:27385427/67‑1 (MQ=255)
                            gTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCttt                                         <  1:18521709/68‑1 (MQ=255)
                            gTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCttt                                         <  1:25373080/68‑1 (MQ=255)
                            gTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCttt                                         <  1:33671704/68‑1 (MQ=255)
                              gATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTta                                       >  1:22099123/1‑68 (MQ=255)
                                    aGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGc                                 >  1:15912701/1‑68 (MQ=255)
                                    aGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGc                                 >  1:36970213/1‑68 (MQ=255)
                                    aGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGc                                 >  1:5924653/1‑68 (MQ=255)
                                    aGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGc                                 >  1:5190494/1‑68 (MQ=255)
                                     gACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg                                >  1:3141131/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:20271730/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:710935/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:21508178/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:36351177/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:22854919/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACTGAATGCTTTTAATAAGCgg                               >  1:26033624/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:24819078/1‑68 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              <  1:30819785/68‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               <  1:5545821/67‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg                                <  1:17329535/66‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg                                <  1:15543692/66‑1 (MQ=255)
                                        acTTTTTCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg                             <  1:28987821/68‑1 (MQ=255)
                                        acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg                             >  1:18173633/1‑68 (MQ=255)
                                        acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg                             <  1:807428/68‑1 (MQ=255)
                                             tgctACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTAcc                        <  1:13153210/65‑1 (MQ=255)
                                              gctACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCg                       >  1:16937627/3‑68 (MQ=255)
                                              gctACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAAAATGCTTTTAATAAGCGGGGTTACCg                       >  1:20537616/3‑68 (MQ=255)
                                                tACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCgg                      <  1:20096444/67‑1 (MQ=255)
                                                tACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGt                     >  1:23155163/1‑68 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGt                     >  1:30954213/1‑68 (MQ=255)
                                                   gcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTg                  <  1:23443202/68‑1 (MQ=255)
                                                        tttGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTa              >  1:9848362/1‑67 (MQ=255)
                                                        tttGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAg             <  1:8505918/68‑1 (MQ=255)
                                                           gTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCga          >  1:33023845/1‑68 (MQ=255)
                                                           gTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCga          >  1:1507270/1‑68 (MQ=255)
                                                             cATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCgaga        >  1:21732832/1‑68 (MQ=255)
                                                             cATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCgag         <  1:25781163/67‑1 (MQ=255)
                                                               tGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaag      >  1:4319431/1‑68 (MQ=255)
                                                                 gCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaaga     >  1:6886739/1‑67 (MQ=255)
                                                                  cTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAagag    >  1:24141713/1‑67 (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: