Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,717 T→C 13.4% intergenic (+145/‑295) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,7170TC13.4% 96.0 / 3.0 60intergenic (+145/‑295)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base T (27/19);  new base C (5/4);  total (34/26)
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.00e+00

CGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACG  >  NC_000913/349577‑349858
                                                                                                                                            |                                                                                                                                             
cGCCCGTTGATGGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGa                                                                                                                                             <  2:620922/143‑1 (MQ=255)
 gCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGAc                                                                                                                                            >  1:531010/1‑143 (MQ=255)
   ccGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGGCGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACtt                                                                                                                                          <  2:311564/143‑1 (MQ=255)
   ccGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGCTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACtt                                                                                                                                          >  1:540096/1‑143 (MQ=255)
     gTTGATTACATACGCGAGCCGACGCCTGATGCGTCAATGACGCGACTTATAACGCCTCCGGGGTGCTCTGACGTTTAGGTCGGTTAAGACGGCTGGCGACGCGTCCCACAACCGATGCCTGATGCGCCGCTGACGTGACTTAt                                                                                                                                        <  2:271300/143‑1 (MQ=9)
             tATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGCACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGACGCCCGGAGCGCCGCCGCCGCGACTCAGCAGcccgc                                                                                                                                >  2:240049/1‑141 (MQ=14)
             tATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGGCGTAGCAGGCCTa                                                                                                                                >  2:633380/1‑143 (MQ=14)
              aTTCGACAGCCGATGCCTGATGCGTCGCTCACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGCCAACCGATCCCCCATGCGCCGCTCACCTCCCTTACCCGGccgcc                                                                                                                               >  1:1037093/1‑140 (MQ=14)
               ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGCCCTACGAGGTGCACGGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGCTGCCTGATGCCCCGCTGACGTCACCTGTCAGGccgccg                                                                                                                              >  1:368046/1‑139 (MQ=14)
                tCGCCAGCCGAGGCCTGATGCGTCGCTGACGCGACTTATCCGGCCTACGAGGGGCACTGAGCTGTAGGGTGGCTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACgg                                                                                                                             <  2:266687/143‑1 (MQ=14)
                 cGACAGCCGATGCCTGATGCTTCTCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTTGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                            <  1:122705/143‑1 (MQ=255)
                 cGACAGCCGATGCCTGATGCGTCGCTGGCGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                            <  2:1073768/143‑1 (MQ=255)
                    cAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTg                                                                                                                         <  1:1014052/143‑1 (MQ=255)
                       ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCAc                                                                                                                      >  1:934107/1‑143 (MQ=255)
                       ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCACCGACGCGAATGATCAGGCCTCCGGCGTGcga                                                                                                                      >  2:351689/1‑141 (MQ=11)
                        cGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGTATGGCGTCGCATCCGACAACAGTTGCCTGATGCGCCGCTGACGTGACTGATCAGGCCTAGgggggacaca                                                                                                                     >  1:75201/1‑137 (MQ=14)
                        cGATGCCTGATGCGTCGCTGACGCGACATATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACt                                                                                                                     >  2:651921/1‑143 (MQ=255)
                           tgacggATGCGTCGCTGACGCGACTGATCAGGCCTCTGGGGCGCACTGAACTGGAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGaa                                                                                                                  <  2:470828/138‑1 (MQ=14)
                           taccTGATGCGTCGCTGACGCTACTTACCAGGCCTACGACGTGCACTGAACTGGAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGTTGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGaa                                                                                                                  <  1:256022/141‑1 (MQ=14)
                           tGCCTGATGCGTCGCTCCCGCGACTTATCAGGCCTACGGGGTGCTCTGAACTGTAGGTCGGATAAGACGGATGTCGTCGCATCCGACAACCGATGCCTGATGCGCCACTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                                                                                                                  <  2:878006/143‑1 (MQ=11)
                             ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGAgcactgcact                                                                                                                >  1:121507/1‑143 (MQ=255)
                             ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCCGGCCTACGGGGCGCCCCGAACt                                                                                                                >  1:762255/1‑143 (MQ=14)
                               tGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCAATGAACTgc                                                                                                              >  1:937842/1‑142 (MQ=255)
                                 aTGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGGGCACTCGACTGtcg                                                                                                            >  1:165889/1‑141 (MQ=255)
                                      tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                                                                                       >  1:480274/1‑143 (MQ=255)
                                              gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTCGGTCGGATAAGGCg                                                                                               >  1:99450/1‑143 (MQ=255)
                                                         ggCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCAATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCCCTGAACTGTAGGTCGGCTAACACCCACAGCGTCGc                                                                                    >  1:248201/1‑143 (MQ=14)
                                                          gCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTGGGATAAGACGCACAGCGTCGCa                                                                                   >  1:419883/1‑143 (MQ=255)
                                                          gCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACAGTAGGTCGGATAAGACGCATAGCGTCGCa                                                                                   >  1:538605/1‑143 (MQ=255)
                                                             tACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATCGCGTCGCATcc                                                                                >  1:124727/1‑143 (MQ=255)
                                                               cGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                                                              >  2:595054/1‑143 (MQ=25)
                                                                 aGGTGCACTGAGCTGGAGGGCGGAGAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACa                                                                            <  2:909206/143‑1 (MQ=255)
                                                                     gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCg                                                                        <  2:1035490/143‑1 (MQ=14)
                                                                       gtggagcTGTAGGGCGGATAAGCCGGATGGCGTCGTATCCGACAGCCGATGCCTGATGCGCCGCTGACGCAGCTTATGAGGCCTACGGGGTTCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                                                                      <  2:377223/137‑1 (MQ=11)
                                                                         tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                    >  1:1050160/1‑143 (MQ=38)
                                                                           aaCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                                                                  <  2:566198/143‑1 (MQ=17)
                                                                            aCTGTAGGTCGGATAAGACGGAGGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                                                                 <  2:47246/143‑1 (MQ=37)
                                                                              tGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                                                               >  1:171547/1‑143 (MQ=38)
                                                                                 aGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTAGCAGGGATAAGGGGTGCGCTGAGATGTAGGGCGGGTGAGAGGGATGGCGTCGGATCCGACAGCCGATGCCTGATgcg                                                            >  2:177386/1‑143 (MQ=14)
                                                                                   gTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                                                                                                          >  2:330043/1‑79 (MQ=255)
                                                                                       gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg                                                      <  1:628888/143‑1 (MQ=37)
                                                                                            gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCGTCCGACAACCGATGCCTGGTGCGCCGCTGAcgcg                                                 >  1:75283/1‑143 (MQ=255)
                                                                                                  tGGCGTCACATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACGGAACTGGAGGTCGGATAAGACGCGTAGCGGCGCATCCCACAACCGATGCCTAAAACACCGCGGACGCGAcagat                                           >  2:958274/1‑139 (MQ=11)
                                                                                                    gcatCGCACGCGACAACAGATGCCTGATGCGCCGCTGCCGTGACTTCTCACGCCACCGGTGTGCACTGGACTGTGGGTAGGATAAGCCGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                         <  1:17537/140‑1 (MQ=11)
                                                                                                     cGTCGCGTCCGCCAACCGATGCCTGATGCGCCGCTGACGTGACTTATCGGGCCTACGGGGTGCACTGAACTGTAGTTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAg                                        <  2:476760/143‑1 (MQ=255)
                                                                                                              cGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTAAAGCGCCGCAGACGCGCCTTATCAGGCCTACggg                               >  2:574854/1‑143 (MQ=255)
                                                                                                               gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAGATGCGTGATGCGCCGCTGACGCGAGTTATGAGGCCTAAgggg                              >  1:748203/1‑143 (MQ=11)
                                                                                                               gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGCATAAGACGCATACCGTCGCATCCGCCAACCGCTACCAGCAGCGCGACTGACGCGACCTCTCAGACCAAcaagg                              >  2:872183/1‑139 (MQ=11)
                                                                                                                 cAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTg                            <  1:651921/143‑1 (MQ=32)
                                                                                                                   aCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCa                          <  1:726944/143‑1 (MQ=33)
                                                                                                                        tGCGTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGACCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                     <  2:236530/143‑1 (MQ=12)
                                                                                                                        tGCCTGATGCGCCGCTGACGTGACTTATCAGGCCCACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                     <  1:302638/143‑1 (MQ=25)
                                                                                                                        tGCCTGATGCGCCGACACAGGGACTTAGCGGGCCTACTGTGTGCACTGCACGGTAGGTCGGATAAGACGCATAGCGGCGCCTCCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                     <  2:648002/143‑1 (MQ=11)
                                                                                                                         gCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGCGCACTGAAc                    <  1:924517/143‑1 (MQ=25)
                                                                                                                           cTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGACCTg                  >  1:1055108/1‑143 (MQ=25)
                                                                                                                                gcgcCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                                                                 >  2:3545/1‑59 (MQ=38)
                                                                                                                                  tccgcTGGCGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg           <  2:484411/142‑1 (MQ=17)
                                                                                                                                    cgcttgcGTGACTTAGCAGGCCTACGGGGTGCACTGAACTGTAGGTGGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa         <  1:142554/137‑1 (MQ=255)
                                                                                                                                           gTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACg  >  2:809666/1‑143 (MQ=33)
                                                                                                                                           gTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAAAACCGATGCCTGATGAACCGCTGAAGCGACTTATCAGGCATAGGGGGAGCACTCAACAGTAGGACGGATAAGACg  >  2:671966/1‑143 (MQ=11)
                                                                                                                                            |                                                                                                                                             
CGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACG  >  NC_000913/349577‑349858

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: