Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,739,103 T→C 7.7% S7G (AGC→GGC)  sdhC ← succinate dehydrogenase membrane b‑556 subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,739,1030TC 7.7% 141.8 / 2.8 52S7G (AGC→GGC) sdhCsuccinate dehydrogenase membrane b‑556 subunit
Reads supporting (aligned to +/- strand):  ref base T (23/25);  new base C (2/2);  total (25/27)
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.15e-01

TCGGAACCCAGCGACTTGCCCAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGCCGAGTATA  >  NC_002947/4738958‑4739250
                                                                                                                                                 |                                                                                                                                                   
tCGGAACCCAGCGACTTGCCCAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGGAGTTGGCtttt                                                                                                                                                  >  8:83753/1‑149 (MQ=255)
   gAACCCAGCGACTTGCCCAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCAc                                                                                                                                               <  3:307838/149‑1 (MQ=255)
         aGCGACTTGCCCAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCCCGGCttt                                                                                                                                         >  5:95948/1‑149 (MQ=255)
          gCGACTTGCCCAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCtttt                                                                                                                                        >  1:348145/1‑149 (MQ=255)
             aCTTGCCCAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCtttttt                                                                                                                                      >  4:297387/1‑148 (MQ=255)
              cTTGCCCAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTca                                                                                                                                    <  5:248915/149‑1 (MQ=255)
               ttGCCCAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTca                                                                                                                                    >  6:35596/1‑148 (MQ=255)
                  cccAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACt                                                                                                                                <  3:222552/149‑1 (MQ=255)
                          aTAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCttt                                                                                                                                                   <  2:226347/122‑1 (MQ=255)
                          aTAAAGCATGATGGCAAGGCCGATGAACAGAATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCttt                                                                                                                                                   >  1:226347/1‑122 (MQ=255)
                             aaGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGcccc                                                                                                                      >  2:156999/1‑148 (MQ=255)
                              aGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGcccc                                                                                                                      >  1:332996/1‑147 (MQ=255)
                              aGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGcccc                                                                                                                      >  1:259324/1‑147 (MQ=255)
                                    atgGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAAt                                                                                                              >  5:202034/1‑149 (MQ=255)
                                     tgGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAAt                                                                                                              <  6:252961/148‑1 (MQ=255)
                                       gCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAg                                                                                                           <  5:123871/149‑1 (MQ=255)
                                       gCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGATTACAGGGCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAg                                                                                                           <  2:88348/149‑1 (MQ=255)
                                        cAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAgg                                                                                                          <  5:378555/149‑1 (MQ=255)
                                        cAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAgg                                                                                                          >  1:317502/1‑149 (MQ=255)
                                            gCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCa                                                                                                       >  2:348082/1‑148 (MQ=255)
                                                ttGAACAGGATGTCGCAGTAGATACTGTGCATGATGGACTTGGACGCAGTGACGGGTAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTAGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAgtt                                                                                                  <  7:155416/148‑1 (MQ=255)
                                                   aaCAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAgttgtt                                                                                               <  5:35596/149‑1 (MQ=255)
                                                    aCAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGttgttg                                                                                              <  7:128778/149‑1 (MQ=255)
                                                      aGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGt                                                                                            <  1:11862/149‑1 (MQ=255)
                                                        gATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTaa                                                                                          >  7:245427/1‑149 (MQ=255)
                                                        gATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTaa                                                                                          <  7:274481/149‑1 (MQ=255)
                                                        gATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTaa                                                                                          >  8:274481/1‑149 (MQ=255)
                                                        gATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTaa                                                                                          <  4:311296/149‑1 (MQ=255)
                                                              gCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCCTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTAc                                                                                    >  8:48920/1‑149 (MQ=255)
                                                              gCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTACGTTTACAGGTCGATGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTGAGTGGTAAGTGTAc                                                                                    >  1:132450/1‑149 (MQ=255)
                                                                cGGAGATACGGTGCACGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTAAGGACTAGGTTTACAGGTCGTTGGCTTTTCACGGCGTTTTTCACACAGAACAGACCCTAGCGAAACAGGGCAAAGTTCTTGGTAAGTGTACTg                                                                                  >  6:251956/1‑149 (MQ=255)
                                                                    gaTACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCa                                                                              <  6:95948/149‑1 (MQ=255)
                                                                     aTACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAg                                                                             <  3:334093/149‑1 (MQ=255)
                                                                        cGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTa                                                                          <  6:202034/149‑1 (MQ=255)
                                                                             cAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACcca                                                                      <  3:297387/148‑1 (MQ=255)
                                                                                    gACGTGTACGCAGTGACGGGGAGCTTGATTGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAg                                                              <  5:284052/149‑1 (MQ=255)
                                                                                       gtgtACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAggg                                                            >  3:145924/1‑148 (MQ=255)
                                                                                              cAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCCTTTCACGGCtttt                                                                                                                                        <  6:369466/65‑1 (MQ=255)
                                                                                              cAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCCTTTCACGGCtttt                                                                                                                                        >  5:369466/1‑65 (MQ=255)
                                                                                                       ggAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAg                                                              >  4:257756/1‑130 (MQ=255)
                                                                                                       ggAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAg                                                              <  3:257756/130‑1 (MQ=255)
                                                                                                                   tCCTTAGGTCTAGGTTTACAGGGCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGGCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGc                               <  7:162427/149‑1 (MQ=255)
                                                                                                                           tCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCt                       >  7:46905/1‑149 (MQ=255)
                                                                                                                           tCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCACTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTCGTCAGGTACACACCACACAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCt                       >  3:310548/1‑149 (MQ=255)
                                                                                                                             tAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCAAACACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTgg                     >  7:219820/1‑149 (MQ=255)
                                                                                                                               ggTTTACACGTCGTTGGCTTTGCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTCAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGGAAGCCCCTgg                     <  2:132450/147‑1 (MQ=255)
                                                                                                                                 tttACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTTAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGt                 >  2:334792/1‑149 (MQ=255)
                                                                                                                                 tttACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGt                 >  7:95755/1‑149 (MQ=255)
                                                                                                                                    aCAGGTCGTTGGCCTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCg               <  7:48920/148‑1 (MQ=255)
                                                                                                                                         tCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGcc         <  2:332996/149‑1 (MQ=255)
                                                                                                                                           gTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGCCg        <  1:348082/148‑1 (MQ=255)
                                                                                                                                            ttGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGCCGAg      >  2:372594/1‑149 (MQ=255)
                                                                                                                                              ggCTTTTCACGGCTTTTTTCACACCGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGCCGAGt     <  5:251956/148‑1 (MQ=255)
                                                                                                                                                cTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTTAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGCCGAGtata  <  1:334792/149‑1 (MQ=255)
                                                                                                                                                cTTTTCACGACTTTTTTCACACTGAAGAGCACCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGCCGAGtata  >  8:277147/1‑149 (MQ=255)
                                                                                                                                                 ttttCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGCCGAGtata  <  2:259324/148‑1 (MQ=255)
                                                                                                                                                 ttttCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGCCGAGta    >  4:56707/1‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCGGAACCCAGCGACTTGCCCAGTGCATAAAGCATGATGGCAAGGCCGATGAACAGGATGACGCCGGAGATACGGTGCAGGATGGACGTGTACGCAGTGACGGGGAGCTTGATGGTCCTTAGGTCTAGGTTTACAGGTCGTTGGCTTTTCACGGCTTTTTTCACACTGAAGAGCCCCTAGCGAATCAGGGCAAAGTTGTTGGTAAGTGTACTGGTCAGGTACCCACCACCCAGGGAACGGCGACCCCCAGCAGTTCGGGCTTGCAAGCCCCTGGCGGTCGGGTGCCGAGTATA  >  NC_002947/4738958‑4739250

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