Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,447,639 G→A 100% H268Y (CAT→TAT)  prmB ← 50S ribosomal subunit protein L3 N(5)‑glutamine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,447,6390GA92.6% 38.5 / ‑4.7 22H268Y (CAT→TAT) prmB50S ribosomal subunit protein L3 N(5)‑glutamine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base A (8/10);  minor base C (2/0);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGT  >  NC_000913/2447499‑2447779
                                                                                                                                            |                                                                                                                                            
cgttgttACCTAATCTGTATAAATCGCGAAGTGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGAtata                                                                                                                                            <  2:360058/138‑1 (MQ=255)
         ttAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGCGAACGGAACATCCGGATATTGTTCCATAAGATCTACCATGCTGt                                                                                                                                   >  2:372145/1‑143 (MQ=255)
          tAAGCTTTATAAATCGCGAGATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATATACCATGCTGtt                                                                                                                                  <  2:419865/143‑1 (MQ=255)
                   taAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAAATGTACCATCCTCTTGCCGACTTc                                                                                                                         >  1:617959/1‑143 (MQ=255)
                         gcgAAATGTTCTCGTGCGGCAATAAGCGGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAAt                                                                                                                   <  1:220115/143‑1 (MQ=255)
                            aaaTGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAATCaa                                                                                                                >  1:128613/1‑143 (MQ=255)
                               tGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATAGCTCCATAAGATATACCATGCTGTTGCCGACCTCACAACTccagcc                                                                                                             >  2:435443/1‑138 (MQ=255)
                               tGTTCTCGTGCGGCAATAAGCTGCTCTGTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAATCAacac                                                                                                             <  1:323202/143‑1 (MQ=255)
                                        gCGGCAATAAGGTGCTCGTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGAAGAACGGAACACCCGGATATTGTTCCATAAGAAATACCACGCTGGTGCCGACTTCACAAATCAACACGCcaccctg                                                                                                    >  2:631641/1‑138 (MQ=255)
                                                   ctgctCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATCTACCATGCTGTTGCCGACTTCCAAAATCAACACGCAATCACCAGCAAGGTacc                                                                                         >  2:543169/1‑141 (MQ=255)
                                                           tGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCcgcg                                                                                 <  1:346355/143‑1 (MQ=255)
                                                                                      cgTTATCAAATGCCAGCCAGTTAAAGGTAACATCCGGATAGTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTc                                                      <  2:380301/143‑1 (MQ=255)
                                                                                       gTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCGGTTTCa                                                     >  1:2148/1‑143 (MQ=255)
                                                                                        ttATCAAACGCCAGCCACGTGAACGGAACCTCCGGATATTGTTCCATAAGATATACCATGCTGTTGCCGACGTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAg                                                    <  2:141302/143‑1 (MQ=255)
                                                                                                            gAACGGAACATCCGGATATTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCAGGAATGCGACGCGGCAGTTTCAGGCGGGCAGTGCCCGAGGCCa                                >  2:537705/1‑143 (MQ=255)
                                                                                                               cGGAACATCCGGATATTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCCGTTTCAGGCCGTCCGTCCCAGATGCCAGGc                             >  2:41376/1‑143 (MQ=255)
                                                                                                               cGGAACATCCGGATATTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGGCAGTTGCAGGCCGTCAGTGCCAGGTGCGAGGc                             >  2:41408/1‑143 (MQ=255)
                                                                                                                ggAACATCCGAATATTGTTCCATAAGATATACCATGCTGCTGCCGACTTCACAAATCACCACGCCAGCATCATCAAGGTAATCTGCCTCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGcc                            <  1:480994/143‑1 (MQ=255)
                                                                                                                     aTCCGGAGATTGTTCCATAAGTGATACCATGCTGTTGCCGACTTCACAAATCAACACGCCAGCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGtt                       <  1:581522/143‑1 (MQ=255)
                                                                                                                       ccGGGTATTGTGCCATCTGGTATTCCATGCTGTTGCCGACTTCCCACATCACCACGCCATCTTCTGCGAGGTGGTCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGTCAGATGCCAGGCCCAGTTcc                     <  2:283053/143‑1 (MQ=255)
                                                                                                                          gATATTGTTCCATAAGATATACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATGTGCCGCGTTACCGAGCATGCGCCGCGGCGGTGTCAGGCCGTCAGTGCGAGATGCCGAGCCCAGTTCCggg                  >  2:207297/1‑142 (MQ=255)
                                                                                                                                  tCCATAAGATATACCATGCTGGTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCg          <  1:66511/143‑1 (MQ=255)
                                                                                                                                          aTGTGCCATGCGTTCGCCCAGTTCACAAGTCCACACGCCCTCCTCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGt  <  2:428776/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
CGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAACATCCGGATATTGTTCCATAAGATGTACCATGCTGTTGCCGACTTCACAAATCAACACGCCATCATCAGCAAGGTAATCTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGT  >  NC_000913/2447499‑2447779

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