Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,326,658 G→T 8.6% V471L (GTG→TTG)  mqo‑III → malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,326,6580GT 8.6% 127.7 / 2.1 43V471L (GTG→TTG) mqo‑IIImalate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (25/14);  new base T (2/2);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTCGGTACCGAGGTGGTGACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGC  >  NC_002947/3326512‑3326802
                                                                                                                                                  |                                                                                                                                                
tcatcaAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTCGGTACCGAGGTGGTGACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGgtg                                                                                                                                                <  3:9119/149‑1 (MQ=255)
  atcaAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTCGGTACCGAGGTGGTGACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGtt                                                                                                                                              >  3:255121/1‑149 (MQ=255)
     aagaatgaCGAAGAGTTGGGCGGCGTGCTCAAGCTCGGTACCGAGGTGGTGACGGCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCaa                                                                                                                                           <  8:193169/149‑1 (MQ=255)
     aagaaagaCGAAGCGTTGGGCGGCGTGCTCAAGCTCGGTACCGATGTGGTGACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGTGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCaa                                                                                                                                           <  8:198216/149‑1 (MQ=255)
           gaCGAAGCGTTGGGCGGCGTGCTCAAGCTCGGTACCGAGGTGGTGACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAAc                                                                                                                                       >  1:47238/1‑147 (MQ=255)
                               gCTCAAGCTCGGTACCGAGGTGGTGACGTCCAGGGACGGCAGTATCGCCCGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTgg                                                                                                                  >  1:45985/1‑148 (MQ=255)
                                              cGAGGTGGTGACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCGAAGATCaaa                                                                                                   >  1:263656/1‑148 (MQ=255)
                                               gagaTGGTGACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCaaa                                                                                                   >  2:229810/1‑147 (MQ=255)
                                                    ggtgACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGatca                                                                                            >  2:267456/1‑149 (MQ=255)
                                                      tgACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTgg                                                                                                                  >  1:156624/1‑125 (MQ=255)
                                                      tgACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTgg                                                                                                                  <  2:156624/125‑1 (MQ=255)
                                                          gTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGt                                                                                      >  1:128712/1‑149 (MQ=255)
                                                                     gCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATc                                                                           <  4:255121/149‑1 (MQ=255)
                                                                           tCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAACCAGATCATCCCGTCCTATGGCATCCGCCTg                                                                     >  1:87682/1‑149 (MQ=255)
                                                                             gCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGaa                                                                   >  8:46904/1‑149 (MQ=255)
                                                                                  cTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGacc                                                              <  2:47238/149‑1 (MQ=255)
                                                                                       tGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATg                                                                                <  6:53622/126‑1 (MQ=255)
                                                                                       tGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATg                                                                                >  5:53622/1‑126 (MQ=255)
                                                                                       tGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATACCGTCCTAAGGCATCCGCCTGAACGACCACCCg                                                         >  6:18949/1‑149 (MQ=255)
                                                                                        gggggCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCg                                                         >  2:19349/1‑148 (MQ=255)
                                                                                           ggCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCTGTACTATGGCATCAGCATGAACGACCACCAGGAt                                                      >  8:134394/1‑148 (MQ=38)
                                                                                               tctcCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGc                                                  >  3:142627/1‑148 (MQ=255)
                                                                                                             gACCGCTGAGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTggg                                   <  2:45985/149‑1 (MQ=255)
                                                                                                                cGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCa                                 >  5:22830/1‑148 (MQ=255)
                                                                                                                  cTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCTTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCatat                              <  5:18949/149‑1 (MQ=255)
                                                                                                                    gcgcCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATAc                            >  6:93513/1‑149 (MQ=255)
                                                                                                                      gcCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATAcc                           >  5:320671/1‑148 (MQ=255)
                                                                                                                       ccTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCg                                                         >  6:255331/1‑117 (MQ=255)
                                                                                                                       ccTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCg                                                         <  5:255331/117‑1 (MQ=255)
                                                                                                                             aTGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAgg                    <  1:267456/148‑1 (MQ=255)
                                                                                                                                   gATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGtgct                <  5:323469/146‑1 (MQ=255)
                                                                                                                                   gATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGtgct                >  6:323469/1‑146 (MQ=255)
                                                                                                                                    tgttcttGAACACCTTGTTTAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGt                                      >  4:173814/6‑123 (MQ=255)
                                                                                                                                    tgttcttGAACACCTTGTTTAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGt                                      <  3:173814/118‑1 (MQ=255)
                                                                                                                                    aTGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAg            >  6:258114/1‑149 (MQ=255)
                                                                                                                                      ttcttGAACACCTTGTTCAAGAACAAGGTGGGGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATAcc                           >  6:180388/4‑132 (MQ=255)
                                                                                                                                      ttcttGAACACCTTGTTCAAGAACAAGGTGGGGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATAcc                           <  5:180388/129‑1 (MQ=255)
                                                                                                                                        gtTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCtt         <  5:93513/148‑1 (MQ=255)
                                                                                                                                        gtTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCtt         >  5:113114/1‑148 (MQ=255)
                                                                                                                                        gtTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCtt         >  4:13885/1‑148 (MQ=255)
                                                                                                                                        gtTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGGGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGc                                                  <  2:59641/107‑1 (MQ=255)
                                                                                                                                         ttGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCtt         >  6:187962/1‑147 (MQ=255)
                                                                                                                                           gAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCGAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAgc     >  8:251245/1‑149 (MQ=255)
                                                                                                                                              cAAGGTGGTCAAGAACAAGGTGGCGGCGCCCGAGGGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAgccgc  <  2:263656/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTCGGTACCGAGGTGGTGACGTCCAGGGACGGCAGTATCGCCGGCTTGCTGGGGGCCTCTCCCGGTGCGTCGACCGCTGCGCCTATCATGCTCGATGTGTTGAACAAGGTGTTCAAGAACAAGGTGGCGTCGCCCGAGTGGCAGGCCAAGATCAAACAGATCATCCCGTCCTATGGCATCCGCCTGAACGACCACCCGGATCAGCTAGAGAAGGAGTGGGCATATACCAACGAGGTGCTGCAGCTTGAGCCGC  >  NC_002947/3326512‑3326802

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: ‑