Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 66,726 T→C 85.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,7260TC85.8% 6343.8 / inf 2369intergenic (‑194/‑148)araD/araCpseudogene, truncated/ara regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (168/168);  new base C (862/1167);  total (1032/1337)
Fisher's exact test for biased strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 2387 total aligned reads displayed.

GCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCC  >  CP009273/66659‑66793
                                                                   |                                                                   
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTgg                                                                 >  1:11491803/1‑67 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                                                <  1:5129710/68‑1 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                                                <  1:11285320/68‑1 (MQ=255)
     taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                                                <  1:21179147/68‑1 (MQ=255)
        tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:8306842/1‑68 (MQ=255)
          aaTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:3122844/1‑66 (MQ=255)
          aaTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                                             >  1:22978841/1‑66 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:15728725/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:14864253/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:19746520/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:26619633/68‑1 (MQ=255)
           aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                                          <  1:12197004/68‑1 (MQ=255)
            tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                                         >  1:18013272/1‑68 (MQ=255)
            tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                                         >  1:17449997/1‑68 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:26359269/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:21482095/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:2771701/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:3994449/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:17205376/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:6148299/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:10234586/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:13784029/68‑1 (MQ=255)
             gtgGACTTTTCTGCCGTGATTAAAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCttt                                                        <  1:23559329/68‑1 (MQ=255)
              tgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:16879116/67‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:9475477/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:7479797/68‑1 (MQ=255)
               ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                                      <  1:17738702/68‑1 (MQ=255)
                 aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGtt                                                     >  1:16585832/1‑67 (MQ=255)
                 aCTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGtt                                                     >  1:22624991/1‑67 (MQ=255)
                  cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                                                   >  1:9734074/1‑68 (MQ=255)
                  cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                                                   >  1:25229687/1‑68 (MQ=255)
                  cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                                                   >  1:3749924/1‑68 (MQ=255)
                   ttttCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACa                                                  <  1:6864805/68‑1 (MQ=255)
                    tttCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAg                                                 >  1:10775285/1‑68 (MQ=255)
                    tttCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAg                                                 >  1:8580192/1‑68 (MQ=255)
                      tCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGaa                                               >  1:19070326/1‑68 (MQ=255)
                        tGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATg                                             <  1:119905/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:19333350/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                                              >  1:18506107/1‑66 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                                              >  1:12477969/1‑66 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                                              >  1:15033960/1‑66 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:18764940/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:20930407/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:13577735/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:28701982/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:7757571/68‑1 (MQ=255)
                         gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                                            <  1:13035846/68‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCtt                                          >  1:8849449/1‑68 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:20546812/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:7156273/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:26934820/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:14026007/67‑1 (MQ=255)
                           cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                                           <  1:15992452/67‑1 (MQ=255)
                            gTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCttt                                         <  1:29252925/68‑1 (MQ=255)
                              gATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtttt                                        <  1:24548462/67‑1 (MQ=255)
                              gATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTta                                       <  1:22376259/68‑1 (MQ=255)
                              gATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTta                                       <  1:14629579/68‑1 (MQ=255)
                                 tataGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                                        <  1:25937390/48‑1 (MQ=255)
                                    aGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGc                                 >  1:27160349/1‑68 (MQ=255)
                                     gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg                                >  1:5830112/1‑68 (MQ=255)
                                     gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg                                >  1:15361292/1‑68 (MQ=255)
                                     gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg                                >  1:11742369/1‑68 (MQ=255)
                                     gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg                                >  1:25423791/1‑68 (MQ=255)
                                     gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg                                >  1:20246963/1‑68 (MQ=255)
                                     gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAAAAAGCg                                >  1:4540554/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:26194542/1‑68 (MQ=255)
                                      acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:28018412/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGGTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:16360433/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:12765787/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:10443130/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:23754688/1‑68 (MQ=255)
                                      acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               >  1:5616862/1‑68 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              <  1:4217980/68‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              <  1:8016367/68‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              <  1:22106566/68‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              <  1:23982504/68‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              >  1:24237837/1‑68 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg                               <  1:20735764/67‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg                                <  1:15536503/66‑1 (MQ=255)
                                       cacTTTTGCTACGCGTTTTTGTCATGGCCGTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg                              <  1:10019690/68‑1 (MQ=255)
                                        acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg                             >  1:4837482/1‑68 (MQ=255)
                                             tgctACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTTcc                        <  1:24861569/65‑1 (MQ=255)
                                              gctACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCg                       >  1:13266624/3‑68 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCgg                      <  1:28390767/67‑1 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCgg                      <  1:92526/67‑1 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCgg                      <  1:6342691/67‑1 (MQ=255)
                                                tACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGt                     >  1:19939704/1‑68 (MQ=255)
                                                   gcgTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTg                  <  1:23276814/68‑1 (MQ=255)
                                                     gTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTgg                 <  1:28922515/67‑1 (MQ=255)
                                                     gTTTTTGGCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTgg                 <  1:14256486/67‑1 (MQ=255)
                                                        tttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAg             <  1:6593092/68‑1 (MQ=255)
                                                        tttGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAg             <  1:19587405/68‑1 (MQ=255)
                                                        tttGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAg             <  1:18173470/68‑1 (MQ=255)
                                                           gTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCGTTTAATAAGCGGGGGTACCGGTTTGGTTAGCgg          >  1:5409266/1‑67 (MQ=255)
                                                           gTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCga          >  1:25672970/1‑68 (MQ=255)
                                                             cATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCgag         <  1:27454786/67‑1 (MQ=255)
                                                             cATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCgaga        >  1:3438759/1‑68 (MQ=255)
                                                             cATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCgag         <  1:9056608/67‑1 (MQ=255)
                                                                ggCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaaga     <  1:1103942/68‑1 (MQ=255)
                                                                ggCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaaga     <  1:21810577/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: