Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,007,942 T→C G33G (GGT→GGC rlmL → fused 23S rRNA m(2)G2445 methyltransferase and 23S rRNA m(7)G2069 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,007,9420TC95.5% 66.3 / ‑4.5 22G33G (GGT→GGCrlmLfused 23S rRNA m(2)G2445 methyltransferase and 23S rRNA m(7)G2069 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (10/11);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCG  >  NC_000913/1007815‑1008077
                                                                                                                               |                                                                                                                                       
cGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTc                                                                                                                             <  1:335443/140‑1 (MQ=255)
 ggTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCa                                                                                                                            <  1:138214/140‑1 (MQ=255)
   tAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAg                                                                                                                          >  2:174841/1‑140 (MQ=255)
                    aGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGt                                                                                                         <  2:109322/140‑1 (MQ=255)
                                         tttGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAg                                                                                    >  2:204446/1‑140 (MQ=255)
                                         tttGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAg                                                                                    >  2:104296/1‑140 (MQ=255)
                                                 tACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGAGGTTCAGGGCGGGGACCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTgg                                                                            >  1:465061/1‑140 (MQ=255)
                                                   cGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGcc                                                                          <  1:160093/140‑1 (MQ=255)
                                                     gCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTc                                                                        <  1:204446/140‑1 (MQ=255)
                                                          tGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGta                                                                   >  2:221242/1‑140 (MQ=255)
                                                           gggCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGtat                                                                  <  2:465061/140‑1 (MQ=255)
                                                                           tAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGc                                                  >  2:254423/1‑140 (MQ=255)
                                                                               aactgaactgGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAgt                                              <  1:254423/140‑1 (MQ=255)
                                                                                    aactgGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAg                                         >  1:148530/1‑140 (MQ=255)
                                                                                              ccTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCg                               >  1:153626/1‑140 (MQ=255)
                                                                                                   gggCCGTTGAATGCCAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTa                          <  1:221242/140‑1 (MQ=255)
                                                                                                                ccAGGTGGTTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCg             <  2:148530/140‑1 (MQ=255)
                                                                                                                      ggtTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTc       <  1:174841/140‑1 (MQ=255)
                                                                                                                      ggtTCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTc       >  1:127415/1‑140 (MQ=255)
                                                                                                                        ttCAGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAg     >  1:177671/1‑140 (MQ=255)
                                                                                                                         tCAGGGTGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAgg    <  1:33590/140‑1 (MQ=255)
                                                                                                                           aGGGCGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCg  <  1:104296/140‑1 (MQ=255)
                                                                                                                               |                                                                                                                                       
CGGTAGGCTCTTAGAGTAAAAGTGACGATATGAATTCTCTGTTTGCCAGTACGGCCCGTGGGCTGGAAGAGCTGTTAAAAACTGAACTGGAAAACCTGGGGGCCGTTGAATGCCAGGTGGTTCAGGGTGGGGTCCATTTCAAGGGCGACACACGGCTTGTTTACCAGAGCCTGATGTGGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCG  >  NC_000913/1007815‑1008077

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