Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,326,411 G→A 8.5% V388V (GTG→GTA)  mqo‑III → malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,326,4110GA 8.5% 129.6 / 2.7 47V388V (GTG→GTA) mqo‑IIImalate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base A (2/2);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

TGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTCG  >  NC_002947/3326264‑3326553
                                                                                                                                                   |                                                                                                                                               
tGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGtg                                                                                                                                                 <  8:336462/148‑1 (MQ=255)
    ccGTTTGCCACCTTCTCGACAAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGCACCTGCCCGATCAGGTACTGCa                                                                                                                                            <  4:21678/149‑1 (MQ=255)
           ccACCTTCTCGACCAAGGTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGGGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGAc                                                                                                                                     <  7:204289/149‑1 (MQ=255)
                  ctcGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGAcc                                                                                                                              >  8:336019/1‑149 (MQ=255)
                     gACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCt                                                                                                                           >  5:29906/1‑149 (MQ=255)
                         aaGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATc                                                                                                                                         <  8:390830/131‑1 (MQ=255)
                         aaGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATc                                                                                                                                         >  7:390830/1‑131 (MQ=255)
                              ccTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGAt                                                                                                                                  >  5:80365/1‑133 (MQ=255)
                              ccTCAAGCAGGGCTCATTGCTTGACCTGCTCACCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGAt                                                                                                                                  <  6:80365/133‑1 (MQ=255)
                                    gCAGGGCTCATTGCTTGACCTGCTCGCCAGAACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCa                                                                                                            >  1:33776/1‑149 (MQ=255)
                                        ggCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTa                                                                                                        <  4:285042/149‑1 (MQ=255)
                                           tCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACtt                                                                                                     >  4:221076/1‑149 (MQ=255)
                                                gCTTGACCTGCTCGCCAGTACCAGCGTGGACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCTGTACCTGATCGGGCAGGtgctg                                                                                                                                              >  4:325355/1‑103 (MQ=255)
                                                gCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCTGTACCTGATCGGGCAGGtgctg                                                                                                                                              <  3:325355/103‑1 (MQ=255)
                                                  ttGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTgc                                                                                                               <  8:134106/132‑1 (MQ=255)
                                                  ttGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTgc                                                                                                               >  7:134106/1‑132 (MQ=255)
                                                          cTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCgaaga                                                                                       >  1:383547/1‑148 (MQ=255)
                                                             gCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGaagaaag                                                                                    >  3:210821/1‑148 (MQ=255)
                                                                   aCCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACt                                                                              <  3:221076/148‑1 (MQ=255)
                                                                     cAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTggc                                                                           >  8:123102/1‑149 (MQ=255)
                                                                       gCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTggcgg                                                                         <  2:33776/149‑1 (MQ=255)
                                                                       cagTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTggcgg                                                                         >  7:267489/3‑149 (MQ=255)
                                                                            cacaATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCtgt                                                                     <  5:415209/148‑1 (MQ=255)
                                                                             acaATTCCTGGGCGATGGTGCGTGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTTCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCtgtg                                                                    <  8:267489/148‑1 (MQ=255)
                                                                             acaATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTgg                                                                   >  8:122560/1‑149 (MQ=255)
                                                                             acaATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTgg                                                                   <  7:336019/149‑1 (MQ=255)
                                                                                   ccTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGc                                                              <  2:383547/148‑1 (MQ=255)
                                                                                       gCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCg‑gg                                                          <  4:210821/147‑1 (MQ=255)
                                                                                       gCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGc                                                         >  4:22316/1‑148 (MQ=255)
                                                                                         cGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAg                                                      >  8:225693/1‑149 (MQ=255)
                                                                                          gATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGc                                                     <  7:122560/149‑1 (MQ=255)
                                                                                                    gTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAat                                           <  3:22316/149‑1 (MQ=255)
                                                                                                           tACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAatcatc                                       >  3:175194/1‑146 (MQ=255)
                                                                                                            aCGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCaag                                    >  4:346069/1‑148 (MQ=255)
                                                                                                                 aGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGAc                              >  8:31715/1‑149 (MQ=255)
                                                                                                                   ttCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACg                             >  3:306417/1‑148 (MQ=255)
                                                                                                                              gcagcaCCTGCCCGATCAGGTACTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGggcg                 >  3:304093/1‑149 (MQ=25)
                                                                                                                              gcagcaCCTGCCCGATCAGGTACTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGggcg                 >  6:386365/1‑149 (MQ=25)
                                                                                                                              gCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTcc                                                                                                   >  7:313456/1‑68 (MQ=255)
                                                                                                                              gCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTcc                                                                                                   <  8:313456/68‑1 (MQ=255)
                                                                                                                                  tACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGggcggc               <  4:254464/147‑1 (MQ=255)
                                                                                                                                  tACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGggcggc               >  3:254464/1‑147 (MQ=255)
                                                                                                                                  tACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGt             >  4:79865/1‑149 (MQ=255)
                                                                                                                                     cTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCt          <  6:29906/149‑1 (MQ=255)
                                                                                                                                     cTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCt          >  3:387123/1‑149 (MQ=255)
                                                                                                                                         tCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAg      <  7:225693/149‑1 (MQ=255)
                                                                                                                                           gggCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGATGCAGGCGATATAAGCGGGTGAAAATCATAAAGAAAGACGAAGCGTTGGGCGGCGTGATAAAGc     >  1:337173/1‑149 (MQ=25)
                                                                                                                                           gggCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGAGCACATACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCt    >  1:154695/1‑149 (MQ=255)
                                                                                                                                            atcAGGTACTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTc   <  4:304093/147‑1 (MQ=255)
                                                                                                                                             gCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCAACTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTCg  <  2:154695/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCG‑GGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTCG  >  NC_002947/3326264‑3326553

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