Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,497,924 C→T 13.3% G191G (GGC→GGT petA → Ubiquinol‑cytochrome c reductase iron‑sulfursubunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,497,9240CT13.3% 60.6 / 5.1 31G191G (GGC→GGTpetAUbiquinol‑cytochrome c reductase iron‑sulfursubunit
Reads supporting (aligned to +/- strand):  ref base C (12/14);  new base T (2/2);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

GCCCGTGCCACGGCTCTCACTATGACCTGGCTGGTCGTGTCTACAAGTCCCAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTATTACGCACCCAAGAACTTCAACTTCCTGTACTTCTT  >  ACB122/1497806‑1498054
                                                                                                                      |                                                                                                                                  
gCCCGTGCCACGGCTCTCACTATGACCTGGCTGGTCGTGTCTACAAGTCCCAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACg                                                                                                                    >  1:5705/1‑135 (MQ=255)
 cccTTGCAACGGCTCTCACTATGACCTGGCTGGTCGTGTCTACAATTCCGAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCAATCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCATGAGAACGc                                                                                                                   <  1:102905/135‑1 (MQ=255)
 cccGTGCCACGGCTCTCACTATGACCTGGCTGGTCGTGTCTACAAGTCCCAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGc                                                                                                                   <  5:203037/135‑1 (MQ=255)
       ccACGGCTCTCACTATGACCTGGCTGGTCGTGTCTACAAGTCCCAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGCCGAGTTCATCGTGATTGGCGGGGACCAGGAGACCGCatgat                                                                                                              >  5:119944/1‑134 (MQ=255)
                            ggctggTCGTGTCTACAAGTCCCAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGAtt                                                                                        >  6:104864/1‑135 (MQ=255)
                               tggTCGTGTCTACAAGTCCCAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGTGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGCCTGGATTGAtt                                                                                    >  6:107607/1‑135 (MQ=255)
                               tggTCGTGTCTACAAGTCCCAGCCGGCGCCTCTCAACCTGCCAGTGCCACCGCACTCTTATGAGTCGGACGAGATCATAGTCATTGGTGTGGCCCTGGACAACGCCAGATGAGCACGTTCATGGACTGGATTGATg                                                                                    >  6:110089/1‑136 (MQ=255)
                               tggTCGTGTCTACAAGTCCCAGCCGGCGCCTCTCAACCAGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGCTg                                                                                    >  3:57249/1‑136 (MQ=255)
                                       tCTACAAGTCCCAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCtt                                                                             <  2:5705/135‑1 (MQ=255)
                                          aCAAGTCCCAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGCCCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCACGCTTcccc                                                                         >  2:111966/1‑136 (MQ=255)
                                                            ctctCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCAGCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGaa                                                       <  2:33470/136‑1 (MQ=255)
                                                                 aaCCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCTTTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGac                                                    <  4:58900/134‑1 (MQ=255)
                                                                 aaCCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGacca                                                  >  5:92678/1‑136 (MQ=255)
                                                                       cccGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTc                                              <  1:154153/134‑1 (MQ=255)
                                                                       cccGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAg                                            <  4:118567/136‑1 (MQ=255)
                                                                         cGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGc                                           >  5:97700/1‑135 (MQ=255)
                                                                           tGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGAGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTCATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAg                                        >  1:102547/1‑136 (MQ=255)
                                                                             ccACCGCACTCTTATGCGTCGGACGAGATCATCTTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTa                                      <  6:5500/136‑1 (MQ=255)
                                                                              cACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTAt                                     >  3:9537/1‑136 (MQ=255)
                                                                               aCCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTAtt                                    <  1:25244/136‑1 (MQ=255)
                                                                               aCCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTAtt                                    <  6:119944/136‑1 (MQ=255)
                                                                                ccGCTCTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAg                                            <  1:131740/127‑1 (MQ=255)
                                                                                ccGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAg                                            >  2:131740/1‑127 (MQ=255)
                                                                                                 ggACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTGCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTATTACGCACCCAAGaacttca                  <  4:160253/136‑1 (MQ=255)
                                                                                                 ggACGAGATCATCGTCATTGGCGGCGAACAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGGTGTGGGAAGACCACCTCAGCAAGTATTACGCACCCTAGaacttca                  <  2:152964/136‑1 (MQ=255)
                                                                                                     gagaTCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTATTACGCACCCAAGaacttcaactt              >  1:9465/1‑136 (MQ=255)
                                                                                                     gagaTCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTATTACGCACCCAAGaact                     >  7:73057/1‑129 (MQ=255)
                                                                                                     gagaTCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTATTACGCACCCAAGaact                     <  8:73057/129‑1 (MQ=255)
                                                                                                                cATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTATTACGCACCCAAGAACTTCAACTTCCTGTActtct   >  2:20016/1‑136 (MQ=255)
                                                                                                                 aTTGGTGTCGACCAGGATAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTATTACGCACCCAAGAACTTCAACTTCCTGTActtctt  <  5:110089/136‑1 (MQ=255)
                                                                                                                 aTTGGTGACGACCATGAGAACGCATGATGACCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTATTACGCACCCAAGAACTTCAACTTCCTGTActtctt  <  5:107607/136‑1 (MQ=255)
                                                                                                                      |                                                                                                                                  
GCCCGTGCCACGGCTCTCACTATGACCTGGCTGGTCGTGTCTACAAGTCCCAGCCGGCGCCTCTCAACCTGCCCGTGCCACCGCACTCTTATGAGTCGGACGAGATCATCGTCATTGGCGTCGACCAGGAGAACGCATGATGAGCAAGTTCATGGACTGGATTGATGCTCGCTTCCCCGCCACCAAGATGTGGGAAGACCACCTCAGCAAGTATTACGCACCCAAGAACTTCAACTTCCTGTACTTCTT  >  ACB122/1497806‑1498054

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