Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,654,515 A→C 32.7% T594P (ACC→CCC)  nuoG → NADH‑quinone oxidoreductase subunit G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,654,5150AC32.7% 50.5 / 6.4 46T594P (ACC→CCC) nuoGNADH‑quinone oxidoreductase subunit G
Reads supporting (aligned to +/- strand):  ref base A (20/10);  new base C (2/14);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 5.80e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGCGGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGCTTCTTCCAGGTGTTCGACCCGCAGTACCTGGACA  >  ACB143/4654382‑4654644
                                                                                                                                     |                                                                                                                                 
aGACGCGGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATAGTTGCCGACCACTCCAACAcc                                                                                                                                 >  6:376678/1‑136 (MQ=255)
    gcgGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCa                                                                                                                             >  1:201630/1‑136 (MQ=255)
    gcgGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCa                                                                                                                             >  3:88085/1‑136 (MQ=255)
         catcaGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGCCCACTCCAAACCCGCCACCGTc                                                                                                                        >  7:188771/1‑136 (MQ=255)
         catcaGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGCCCACTCCAAAACCGCCACCGTc                                                                                                                        >  4:319415/1‑136 (MQ=255)
         catcaGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTc                                                                                                                        >  5:274785/1‑136 (MQ=255)
         catcaGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCAGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCgac                                                                                                                        >  6:219283/1‑134 (MQ=255)
          atcaGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTc                                                                                                                        >  4:165314/1‑135 (MQ=255)
              gCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACcg                                                                                                                   <  3:146081/136‑1 (MQ=255)
                     aGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCAc                                                                                                            >  1:48679/1‑136 (MQ=255)
                               aTCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGtgctgc                                                                                                  <  7:245360/136‑1 (MQ=255)
                                           gAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGt                                                                                      >  1:209100/1‑136 (MQ=255)
                                             aaaCGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAACCCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCGCGtac                                                                                    >  5:63544/1‑134 (MQ=255)
                                              aaCGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAACCCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCg                                                                                   <  3:289079/136‑1 (MQ=255)
                                              aaCGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCg                                                                                   <  7:3559/136‑1 (MQ=255)
                                                ggACCTGTACGCCCGCGTAGCGGCCGCCTTGGGGGAGGCGGCCCTGGGGGCGGCCAAGGTGGGGATGGTTGCGGCCCCCTCAAAAACCGCCACCGTCGACCTCGCCCTCCTGGGGGTGCCGGCCGCGTCGTTCGcc                                                                                 <  4:67010/135‑1 (MQ=255)
                                                  actggTACCCCCGCGTCCGGGCCGCAAAGGTGGTTGGTGCCGTGGGTGCGGCCATGGTGGTGATGTTTGCCGACCACTCAAAACCCGCCACCGTCTACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGccta                                                                               <  6:391974/132‑3 (MQ=255)
                                                  aCCTGTCCGCCCGCGAACCGGCCGCCATGGTCGATGCGGCCCGGGCTGGGGCCATGGTGGTGAGCGTTGCCGACCACTCCAAACCCGCCCCCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGa                                                                               <  8:303911/136‑1 (MQ=255)
                                                     tgaaCGCCCGGGAACCGGCCGCAAAGGGCGTTGCTCCCCTGGGTGCGGCCATGGTGGTTTTCGTGGCCTACCACTCCAAAACCGCCACGGTCGACCGCGCACACCTGGTGCGGCCGGCCGCCTCGTGCGCCGAAgg                                                                            <  8:10524/133‑1 (MQ=255)
                                                     ggtccGCCCGCGTACGGGCCGCCACGGGCTTTGCTGCCCTGGCGGCGGCCACGGGGGGTTTCGTTGCCCACCACTCCAACCCCGCCCCCGTCGACCGCGCCCACCTGGGGCTGCCGGCCGCCGCGTTCGCCGACgg                                                                            <  4:64982/132‑1 (MQ=255)
                                                       tACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTAGCCGAAGGCg                                                                          <  7:53750/136‑1 (MQ=255)
                                                        aCGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGCCCGCGCCCACCTGGTGCTGCCGGCCGCCGCGTTCGCCGAAGGCGa                                                                         >  1:131037/1‑136 (MQ=255)
                                                        aCGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGa                                                                         >  1:212078/1‑136 (MQ=255)
                                                          gCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACg                                                                       >  6:301214/1‑136 (MQ=255)
                                                            ccGCGTACCGGCCGCCACGGTCGATGCGGCCCTGGCTGCGGCCAAGGTGGGGATCTTGGCCCACCACCCCAAACCCGCCACCGTCCCCCGCGCCCACCTGGGGCGGCCGGCCGCCTCGTTCGCCGACGGCGACGGt                                                                     <  2:26625/136‑1 (MQ=255)
                                                               cgTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCGAGTTCGCCGAAGGCGACGGTAcc                                                                  >  4:96069/1‑136 (MQ=255)
                                                                 tACCGGCCGCCAAGGGCGATGCGGCCCTGGCGGCGGCCAAGGGGGTGATCGTTGCCCACCACTCCAAACCCGCCACCGGCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTGCGCCGAAGGCGACGGTACCCt                                                                <  5:298501/136‑1 (MQ=255)
                                                                    cGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGt                                                             >  4:201940/1‑136 (MQ=255)
                                                                          ccACGGGCGATGCTGCCCTGGCGGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAACCCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGACGGCGACGGTACCCTGGTCAGCCa                                                       <  2:71799/136‑1 (MQ=255)
                                                                             aGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCCGCCAgg                                                     >  5:394340/1‑135 (MQ=255)
                                                                               gccGAAGCGGCCCTGGCTGCGGCCAAGGGGGGTAGCGTGGCCGCCCACTCCAAAACCCCCCCCGGCGCCCGCGCCCACCTGGTGCTGCCGGCCGCCGCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAg                                                  <  7:79297/134‑1 (MQ=255)
                                                                               gacGATGCTGCCCTGGCTGCGCCCAAGGTGGTGAGCGTTGCCCACCCCGCAAACCCCGCCACCGGCGCCCGCGCCCACCGGGGGCGGCGGGCCGCCTCGTTCGCCGACGGCGACGGTACCCTGGTCAGCCAGGAAg                                                  <  5:221973/134‑1 (MQ=255)
                                                                                      ctccccTGGGTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGCCCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTg                                            <  1:18645/132‑1 (MQ=255)
                                                                                         cccTGGCTGCGGCCATGGTGGTGATCGTTGCCCACCACGCCAAACCCGCCCCCGTCGACCGCGCCCACCGGGTGCTGCCGGCCGCCGCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCa                                        <  6:274785/136‑1 (MQ=255)
                                                                                            tGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAgcg                                     >  2:49894/1‑136 (MQ=255)
                                                                                                 ggggCCAAGGGGGTGATCGTTGCCGACCACTCAAAACCCGCCCCCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGcttct                                <  5:394086/134‑1 (MQ=255)
                                                                                                 gCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGcttct                                >  4:132268/1‑136 (MQ=255)
                                                                                                 gCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGcttct                                >  8:86714/1‑136 (MQ=255)
                                                                                                  cGGCCATGGGGGTGAGCGTTGCCCACCACCCAAAAACCCCCCCCGGCGCCCCCGCCCACGTGGGGCTGCCGGCCGCGGCGTCCGCCCACGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGcttctt                               <  7:149579/136‑1 (MQ=255)
                                                                                                                agcGTTGCCGACCACTCCAAACCCGCCACCGTCCACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGCTTCTTCCAGGTGTTCGAcc                 <  8:252031/134‑1 (MQ=255)
                                                                                                                agcGTTGCCGACCACTCCAAAACCGCCACCGTCGCCCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTGCGCCGACGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGCTTCTTCCAGGTGTTCGAcc                 <  5:301214/134‑1 (MQ=255)
                                                                                                                agcGTTGCCGACCACGCAAACCCCGCCCCCGGGCCCCGCGCCCACCGGGTGCTGCCGGCCGCCGGGTGCGCCGACGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGCTTCTTCCAGGTGTTCGAcc                 <  4:88085/134‑1 (MQ=255)
                                                                                                                      gCCGACCACTCCAAACCCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGCTTCTTCCAGGTGTTCGACCCGCAg            <  7:201191/135‑1 (MQ=255)
                                                                                                                       ccGACCACTCCAAAACCGCCACCGTCGACCGCGCCCCCCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCCGCGCTTCTTCCAGGTGTTCGACCCGCAGTa          >  8:146307/1‑136 (MQ=255)
                                                                                                                          cccACTCCAAACCCGCACCCGTCGACCGCGCCCACCTGGGGCGGCCGGCCGCCTCGTGCCCCGACGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGCTTCTTCCAGGTGTTCGACCCGCAGTACCt       <  3:371863/135‑1 (MQ=255)
                                                                                                                             aCTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGCTTCTTCCAGGTGTTCGACCCGCAGTACCTGGa    >  1:309950/1‑136 (MQ=255)
                                                                                                                               ccaaaacCCGCCCCCGGCGACCGCGCCCACCTGGTGGGGCCGGCCGCCTCGTGCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGCTTCTTCCAGGTGTTCGACCCGCAGTACCTGGAca  <  8:373245/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                 
AGACGCGGTCATCAGCGGCAAAGCCGATGCCATCGTGGTGCTGGAAAACGACCTGTACGCCCGCGTACCGGCCGCCAAGGTCGATGCTGCCCTGGCTGCGGCCAAGGTGGTGATCGTTGCCGACCACTCCAAAACCGCCACCGTCGACCGCGCCCACCTGGTGCTGCCGGCCGCCTCGTTCGCCGAAGGCGACGGTACCCTGGTCAGCCAGGAAGGCCGTGCCCAGCGCTTCTTCCAGGTGTTCGACCCGCAGTACCTGGACA  >  ACB143/4654382‑4654644

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