Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 130,455 Δ1 bp 11.8% coding (2354/2598 nt) acnB → bifunctional aconitate hydratase 2/2‑methylisocitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273130,4550T.11.8% 75.0 / 8.1 34coding (2354/2598 nt)acnBbifunctional aconitate hydratase 2/2‑methylisocitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base . (2/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

GGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGTCACTTCCGTGCTGCGGGTAAACTGCTGGATGCGCATAAAGGTCAGTTGCCGACCCGCCTGTGGGTGGCACCGCCAACCCGTATGGACGCCGCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTG  >  CP009273/130196‑130718
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                       
cgTGAGAAGATCGACGCCCTGTCTGTCGGTTCCTGCATGACCAACATCCGTCACTTCGGTACTCCGGGTAAACTCCTTGTCGTGCATAAATCTCACTCTCCGACCGGCCTGTTGGTGGTACCGCCAACCCGTATGGACGCCACCAAATTGACCGAAGACGGCAATCACAACGTCTTCGTGAAGAGTGGTGCGCGTAAAGAGATACCTAGCTGTTCCCTGGGTAGGGGTAAGCAGGCGCGGGTGGGGGACTGTGCAAGGGTGGGTTctgcccaag                                                                                                                                                                                                                                                           >  2:127335/2‑266 (MQ=255)
     gaagaTCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGTCACTTCCGTGCTGCGGGTAAACTGCTGGATGCTCATAAAGGTCAGTTGCCGACCCGCCTGTGGGTGGCACCGCCAACCCGGATGGACGCCGCACAGTTGACCTAAGAAGGCTAATACAGCGTCTCAGGTATGAGTGATGCGCGCATCGAGACGCCGGCCTGTTCCCGTGGTAAGGCTAACCAGGGCCGTGGGCGGGGCGGGGCAGTGGTGGTTTNNANctctttccc                                                                                                                                                                                                                                                       >  2:391505/1‑268 (MQ=255)
     gaagaTCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGTCAATTCCGTGCTGCGGGTAAACTGCTGAATGCGCACAAAGGTCACTTGCCGACCCGCCTGTGGGTGGCACCGCCAACCCGTATGTACGCCCCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTATGAGTGGTGCGCGTATCGAGATCCCTGACCGTTCCCTGTGCATGGGCAACCCGGAGCGCGAGGCGGCGGGTGCCNCTTTTTTTTCCCCCTCTAccgtt                                                                                                                                                                                                                                                      >  2:293690/1‑271 (MQ=255)
                                                           tgctgcGGGTAAACTGCTGGATGCGCATAAAGGTCAGTTGCCGACCCGCCTGTGGGTGGCACCGCCAACCCTTATGGACGCCGCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGAGTAACCAGGCGCGTGTGGAGGACGGTGCAACGGTGGTTTCCACCTCTACACGTAACTCCACGACCTGTCAGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCTGAAc                                                                                                                                                                                               >  1:54583/1‑275 (MQ=255)
                                                                                               gTTGCCGGCCCCTCCTGTGGGTGGCCCTGTCGACCCGTATTTAGTCCGCACATTTGACCGATGAAGGCTACGTCCTCGTCTTCCGTAATGTTGGTTCGCGTATCGAGATCCCTGGCTGTGCCCTTTGTATTGGTAACCAAGCGCTGGTGGCGGACGGTGCAACGGTGGTGTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGTCGCGAATGTCTTCCTGGCTGCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGcc                                                                                                                                                             <  2:58558/274‑1 (MQ=255)
                                                                                                  gCCCACCCGCCTGTGTGTGGCACCGCCAACCCGTATGGACGCCGCACAGTTGACCGAAGAAGGCTCCTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTATTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGc                                                                                                                                                        <  1:127341/275‑1 (MQ=255)
                                                                                                        aagggaccggggaggCACCTCCAACCCGTATGGACGCCGCATAGTTGACCGAAGACGGCTACCACAGCGTCTTCGGTAATCGTGGTGCGGGTATCGAGATCCCCGGCTGTTCCCTGTGTAGGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACGTCCCGAGCCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGaa                                                                                                                                                   <  1:247067/266‑1 (MQ=255)
                                                                                                                    tggCACCGCCAACCCGTATGGACGCCGCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACt                                                                                                                                                                                                                              >  1:183832/1‑188 (MQ=255)
                                                                                                                        aCCGCCAACCCGTATGGACGCCGCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGgcg                                                                                                                                                                                                                          <  2:183826/188‑1 (MQ=255)
                                                                                                                         ccGCCAACCCGTATGGACGCCGCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACggt                                                                                                                                                                                                                                                                         <  2:102121/140‑1 (MQ=255)
                                                                                                                          cGCCAACCCGTATGGACGCCGCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTggg                                                                                                                                                                                                                                  >  2:439416/1‑178 (MQ=255)
                                                                                                                              aaCCCGTATGGACGCCGCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACt                                                                                                                                                                                                                              <  1:439429/178‑1 (MQ=255)
                                                                                                                                    tATGGACGCCGCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTGCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAg                                                                                                                       >  1:437396/1‑275 (MQ=255)
                                                                                                                                           gccgcACAGTTGACCGAAGAAGGCTAGTACAGCGCCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCTAATGTCTTCCTGNCTTCTGCGNAACTgg                                                                                                                                                                                            >  1:330592/1‑199 (MQ=255)
                                                                                                                                            ccgcACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCGGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACa                                                                                                               <  1:292247/275‑1 (MQ=255)
                                                                                                                                                    ttGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGtt                                                                                                         <  1:212931/273‑1 (MQ=255)
                                                                                                                                                                                tCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGgcgc                                                                                                                                                                             >  2:487005/1‑177 (MQ=255)
                                                                                                                                                                                 cGGTAAGGGTGGTGCGCGTATCGTGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGc                                                                          <  1:91226/275‑1 (MQ=255)
                                                                                                                                                                                   gTAAGAGTGGTGCGCGTATCGATCTCCCTGTCTGTTCCCTGTGTATGGGTACCCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCtt                                                                        <  1:391518/275‑1 (MQ=255)
                                                                                                                                                                                    tAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTACGGCGCTGAt                                                                                                                                                                         <  1:487018/177‑1 (MQ=255)
                                                                                                                                                                                       gggttgttctcgTATCGAGATCCCTGGCTGTTCCCTGTGTAGGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTtctc                                                                    <  1:293698/265‑1 (MQ=255)
                                                                                                                                                                                             tgCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACGGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAaga                                                                                                                                                 >  1:453485/1‑192 (MQ=255)
                                                                                                                                                                                             tgCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCaa                                                                                                                                                                   >  2:417426/1‑173 (MQ=255)
                                                                                                                                                                                                 cgTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTAc                                                                                                                                             <  2:453472/192‑1 (MQ=255)
                                                                                                                                                                                                 cgTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTg                                                                                                                                                               <  1:417439/173‑1 (MQ=255)
                                                                                                                                                                                                                      tCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGaa                                                                                                                                                                                                                                           >  1:409418/1‑77 (MQ=255)
                                                                                                                                                                                                                      tCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTg                                                                                                                                 >  1:62338/1‑183 (MQ=255)
                                                                                                                                                                                                                          tgtgTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGt                                                                                                                                                                                                                                       <  2:409405/77‑1 (MQ=255)
                                                                                                                                                                                                                          tgtgTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGgcgc                                                                                                                             <  2:62336/183‑1 (MQ=255)
                                                                                                                                                                                                                                      aCCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTgcggcg                                                                                                                                                                              >  1:333984/1‑122 (MQ=255)
                                                                                                                                                                                                                                      aCCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGAACTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCCACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGATTa                      >  1:173270/1‑274 (MQ=255)
                                                                                                                                                                                                                                      aCCAGGCGCGTGTGGCGGACGGTGCAACGGGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGaa                                                                                                                                                   >  1:289527/1‑148 (MQ=255)
                                                                                                                                                                                                                                          ggCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGa                                                                                                                                                                          <  2:333973/122‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCGCGTGTGGCGGACGGTGCAACGGGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGt                                                                                                                                               <  2:289519/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            tGCAACGGTGGTTTCCACCGCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCATGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACCCCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAACACCCGCGCTg  >  1:402204/1‑272 (MQ=255)
                                                                                                                                                                                                                                                                    tggtTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCAGGAAGAGTACCAGAcc                                                                                                                                       >  1:13895/1‑131 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                       
GGTGAGAAGATCGACGAAGTGTTTATCGGTTCCTGCATGACCAACATCGGTCACTTCCGTGCTGCGGGTAAACTGCTGGATGCGCATAAAGGTCAGTTGCCGACCCGCCTGTGGGTGGCACCGCCAACCCGTATGGACGCCGCACAGTTGACCGAAGAAGGCTACTACAGCGTCTTCGGTAAGAGTGGTGCGCGTATCGAGATCCCTGGCTGTTCCCTGTGTATGGGTAACCAGGCGCGTGTGGCGGACGGTGCAACGGTGGTTTCCACCTCTACCCGTAACTTCCCGAACCGTCTGGGTACTGGCGCGAATGTCTTCCTGGCTTCTGCGGAACTGGCGGCTGTTGCGGCGCTGATTGGCAAACTGCCGACGCCGGAAGAGTACCAGACCTACGTGGCGCAGGTAGATAAAACAGCCGTTGATACTTACCGTTATCTGAACTTCAACCAGCTTTCTCAGTACACCGAGAAAGCCGATGGGGTGATTTTCCAGACTGCGGTTTAAAAAGTCAGCGCACGCGCTG  >  CP009273/130196‑130718

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: