Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,664,179 T→C 7.8% intergenic (+161/‑2) prpC → / → acnA‑II methylcitrate synthase/aconitate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,664,1790TC 7.8% 145.3 / 2.2 52intergenic (+161/‑2)prpC/acnA‑IImethylcitrate synthase/aconitate hydratase
Reads supporting (aligned to +/- strand):  ref base T (26/21);  new base C (2/2);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GAGGCAACCCGTTTCCCTTGAGGGAGCGGGCGTGCCCGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGCCATACACGTCCCGCGTGCTCGCCGAGAACCTGGTGCG  >  NC_002947/2664030‑2664320
                                                                                                                                                     |                                                                                                                                             
gagGCAACCCGTTTCCCTTGAGGGAGCGGGCGTGCCCGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCat                                                                                                                                               <  2:512806/150‑1 (MQ=255)
  ggCAACCCGTTTCCCTTGAGGGAGCGGGCGTGCCCGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCatga                                                                                                                                             <  1:1345425/150‑1 (MQ=255)
  ggCAACCCGTTTCCCTTGAGGGAGCGGGCGTGCCCGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGACTGATACCatga                                                                                                                                             <  2:708266/150‑1 (MQ=255)
   gCAACCCGTTTCCCTTGAGGGAGCGGGCTTGCCCGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGGGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCaggat                                                                                                                                            >  2:560787/1‑150 (MQ=255)
                 ttGAGGGAGCGGGCGTGCCCGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATAcc                                                                                                                              <  1:220480/150‑1 (MQ=255)
                     gggAGCGGGCGTGCCCGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCaa                                                                                                                          <  2:1054954/150‑1 (MQ=255)
                                  cccGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGc                                                                                                             >  1:893580/1‑150 (MQ=255)
                                   ccGCGAATGCACCCGGGCTTGTAACCACTGTGCGCCTGCTGACGGATTGGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCa                                                                                                            <  2:1424003/150‑1 (MQ=255)
                                   ccGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCa                                                                                                            >  1:1345314/1‑150 (MQ=255)
                                   ccGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCACGATGAACACTGCATACCGCAAGCACCTGCCAGGCa                                                                                                            >  2:1161775/1‑150 (MQ=255)
                                   ccGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATAACGGGAGGAACACTGCATACCGCAAGCACCTGCCAGGcc                                                                                                            >  2:690609/1‑149 (MQ=255)
                                    cgcgAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCAc                                                                                                           >  1:163729/1‑150 (MQ=255)
                                       gAATGCGCCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGa                                                                                                        <  1:1133182/150‑1 (MQ=255)
                                            cACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTgg                                                                                                   <  1:534930/150‑1 (MQ=255)
                                             aCCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGa                                                                                                  >  1:735737/1‑150 (MQ=255)
                                              ccgggCCTTGTAACCACTGTGCACCTGCTGACGCATCAGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTGCGACTCCTGTGACCGAGCCTGATGCGATGATGAACACTGCATACCGCAAGCACCTCCCCGGCACCGACCTggcc                                                                                                 >  2:65715/4‑148 (MQ=255)
                                               ccGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCCGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGcc                                                                                                                 <  2:826575/133‑1 (MQ=255)
                                               ccGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCCGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGcc                                                                                                                 >  1:826575/1‑133 (MQ=255)
                                                   ccTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACtt                                                                                            <  2:144542/150‑1 (MQ=255)
                                                   ccTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACtt                                                                                            <  2:893580/150‑1 (MQ=255)
                                                     ttGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCg                                                                                          >  1:1111928/1‑150 (MQ=255)
                                                                   cagctgctgACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACAGCAAACACCTGCAAGGCACCGACCTGGACTACTTCGACGCCCGCgcggcg                                                                            >  1:1119217/4‑150 (MQ=255)
                                                                       tgctgACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTggca                                                                                                 >  1:238490/1‑123 (MQ=255)
                                                                       tgctgACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTggca                                                                                                 <  2:238490/125‑3 (MQ=255)
                                                                            aCGCATTCGCGGGCACGCCCGCTCCAACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCg                                                                   >  1:634727/1‑150 (MQ=255)
                                                                              gCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGAt                                                                 >  2:923298/1‑150 (MQ=255)
                                                                                aTTCGCGGGCACGCCCGCTCCCACAGGAACGGTGAAACCTTGGGAAGCCTGTGACCGAGCCTGTTTCCATCATGAACACTTCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCa                                                               >  2:162859/1‑150 (MQ=255)
                                                                                 ttCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTg                                                                                                    >  1:184887/1‑112 (MQ=255)
                                                                                 ttCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTg                                                                                                    <  2:184887/112‑1 (MQ=255)
                                                                                     cgGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGccc                                                          <  1:253140/150‑1 (MQ=255)
                                                                                            gCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCt                                                   <  1:923298/150‑1 (MQ=255)
                                                                                             cccGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCACGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGa                                                                  <  1:842893/134‑1 (MQ=255)
                                                                                             cccGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCACGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGa                                                                  >  2:842893/1‑134 (MQ=255)
                                                                                              ccGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTac                                                 >  2:343251/1‑150 (MQ=255)
                                                                                                 cTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTacgac                                              <  1:417520/150‑1 (MQ=255)
                                                                                                   cccACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACgg                                            <  2:1345314/150‑1 (MQ=255)
                                                                                                    ccACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCg                                                         <  1:1294984/136‑1 (MQ=255)
                                                                                                    ccACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCg                                                         >  2:1294984/1‑136 (MQ=255)
                                                                                                      acaGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCaa                                                                                                                          >  2:977153/1‑69 (MQ=255)
                                                                                                      acaGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCaa                                                                                                                          <  1:977153/69‑1 (MQ=255)
                                                                                                         ggTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCACGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGcc                                      <  1:1161775/150‑1 (MQ=255)
                                                                                                          gTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGAc                                                                                                       <  2:880169/84‑1 (MQ=255)
                                                                                                          gTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGAc                                                                                                       >  1:880169/1‑84 (MQ=255)
                                                                                                           tACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGACGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGCCAt                                    >  1:435705/1‑150 (MQ=255)
                                                                                                               gTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGCCATacac                                >  1:7883/1‑150 (MQ=255)
                                                                                                                   aaCCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATCCCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGCCATACACGTcc                            >  1:1350461/1‑150 (MQ=255)
                                                                                                                       ttCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGcc                                                    <  2:1151893/122‑1 (MQ=255)
                                                                                                                       ttCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGCCATACACGTCCcgcg                        >  1:491555/1‑150 (MQ=255)
                                                                                                                       ttCCGAATCCTGTGACCGAGCCTGATTCCATGATGAAAACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGcc                                                    >  1:1151893/1‑122 (MQ=255)
                                                                                                                               ccTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGCCATACACGTCCCGCGTGCTCGcc                <  1:343251/150‑1 (MQ=255)
                                                                                                                                           ccTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGAc                                                                                        >  1:887248/1‑66 (MQ=255)
                                                                                                                                           ccTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGAc                                                                                        <  2:887248/66‑1 (MQ=255)
                                                                                                                                            cTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGCCATACACGTCCCGCGTGCTCGCCGAGAACCTGGTgc   >  1:495138/1‑150 (MQ=255)
                                                                                                                                             tGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGCCATACACGTCCCGCGTGCTCGCCGAGAACCTGGTgcg  >  1:408893/1‑150 (MQ=255)
                                                                                                                                                   cATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGAt                                                                 >  1:1272888/1‑81 (MQ=255)
                                                                                                                                                   cATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGAt                                                                 <  2:1272888/81‑1 (MQ=255)
                                                                                                                                                    atgatgAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAgg                                                                     >  1:1326829/1‑76 (MQ=255)
                                                                                                                                                    atgatgAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAgg                                                                     <  2:1326829/76‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GAGGCAACCCGTTTCCCTTGAGGGAGCGGGCGTGCCCGCGAATGCACCCGGCCTTGTAACCACTGTGCACCTGCTGACGCATTCGCGGGCACGCCCGCTCCCACAGGTACGGTGAAACCTTCCGAATCCTGTGACCGAGCCTGATTCCATGATGAACACTGCATACCGCAAGCACCTGCCAGGCACCGACCTGGACTACTTCGACGCCCGCGCGGCGGTCGAGGCGATCAAGCCCGGCGCCTACGACGGCTTGCCATACACGTCCCGCGTGCTCGCCGAGAACCTGGTGCG  >  NC_002947/2664030‑2664320

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