Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,310,190 T→G 55.6% noncoding (2370/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,1900TG55.6% ‑6.1 / 20.3 18noncoding (2370/2902 nt)PP_23SG23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base G (2/8);  total (5/13)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGTTACTCTTTGGGAGGAGACCGCCCCAGTCAAACTAC  >  NC_002947/5310047‑5310332
                                                                                                                                               |                                                                                                                                              
cctggATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTc                                                                                                                                          <  1:2536936/147‑1 (MQ=1)
cctggATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTc                                                                                                                                          >  2:2519133/4‑150 (MQ=1)
cctggATGTGATGAGCCGACAGCGAGGTGCCACACACCGCCGTCGATATGAACTCTTGGGCGGTTTCTGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTc                                                                                                                                          <  2:5649341/147‑1 (MQ=1)
    gttgcgatgaGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGAGTGTTGGGGGTGATCAGCCTGTTATCCCCGGCGTAGCTTTTATCCGTTGAGCGATGGGCCTTGCACGCGGGACCACCGGATGACTATGTCCGTCTTTCGtct                                                                                                                                      >  2:6580598/6‑148 (MQ=0)
      tgcgatgaGTCGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGtctct                                                                                                                                    <  1:5657304/147‑5 (MQ=1)
             agtcGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGAc                                                                                                                             <  1:1060911/147‑1 (MQ=1)
             agtcGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGAc                                                                                                                             <  2:4584520/147‑1 (MQ=1)
             agtcGACATCGAGGTGCCAAACAACCCCGTCGATATGGACTCTTGGGGGTCATCAGCCTGTTATCCCCGGCGTACCTTTTATCCGTTGAGCGATGGCCCTTCCACGCGGGACCACCGGATCACTATGACCGACTTTCGTCTCTGCTCGAc                                                                                                                             <  2:5724984/147‑1 (MQ=1)
                                                 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:3003291/147‑1 (MQ=0)
                                                                                                                                        tatgACCGACTTTCGTCTCTGCTCGACTTGTCAGTCTTGCAGTCAGGCAGGCTTATGCCATTGCACTCGACGAACGATTTCCGACCGTTCTGAGCCTACCATCGCGCGCCTCCGTTACTCTTTAGGAGGCGACCGCCCCAGTCAAACTAc  <  1:1892727/147‑1 (MQ=0)
                                                                                                                                               |                                                                                                                                              
CCAGGATGTGATGAGCCGACATCGAGGTGCCAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATACAGAACCACCGGATCACTAAGACCTACTTTCGTACCTGCTCGACGTGTTTGTCTCGCAGTCAAGCGCGCTTTTGCCTTTATACTCTACGACCGATTTCCGACCGGTCTGAGCGCACCTTCGTACTCCTCCGTTACTCTTTGGGAGGAGACCGCCCCAGTCAAACTAC  >  NC_002947/5310047‑5310332

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