Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,309,835 T→C 46.2% noncoding (2725/2902 nt) PP_23SG ← 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,309,8350TC46.2% ‑5.1 / 12.3 13noncoding (2725/2902 nt)PP_23SG23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base C (2/4);  total (5/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CGCCTCACAGCGCTTACACACCCAACCTATCAACGTCGTAGTCTTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCT  >  NC_002947/5309707‑5309973
                                                                                                                                |                                                                                                                                          
cgacTCACAGCGCTTACACACCCAACCTATCAACGTCGTAGTCTTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCgg                                                                                                                       <  1:6382829/147‑1 (MQ=1)
               acccacCCAACCTATCAACGTCGTAGTCTTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCGTAGATGCTTTCAGCGGTTATCTCTTCCGACCATAGCTACCCGGCAATGCCACTGGCGt                                                                                                        <  2:244265/147‑1 (MQ=0)
                                           tttgACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCa                                                                            >  2:3157443/4‑150 (MQ=1)
                                           ttagACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCa                                                                            >  2:4577448/4‑150 (MQ=1)
                                                   cccataaGGGGATTCAAGATCCCAGAGAGATCTCATCTTGAGGCAAGTTTCCCGCATAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTACGTCCACTCCGGTc                                                                    >  2:720086/7‑150 (MQ=1)
                                                                               gaactcATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCtctcaa                                        <  2:6834304/147‑1 (MQ=1)
                                                                                       ctagAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCtctcaaatctcaaa                                <  2:6687002/147‑1 (MQ=1)
                                                                                                                gccttCAGCGGTTATCCCGTCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc       <  1:627228/147‑1 (MQ=0)
                                                                                                                gccttCAGCGGTTATCCCGTCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc       >  2:7953280/4‑150 (MQ=0)
                                                                                                                gccttCAGCGGTTATCCCGTCCGTATATAGCTAACCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc       >  2:3337634/4‑150 (MQ=0)
                                                                                                                     cagcacttATCCCTTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct  <  1:7072608/144‑1 (MQ=0)
                                                                                                                     cagcacttATCCCTTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct  <  2:6230250/144‑1 (MQ=0)
                                                                                                                     cagcacttATCCCTTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct  <  2:6230280/144‑1 (MQ=0)
                                                                                                                                |                                                                                                                                          
CGCCTCACAGCGCTTACACACCCAACCTATCAACGTCGTAGTCTTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCT  >  NC_002947/5309707‑5309973

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