Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,310,242 T→A 65.5% noncoding (2318/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,310,2420TA65.5% ‑5.9 / 7.8 17noncoding (2318/2902 nt)PP_23SG23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base A (7/4);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGTTACTCTTTGGGAGGAGACCGCCCCAGTCAAACTACCCACCATACACTGTCCTCGATCCGGATAACGGACCTGAGTTAGAACCTCAAA  >  NC_002947/5310096‑5310384
                                                                                                                                                  |                                                                                                                                              
gatctctTGGGCGGTATCAGCCTGTTATCCCCGGAGTAGCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTta                                                                                                                                             >  2:1919709/4‑150 (MQ=1)
                  agacTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTc                                                                                                                           <  1:7651898/147‑1 (MQ=1)
                                ggcgTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCTGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAgc                                                                                                             <  2:8243483/147‑1 (MQ=0)
                                ggcgTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCt                                                                                                                 >  1:4385741/4‑146 (MQ=0)
                                ggcgTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCt                                                                                                                 <  2:4385741/143‑1 (MQ=0)
                                  aggaCCTTTTAGCCGTTGCGCGACGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAgcgc                                                                                                           <  2:5133511/147‑1 (MQ=0)
                                                   gaccGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATAATCTACGTCCGATTTCCGACCGGTCTGAGCGCTCCTTCGTACTCCTCCg                                                                                          >  2:8216657/4‑150 (MQ=0)
                                                             ccgtCCAGACAGCACCCCCGGATCACGAAGACCTACTTTCGTACCGGCTCGACGTGTATGTCTCGCAGTCAAGCGCGCGTTTGCCTTCATACTCTCCGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGGTACTCttgg                                                                                <  2:3505395/147‑3 (MQ=1)
                                                                                       tatgACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTAc                                                      <  1:1892727/147‑1 (MQ=0)
                                                                                       tatgACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTAc                                                      <  1:3003291/147‑1 (MQ=0)
                                                                                            ccgaCTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACccacc                                                 >  2:2530519/4‑150 (MQ=0)
                                                                                            ccgaCTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACccacc                                                 >  1:6293344/4‑150 (MQ=0)
                                                                                            ccgaCTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACccacc                                                 >  1:3207860/4‑150 (MQ=0)
                                                                                            ccgaCTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACccacc                                                 >  1:1111024/4‑150 (MQ=0)
                                                                                            ccgaCTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACccacc                                                 >  2:827206/4‑150 (MQ=0)
                                                                                            ccgaCTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTACccacc                                                 >  2:833160/4‑150 (MQ=0)
                                                                                                                                           ttgtGCCTTTATACTCTTCGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGTTACTCTTTGGGAGGAGACCGCCCCAGTCAAACTACCCACCATACACTGTCCTCGATCCGGATAACGGACCTGAGTTAGAACCTCaaa  <  1:7646003/147‑1 (MQ=1)
                                                                                                                                                  |                                                                                                                                              
GAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGTTACTCTTTGGGAGGAGACCGCCCCAGTCAAACTACCCACCATACACTGTCCTCGATCCGGATAACGGACCTGAGTTAGAACCTCAAA  >  NC_002947/5310096‑5310384

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: