Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,309,871 T→C 73.0% noncoding (2689/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,8710TC73.0% ‑5.5 / 16.8 32noncoding (2689/2902 nt)PP_23SG23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (3/2);  major base C (13/7);  minor base A (2/5);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

ACACACCCAACCTATCAACGTCGTAGTCTTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCAC  >  NC_002947/5309722‑5310004
                                                                                                                                                     |                                                                                                                                     
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                                      >  2:7953280/4‑150 (MQ=0)
                                                                                                 gccttCAGCGGTTATCCCGTCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc                                      <  1:627228/147‑1 (MQ=0)
                                                                                                 gccttCAGCGGTTATCCCGTCCGTATATAGCTAACCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAAc                                      >  2:3337634/4‑150 (MQ=0)
                                                                                                      cagcacttATCCCTTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct                                 <  1:7072608/144‑1 (MQ=0)
                                                                                                      cagcacttATCCCTTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct                                 <  2:6230280/144‑1 (MQ=0)
                                                                                                      cagcacttATCCCTTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGtct                                 <  2:6230250/144‑1 (MQ=0)
                                                                                                              atcccgtCCGTATATAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGt                         >  2:7348559/7‑150 (MQ=0)
                                                                                                               tccctTCCGCACATAGCTACCCAGCGATGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGtt                        <  2:4545261/147‑1 (MQ=0)
                                                                                                                     ccgtataTAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACt                                     <  1:6464269/125‑1 (MQ=0)
                                                                                                                     ccgtataTAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACt                                     >  2:6464269/7‑131 (MQ=0)
                                                                                                                     ccgtataTAGCTACCCTGCTATGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAAc                  <  1:6886888/144‑1 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  2:7255437/7‑150 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  2:5946123/7‑150 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    <  2:3332196/144‑1 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  2:5842217/7‑150 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  2:5408710/7‑150 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    <  1:836172/144‑1 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    <  1:8223104/144‑1 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  1:8143102/7‑150 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  1:4162928/7‑150 (MQ=0)
                                                                                                                                   cctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAcc    >  1:154373/7‑150 (MQ=0)
                                                                                                                                   ccagcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTAcc    >  1:1469488/7‑150 (MQ=0)
                                                                                                                                   ccagcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTAcc    >  2:281772/7‑150 (MQ=0)
                                                                                                                                    ctgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAccg   >  1:4177393/6‑149 (MQ=0)
                                                                                                                                     tgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAccgc  >  1:852099/5‑148 (MQ=0)
                                                                                                                                     tgctaTGCCGTTGGCACGACAACAGGTCCACCAGAGATACGTCCATCCCGGTCCTCTCGTACTAGGGACAGATCCTATCAATATTCCTACACCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAACTCACGTAccgc  <  2:7050878/146‑3 (MQ=0)
                                                                                                                                     agcgaTGCCTTTGGCAAGACAACTGGTACACCAGCGGTGCGTCCATCCCGGTCCTCTCGTACTAAGGACAGCTCCTCTCAAATTTCCTACGCCCACGACGGATAGGGACCGAACTGTCTCACGACGTTCTGAACCCAGCTCGCGTAccgc  <  1:410333/146‑3 (MQ=0)
                                                                                                                                                     |                                                                                                                                     
ACACACCCAACCTATCAACGTCGTAGTCTTCGACGGCCCTTTAGGGGATTCAAGATCCCAGTGAGATCTCATCTTGAGGCAAGTTTCCCGCTTAGATGCTTTCAGCGGTTATCTCTTCCGAACATAGCTACCCGGCAATGCCACTGGCGTGACAACCGGAACACCAGAGGTTCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCTCTCAAATCTCAAACGTCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCAC  >  NC_002947/5309722‑5310004

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: