Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,431,349 C→A 26.8% V460V (GTC→GTA PP_1252 → group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,431,3490CA26.8% 13.0 / 4.9 15V460V (GTC→GTAPP_1252group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base C (8/3);  new base A (2/2);  total (10/5)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CCAAAGAAACTGTGGGACCGGCTCGGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCAGCGCAACCCTCAGGCTTAT  >  NC_002947/1431317‑1431492
                                |                                                                                                                                               
ccAAAGAAACTGTGGGACCGGCTCGGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGcc                            <  2:2095732/150‑1 (MQ=11)
 cAAAGAAACTGTGGGACCGGCTCGGGCTGGTCTCGATACTGGATACGATTAACCCGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCa                           >  1:1488581/1‑150 (MQ=14)
  aaaGAAACTGTGGGACCGGCTCGGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCaa                          >  2:1246992/1‑150 (MQ=32)
  aaaGAAACTGTGGGACCGGCTCGGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGAACGGTGGTGTGAGAGAACTGCGAGTGCGAATCCGCCTCCTACTCGATTGccca                          >  2:623864/1‑148 (MQ=11)
      aaaCTGTGGGACCGGCTCGGGCTGGTATCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGccc                      >  2:571935/1‑150 (MQ=25)
      aaaCTGTGGGACCGGCTCGGGCTGGTATCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGccc                      >  2:580551/1‑150 (MQ=25)
       aaCTGTGGGACCGGCTCGGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCa                     <  1:865115/150‑1 (MQ=35)
                  cGGCTCGGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGCGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCAGCGCAACCCTc          >  2:1441057/1‑150 (MQ=255)
                    gCTCGGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCCGCGCAACCCCCac        >  2:60179/1‑149 (MQ=255)
                    gCTCGGGCTGGTATCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCAGCGCAACCCTCAg        <  1:571935/150‑1 (MQ=255)
                    gCTCGGGCTGGTATCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCAGCGCAACCCTCAg        <  1:580551/150‑1 (MQ=255)
                        gggCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCAGCGCAACCCTCAGGCtt    >  1:2037412/1‑150 (MQ=255)
                          gCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCAGCGCAACCCTCAGGCTTAt  >  1:2215902/1‑150 (MQ=255)
                              gTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCAGCGCAAcc             <  1:2111720/135‑1 (MQ=255)
                              gTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCAGCGCAAcc             >  2:2111720/1‑135 (MQ=255)
                                |                                                                                                                                               
CCAAAGAAACTGTGGGACCGGCTCGGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCCAAGCCCAGCGCAACCCTCAGGCTTAT  >  NC_002947/1431317‑1431492

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