Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 66,706 T→C 91.8% 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,7060TC91.8% 20225.6 / inf 6430intergenic (‑174/‑168)araD/araCpseudogene, truncated/ara regulon transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (297/228);  new base C (3155/2727);  total (3464/2966)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

Only 100 of 6439 total aligned reads displayed.

ATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGG  >  CP009273/66639‑66773
                                                                   |                                                                   
    cAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                  >  1:32847280/1‑67 (MQ=255)
    cAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAc                                                                  >  1:36007787/1‑67 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                                >  1:32478388/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                                >  1:40828991/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                                >  1:7250261/1‑68 (MQ=255)
     aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgc                                                                >  1:20335376/1‑68 (MQ=255)
      aaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTAcgcg                                                               <  1:29476878/68‑1 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                            >  1:31376317/1‑68 (MQ=255)
         tgtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGttt                                                            >  1:23624588/1‑68 (MQ=255)
          gtgACGCCGTGCTAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtt                                                             <  1:8413708/66‑1 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtttt                                                           >  1:6410649/1‑68 (MQ=255)
          gtgACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGtt                                                             <  1:42392183/66‑1 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGt                                                        <  1:26472569/67‑1 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGt                                                        <  1:41953584/67‑1 (MQ=255)
              cGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGt                                                        <  1:31011001/67‑1 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                      >  1:38358745/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                      >  1:24034086/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                      >  1:2768385/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                      >  1:43134293/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                      >  1:9583555/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                      >  1:10003308/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                      >  1:8846401/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCa                                                      >  1:12737914/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACCCGTTTTTGTCa                                                      >  1:10826220/1‑68 (MQ=255)
               gCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTCGCTACGCGTTTTTGTCa                                                      >  1:4153520/1‑68 (MQ=255)
                 cGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATg                                                    >  1:16459385/1‑68 (MQ=255)
                  gTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATgg                                                   <  1:603816/68‑1 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGcc                                                 <  1:14435099/68‑2 (MQ=255)
                    gCAAATAATCAATGTGGACTTTTCTGCCGGGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGcc                                                 <  1:15311188/68‑2 (MQ=255)
                        ataatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                             <  1:39991290/68‑1 (MQ=255)
                        ataatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                             <  1:14878135/68‑1 (MQ=255)
                        ataatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                             <  1:19848984/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGTTTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                            <  1:8021679/68‑1 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                             >  1:3752045/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                             >  1:13163670/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTgg                                             >  1:1712444/1‑67 (MQ=255)
                         taatCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGt                                            <  1:28366106/68‑1 (MQ=255)
                           atCAATGTGGACTTGTCTGCGGTGATTATAGAATCTTTTGCTACGCGTTTTTGTCATGGCCTTGGTc                                           >  1:20691861/1‑67 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                         >  1:4585863/1‑68 (MQ=255)
                            tCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                         >  1:40388292/1‑68 (MQ=255)
                              aaTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTccc                                         >  1:9243336/1‑66 (MQ=255)
                               aTGTGTACGTTTCTGACGGGATTCTAGACACTTTTTTTACGCGGTGTTTTCATGGCTTTTGTCCCGCt                                      <  1:19498095/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                      <  1:10398272/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                      <  1:26894822/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                      <  1:11236248/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                      <  1:17488506/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                      <  1:2725233/68‑1 (MQ=255)
                               aTGTGGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCt                                      <  1:26034070/68‑1 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                     >  1:24802663/1‑68 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                     >  1:37943668/1‑68 (MQ=255)
                                tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCtt                                     >  1:35583361/1‑68 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCttt                                    <  1:42789696/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                    <  1:38806066/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                    <  1:15715494/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                    <  1:43587262/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                    <  1:17963272/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                    <  1:18693203/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                    >  1:39625736/1‑68 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                    <  1:22260779/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCttt                                    <  1:30270456/68‑1 (MQ=255)
                                 gtgGACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGGTTTTGGCATGGCCTTGGTCCCGCttt                                    <  1:37115703/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGt                                  <  1:6830616/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGt                                  <  1:2562756/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCTCGCTTTGt                                  <  1:34001683/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                  <  1:36386022/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                  <  1:25196997/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                  <  1:21478544/68‑1 (MQ=255)
                                   ggACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGt                                  >  1:16096802/1‑68 (MQ=255)
                                    gACTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGtt                                 >  1:18292830/1‑68 (MQ=255)
                                     aCTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGtt                                 >  1:31776883/1‑67 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                               >  1:36762530/1‑68 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                               <  1:28005315/68‑1 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                               >  1:33594431/1‑68 (MQ=255)
                                      cTTTTCTGCCGTGATTATAGACACTTTTGCTACGAGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTAc                               >  1:29934573/1‑68 (MQ=255)
                                        tttCTGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAg                             >  1:5249060/1‑68 (MQ=255)
                                            tGCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATg                         <  1:28695809/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAAt                          >  1:5671832/1‑66 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                        <  1:19068880/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                        <  1:21899217/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                        <  1:16826121/68‑1 (MQ=255)
                                             gCCGTGATTATAGACACTTTTGCTACGCGGTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGc                        <  1:35154929/68‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                      >  1:11920656/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                      >  1:34403049/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCtt                      >  1:34755481/1‑68 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                       <  1:3318405/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                       <  1:1403372/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                       <  1:24436396/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                       <  1:30645636/67‑1 (MQ=255)
                                               cGTGATTATAGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCt                       <  1:21019353/67‑1 (MQ=255)
                                                        aGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGc             >  1:7663465/1‑68 (MQ=255)
                                                        aGACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCACGCTTTGTTACAGAATGCTTTTAATAAGc             >  1:29129313/1‑68 (MQ=255)
                                                         gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg            >  1:23128188/1‑68 (MQ=255)
                                                         gACACTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCg            >  1:13627356/1‑68 (MQ=255)
                                                          acacTTTTGCTTCGCGTTTTTGTCATGGCCTTGGTCCCGCCTTGTTACAGAATGCTTTTAATAAGCgg           >  1:41143560/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTa                                >  1:11572917/1‑47 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg           >  1:32109669/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg           >  1:12288498/1‑68 (MQ=255)
                                                          acacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg           >  1:39251763/1‑68 (MQ=255)
                                                           cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCggg          <  1:4849399/68‑1 (MQ=255)
                                                           cacTTTTGCTACGCGTTTTTGTCATGGCCTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCgg           <  1:2263881/67‑1 (MQ=255)
                                                                   |                                                                   
ATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGG  >  CP009273/66639‑66773

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: