Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,429,568 G→A 20.0% intergenic (+78/‑402) mqo‑II → / → PP_1252 malate:quinone oxidoreductase/group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,429,5680GA20.0% 44.3 / 7.6 25intergenic (+78/‑402)mqo‑II/PP_1252malate:quinone oxidoreductase/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base G (8/12);  new base A (2/3);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

CAGCGTTGAATTCGGTGGTGCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGG  >  NC_002947/1429421‑1429674
                                                                                                                                                   |                                                                                                          
cAGCGTTGAATTCGGTGGTGCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTa                                                                                                          >  2:686269/1‑150 (MQ=255)
     ttGAATTCGGTGGTGCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCGTTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTGTGGGAGCGCCTGTCGTGTTGCGTTGCGCGTAtcagc                                                                                                     >  2:817107/1‑146 (MQ=255)
          ttCGGTGGTGCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAcgc                                                                                                <  2:167396/150‑3 (MQ=255)
                      gcaagcaaGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTCgtgg                                                                                    <  1:200380/150‑1 (MQ=255)
                      gcaagcaaGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGAGGCGCTGGTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtgg                                                                                    >  1:574362/1‑150 (MQ=37)
                        aagcaagcCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                                  <  1:817107/150‑1 (MQ=255)
                                 ggTCGAGAGAACGCCTGAGAACGCCATGTCGCTTGAAGCGGCCATTGTTGGCAAAGAAAGCCACGGGGGGGTCGTCGTGGCTTTTTTGTGTGCGCCAGGCTTGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                         <  2:227670/150‑1 (MQ=18)
                                        agCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                         <  1:44686/111‑1 (MQ=255)
                                        agCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                         >  2:44686/1‑111 (MQ=255)
                                           aaGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                               >  2:362454/1‑150 (MQ=255)
                                           aaGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                               <  1:362454/150‑1 (MQ=255)
                                            aGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTTCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                              <  2:276326/150‑1 (MQ=255)
                                             gCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAcgc                                                                                                <  1:577930/115‑3 (MQ=255)
                                             gCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAcgc                                                                                                >  2:577930/1‑113 (MQ=255)
                                                       gACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                                   <  1:664646/150‑1 (MQ=255)
                                                        aCATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                                  <  1:686269/150‑1 (MQ=255)
                                                                              tgttgtGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt                            <  2:287051/150‑1 (MQ=255)
                                                                                     aaagaaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGcc                     >  1:280639/1‑150 (MQ=255)
                                                                                        gaaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAc                                                                                                  <  1:754435/70‑1 (MQ=21)
                                                                                        gaaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAc                                                                                                  >  2:754435/1‑70 (MQ=21)
                                                                                              cACGGGGGAGACCTCGTGGCTTTGTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAg            >  2:672830/1‑150 (MQ=255)
                                                                                                       ccccccGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAGCCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTcgcg   <  1:258994/145‑1 (MQ=19)
                                                                                                        aCCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg  <  2:574362/150‑1 (MQ=38)
                                                                                                                gCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCg                                        >  2:182510/1‑104 (MQ=32)
                                                                                                                gCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCg                                        <  1:182510/104‑1 (MQ=32)
                                                                                                                                                   |                                                                                                          
CAGCGTTGAATTCGGTGGTGCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGG  >  NC_002947/1429421‑1429674

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: