Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 66,706 T→C 86.2% 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,7060TC86.2% 10441.7 / inf 3865intergenic (‑174/‑168)araD/araCpseudogene, truncated/ara regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (322/212);  new base C (1807/1516);  total (2135/1730)
Fisher's exact test for biased strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

Only 100 of 3872 total aligned reads displayed.

ATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCG  >  CP009273/66639‑66772
                                                                   |                                                                  
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:10478124/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:5741471/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               >  1:32243675/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                               <  1:27870185/68‑1 (MQ=255)
       aGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGt                                                             >  1:18945206/1‑68 (MQ=255)
       aGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGt                                                             <  1:35813613/68‑1 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                           >  1:37928050/1‑68 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTATGCCGTGAATATAGACACTTTTGCTACGCGttt                                                           >  1:1222390/1‑68 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtttt                                                          >  1:20757010/1‑68 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtt                                                            <  1:13981620/66‑1 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtt                                                            <  1:1508847/66‑1 (MQ=255)
             aCGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTg                                                        >  1:24380935/1‑67 (MQ=255)
             aCGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTg                                                        >  1:37144671/1‑67 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGt                                                       <  1:22151632/67‑1 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGt                                                       <  1:16028769/67‑1 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCa                                                     >  1:23994840/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCa                                                     >  1:30038287/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:6616946/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:17337848/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:15351774/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:35486612/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:13153875/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                     >  1:11875162/1‑68 (MQ=255)
                 cGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATg                                                   >  1:6355140/1‑68 (MQ=255)
                  gTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATgg                                                  <  1:36294545/68‑1 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGcc                                                <  1:23625026/68‑2 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGcc                                                <  1:38529668/68‑2 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGc                                                 <  1:36006590/67‑1 (MQ=255)
                    gCAAATAAACAAGGTGGACTTTTCTGCCGTGATTATAGAGACTTTTGCTACGCGTTTTTGTCATgg                                                  >  1:7338287/1‑66 (MQ=255)
                      aaaTAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCctt                                              >  1:21844949/1‑68 (MQ=255)
                       aataatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTg                                             >  1:27550936/1‑68 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGTCTTTgg                                            >  1:16438753/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTgg                                            >  1:38828168/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                            >  1:37587425/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                            >  1:29051801/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                            >  1:15707860/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                           <  1:17973563/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                           <  1:34540712/68‑1 (MQ=255)
                           atCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTc                                          >  1:3422665/1‑67 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:2041916/1‑68 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:272618/1‑68 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:33881317/1‑68 (MQ=255)
                              aaTGTGGACGTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                        >  1:26930797/1‑66 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATTGACACTTTTGCTACGGGTTTTTGTCATGGCCTTGGTCCctct                                     <  1:34792390/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCt                                     <  1:30467299/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCt                                     <  1:22467469/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:31510496/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     >  1:30815287/1‑68 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:7742594/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                     <  1:29342506/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:5489151/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:9293997/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   <  1:38269856/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTATGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                   >  1:18597123/1‑68 (MQ=255)
                                   ggAGTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGt                                 >  1:8571741/1‑68 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGt                                 <  1:23204298/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:19658494/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:4696293/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                 <  1:6995845/68‑1 (MQ=255)
                                     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGtt                                <  1:10047538/67‑1 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGATAc                              >  1:11570755/1‑68 (MQ=255)
                                           cTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAAt                         <  1:31913370/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGc                       <  1:2969796/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                         >  1:211444/1‑66 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:36745571/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:81273/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:14751009/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                       <  1:14323373/68‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGcc                      <  1:8849453/67‑2 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                     >  1:12714147/1‑68 (MQ=255)
                                                    ttATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaat                >  1:28442019/1‑68 (MQ=255)
                                                    ttATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaat                >  1:19947477/1‑68 (MQ=255)
                                                     tataGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTtaata               <  1:22872045/68‑1 (MQ=255)
                                                        aGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGc            >  1:26189679/1‑68 (MQ=255)
                                                        aGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGc            >  1:19275170/1‑68 (MQ=255)
                                                        aGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGc            >  1:21574776/1‑68 (MQ=255)
                                                         gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg           >  1:9705280/1‑68 (MQ=255)
                                                         gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg           >  1:5166038/1‑68 (MQ=255)
                                                          acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:4051150/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:33432371/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:13595701/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:16710873/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:17032287/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:17651015/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:6111707/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:24739504/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:26639767/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:31152270/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg          >  1:35145272/1‑68 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:25472434/68‑1 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:28157470/68‑1 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:2877586/68‑1 (MQ=255)
                                                            acTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgggg        <  1:10836129/68‑1 (MQ=255)
                                                             cTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGt       <  1:18304563/68‑1 (MQ=255)
                                                             cTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGt       <  1:32684599/68‑1 (MQ=255)
                                                             cTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGt       <  1:25113398/68‑1 (MQ=255)
                                                             cTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGt       >  1:4299968/1‑68 (MQ=255)
                                                             cTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGt       <  1:25835689/68‑1 (MQ=255)
                                                             cTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGt       <  1:33077520/68‑1 (MQ=255)
                                                               tttGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGtt      >  1:11193480/1‑67 (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: