Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,729,508 T→C 69.3% noncoding (306/1542 nt) rrsG ← 16S ribosomal RNA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077792,729,5080T→C69.3% 5.6 0.0 ‑0.213noncoding (306/1542 nt)rrsG16S ribosomal RNA
Reads supporting (aligned to +/- strand):  new base (3/6):  ref base (2/2):  total (5/8)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.07e-01

GTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGG  >  NC_007779/2729442‑2729608
                                                                  |                                                                                                    
gtgCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAATGTGGCCGTCCACCCTCTCAGGCCGGCTACCCGTCGCCGCCTTg                                                                    <  1:12722532/101‑1 (MQ=0)
gtgCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAATGTGGCCGTCCACCCTCTCAGGCCGGCTACCCGTCGCCGCCTTg                                                                    >  1:2674951/1‑101 (MQ=0)
gtgCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAATGTGGCCGTCCACCCTCTCAGGCCGGCTACCCGTCGCCGCCTTg                                                                    >  1:3787293/1‑101 (MQ=0)
gtgCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAATGTGGCCGTCCACCCTCTCAGGCCGGCTACCCGTCGCCGCCTTg                                                                    <  1:6605986/101‑1 (MQ=0)
gtgCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAATGTGGCCGTCCACCCTCTCAGGCCGGCTACCCGTCGCCGCCTTg                                                                    >  1:8976554/1‑101 (MQ=0)
gtacAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAATGTGGCCGTCCACCCTCTCAGGCCGGCTACCCGTCGCCGCCTTg                                                                    <  1:4958500/98‑1 (MQ=0)
          ccccACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAATGTGGCCGTCCACCCTCTCAGGCCGGCTACCCGTCGCCGCCTTGGtaggccatt                                                          <  1:13771286/101‑9 (MQ=0)
           cccACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAATGTGGCCGTCCACCCTCTCAGGCCGGCTACCCGTCGCCGCCTTGGTAGGCCATTa                                                         <  1:4102141/101‑1 (MQ=0)
           cccACTGCTGCCTCCCGTAGGAGTCTGGGCCGTGTCTCAGTCCCAATGTGGCCGTCCACCCTCTCAGGCCGGCTACCCGTCGCCGCCTTGGTAGGCCATTa                                                         <  1:5622340/101‑1 (MQ=0)
                                                                ggTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa    <  1:11268751/101‑1 (MQ=11)
                                                                 gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg   <  1:14147596/101‑1 (MQ=31)
                                                                 gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAg   >  1:3563379/1‑101 (MQ=31)
                                                                 gTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTAGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGaa  >  1:7683930/1‑100 (MQ=11)
                                                                  tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTCCCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAgg  <  1:2572836/101‑1 (MQ=17)
                                                                  tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAgg  <  1:9782546/101‑1 (MQ=32)
                                                                  |                                                                                                    
GTGCAATATTCCCCACTGCTGCCTCCCGTAGGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGGCAAGAGGCCCGAAGG  >  NC_007779/2729442‑2729608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.