Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,533,413 A→G 6.8% E16E (GAA→GAG atoA → 3‑oxoacid CoA‑transferase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,533,4130AG 6.8% 175.0 / 2.6 59E16E (GAA→GAGatoA3‑oxoacid CoA‑transferase subunit A
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base G (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

ATGACGAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGTTTTGGTCT  >  NC_002947/3533264‑3533556
                                                                                                                                                     |                                                                                                                                               
aTGACGAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGaa                                                                                                                                                 >  1:799624/1‑150 (MQ=255)
aTGACGAAAAAAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGaa                                                                                                                                                 >  2:310238/1‑150 (MQ=255)
         aCAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAGGGCCTGAcc                                                                                                                                        >  2:75692/1‑150 (MQ=255)
              aTTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGacaa                                                                                                                                   <  1:644415/150‑1 (MQ=255)
              aTTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGacaa                                                                                                                                   <  2:1329299/150‑1 (MQ=255)
                    aaaTCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGAc                                                                                                                             <  2:1081623/150‑1 (MQ=255)
                           tCGCTTGCGGTGAAGATGTCTGCACAGGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTg                                                                                                                      >  1:630529/1‑150 (MQ=255)
                           tCGCTTGCGGTGAAGATGTCTGCACACGATCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTg                                                                                                                      <  2:1039170/150‑1 (MQ=255)
                                tGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAGGGCCTGACCGACAACATGACGGTACTGGCCgg                                                                                                                 <  1:75692/150‑1 (MQ=255)
                                  cGGTGAAGATGTCTGCACCCGCTCACCGACAGCCTCGTCGGTTAAACACCTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGTTCTAACCGACAACATGACGGTACTGGCCcgtg                                                                                                               <  2:474653/150‑1 (MQ=255)
                                          aTGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAg                                                                                                                                                <  2:258834/109‑1 (MQ=255)
                                          aTGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAg                                                                                                                                                >  1:258834/1‑109 (MQ=255)
                                           tGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGcc                                                                                                      <  2:1304718/150‑1 (MQ=255)
                                           tGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGcc                                                                                                      >  2:716975/1‑150 (MQ=255)
                                                 cacaCGCTCACCGACAGCCTCGTCGGCTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGCCGGTACTGGCCGGTGGTTTCCGCCTGAGCg                                                                                                >  2:157199/1‑150 (MQ=255)
                                                            cGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACg                                                                                                                            <  2:1150519/111‑1 (MQ=255)
                                                            cGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACg                                                                                                                            >  1:1150519/1‑111 (MQ=255)
                                                               cAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAAc                                                                                  <  1:716975/150‑1 (MQ=255)
                                                                  ccTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCt                                                                                                                                            <  2:1361120/89‑1 (MQ=255)
                                                                  ccTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCt                                                                                                                                            >  1:1361120/1‑89 (MQ=255)
                                                                      gtcgGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCa                                                                           <  1:310238/150‑1 (MQ=255)
                                                                          gTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGa                                                                                                                              <  1:349988/95‑1 (MQ=255)
                                                                          gTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGa                                                                                                                              >  2:349988/1‑95 (MQ=255)
                                                                               acaacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGt                                                                                                                          <  2:964037/94‑1 (MQ=255)
                                                                               acaacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGt                                                                                                                          >  1:964037/1‑94 (MQ=255)
                                                                               acaacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCa                                                                  >  2:742251/1‑150 (MQ=255)
                                                                               acaacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCAGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCa                                                                  >  1:1441620/1‑150 (MQ=255)
                                                                                caacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAg                                                                                                                                                >  2:437340/1‑71 (MQ=255)
                                                                                caacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAg                                                                                                                                                <  1:437340/71‑1 (MQ=255)
                                                                                 aacTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAg                                                                >  1:808799/1‑150 (MQ=255)
                                                                                  acTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAgc                                                               <  2:808799/150‑1 (MQ=255)
                                                                                    tACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGa                                                                                                                                          <  1:1530684/73‑1 (MQ=255)
                                                                                    tACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGa                                                                                                                                          >  2:1530684/1‑73 (MQ=255)
                                                                                    tACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATccc                                                                                         <  1:212487/122‑1 (MQ=255)
                                                                                    tACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATccc                                                                                         >  2:212487/1‑122 (MQ=255)
                                                                                       aaaaGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAgcggcgc                                                          >  1:66267/1‑150 (MQ=255)
                                                                                                     cATGGCCGGACTGGACAAGCGCGTTGCNACCTATGAGCAGNCCCTCGNAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCt                                            >  1:1296756/1‑150 (MQ=255)
                                                                                                         gCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGc                                                                                                 <  1:495971/93‑1 (MQ=255)
                                                                                                         gCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGc                                                                                                 >  2:495971/1‑93 (MQ=255)
                                                                                                         gCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCATCGAAATCAAGCGGCGCGGCGTCaa                                                  <  1:851325/140‑1 (MQ=255)
                                                                                                         gCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCATCGAAATCAAGCGGCGCGGCGTCaa                                                  >  2:851325/1‑140 (MQ=255)
                                                                                                         gCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGAcc                                        <  1:157199/150‑1 (MQ=255)
                                                                                                            ggACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTg                                     <  1:742251/150‑1 (MQ=255)
                                                                                                                  gACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTc                                >  2:963345/1‑149 (MQ=255)
                                                                                                                  gACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTc                                <  1:963345/149‑1 (MQ=255)
                                                                                                                       ccgcgTTGCAACATATGACCAGACACTCGAAGGCCTGACCGACAACATGACGGTACGGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGACCGGACCGTGGTTTCcaacaa                          >  2:1504430/2‑150 (MQ=255)
                                                                                                                              gCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTgcggc                   <  2:66267/150‑1 (MQ=255)
                                                                                                                              gCAACCTAGGAGCAGGCCCTCGAAGGCCTGACCGACAACATGCAAGGACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTgcggc                   <  1:589739/150‑1 (MQ=255)
                                                                                                                                aaCCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGt                 >  1:856763/1‑150 (MQ=255)
                                                                                                                                aaCCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGt                 >  2:825513/1‑150 (MQ=255)
                                                                                                                                    gtcggtcAGGCCTTCGAGGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGAc             >  2:838197/7‑150 (MQ=255)
                                                                                                                                    gtcggtcAGGCCTTCGAGGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGAc             <  1:838197/144‑1 (MQ=255)
                                                                                                                                      tGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACgg           <  1:825513/150‑1 (MQ=255)
                                                                                                                                        aGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGtt         >  1:990890/1‑150 (MQ=255)
                                                                                                                                          cAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGtttt       <  2:630529/150‑1 (MQ=255)
                                                                                                                                           aGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGTTTTg      <  1:831414/150‑1 (MQ=255)
                                                                                                                                              cccTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGTTTTGGTc   <  2:856763/150‑1 (MQ=255)
                                                                                                                                               ccTCGAAGGCCTGACCGACAACATGACNGTACTGGCCGGTNGTTTCGNCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGTTTTGGTCt  >  1:1329178/1‑150 (MQ=255)
                                                                                                                                               ccTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGTTTTGGTCt  <  1:68478/150‑1 (MQ=255)
                                                                                                                                                 tCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTgcggc                   <  1:1523437/131‑1 (MQ=255)
                                                                                                                                                 tCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTgcggc                   >  2:1523437/1‑131 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ATGACGAAAACAAAATTAACAAATCAGTCGCTTGCGGTGAAGATGTCTGCACACGCTCACCGACAGCCTCGTCGGTTAAACAACTACAAAAGGCGGATCAGCATGGCCGGACTGGACAAGCGCGTTGCAACCTATGAGCAGGCCCTCGAAGGCCTGACCGACAACATGACGGTACTGGCCGGTGGTTTCGGCCTGTGCGGCATCCCGGAAAACCTTATCAGCGAAATCAAGCGGCGCGGCGTCAAGGGCCTGACCGTGGTTTCCAACAACTGCGGCGTGGACGGTTTTGGTCT  >  NC_002947/3533264‑3533556

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