Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,180,277 T→C R94R (CGT→CGC rplJ → 50S ribosomal subunit protein L10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,180,2770TC97.1% 111.1 / NA 36R94R (CGT→CGCrplJ50S ribosomal subunit protein L10
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/15);  minor base G (1/0);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTTCGTAACACCCTGCTGCGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGTCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGACCTACGAAGAAGCAATTGCACGC  >  NC_000913/4180155‑4180409
                                                                                                                          |                                                                                                                                    
gttgttCGTAACACCCTGCTGCGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTcgc                                                                                                                     >  1:772868/1‑140 (MQ=255)
        tAACACCCTGCTGCGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCaaa                                                                                                                             <  2:572619/124‑1 (MQ=255)
        tAACACCCTGCTGCGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCaaa                                                                                                                             >  1:572619/1‑124 (MQ=255)
              ccTGCTGCGCAGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACTCGGGCGCTGCTGATCGGCTGTTCAACGAGTACGCGAAAGCGAATGCaa                                                                                                       >  2:308305/1‑140 (MQ=255)
                tgctgcGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTGGCCAGTTCAAAGAGTTCTCGAAAGCGAATGCaaaa                                                                                                     >  2:313692/1‑140 (MQ=255)
                  ctgcGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGGTCGACTGTTCAAAGAGTTCGCAAAAGCGAATGCAAAAtt                                                                                                   >  2:273432/1‑140 (MQ=255)
                  ctgcGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAAtt                                                                                                   >  2:98946/1‑140 (MQ=255)
                  ctgcGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAAtt                                                                                                   >  1:776715/1‑140 (MQ=255)
                   tgcGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAAttt                                                                                                  <  1:98946/140‑1 (MQ=255)
                       ccGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAgg                                                                                              >  2:719329/1‑140 (MQ=255)
                       ccGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGt                                                                                             >  1:706337/1‑140 (MQ=255)
                          tGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCa                                                                                           >  2:670573/1‑140 (MQ=255)
                                tGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCg                                                                                     <  1:493504/140‑1 (MQ=255)
                                tGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCg                                                                                     <  2:772868/140‑1 (MQ=255)
                                  aaGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCt                                                                                   <  2:104864/140‑1 (MQ=255)
                                       aCTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAGAGAGTTCGAGAAGGCGGATGCAAAAGTTGACGTCAAAGCCGCTGCCtt                                                                              >  2:214353/1‑140 (MQ=255)
                                            gTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAg                                                                         <  2:776715/140‑1 (MQ=255)
                                                       tGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATc                                                              >  2:636761/1‑140 (MQ=255)
                                                          aaGACGCGTTTGTTGCTCCGACGCTCATTACATACTCCATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCg                                                           <  1:214353/140‑1 (MQ=255)
                                                          aaGACGCGTTTGTTGCTCCGACCCTAATTCCATACTCGATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCg                                                           <  1:313692/140‑1 (MQ=255)
                                                            gACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGc                                                         <  1:670573/140‑1 (MQ=255)
                                                                           ccGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCg                                          <  1:719329/140‑1 (MQ=255)
                                                                           ccGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCg                                          >  1:668449/1‑140 (MQ=255)
                                                                                       gCATACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTg                             <  1:461486/140‑1 (MQ=255)
                                                                                        cATACTCTATGGAACGCCCGGGCGCCGCTGCTCGCCTCTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTg                             <  1:273432/140‑1 (MQ=255)
                                                                                         aTACTCTATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGc                            >  1:81021/1‑140 (MQ=255)
                                                                                            ctctATGGAACACCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGAc                        <  2:706337/140‑1 (MQ=255)
                                                                                              ctATGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCACAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGAcc                       >  2:763245/1‑140 (MQ=255)
                                                                                                 tGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAACATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGGGTCTCAGTTCGGCCACCTGGCCACTCTGCCGACCTAc                    >  1:247761/1‑140 (MQ=255)
                                                                                                 tGGAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGACCTAc                    >  1:631317/1‑140 (MQ=255)
                                                                                                  ggAACACCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGACCTACg                   >  2:468580/1‑140 (MQ=255)
                                                                                                       acCCGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGACCTACgaagaa              >  2:637448/1‑140 (MQ=255)
                                                                                                         ccGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGACCTACGAAGAAGc            >  1:195683/1‑140 (MQ=255)
                                                                                                         ccGGGCGCTGCTGCTCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCAGGCGTCTCAGATTGACCGCCTGGCACCTCTACCGACCTACGAAGAAGc            >  2:303188/1‑140 (MQ=255)
                                                                                                               gctgctgctCGCCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGACCTACGAAGAAGCAATTGc      <  2:100978/140‑1 (MQ=255)
                                                                                                                   ctgctCGCCTGTTCAAAGAGTTCGCGAAATCGAATGCCAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGACCTACGAAGAAGCAATTGCACGc  <  1:308305/140‑1 (MQ=255)
                                                                                                                          |                                                                                                                                    
GTTGTTCGTAACACCCTGCTGCGCCGTGCTGTTGAAGGTACTCCGTTCGAGTGCCTGAAAGACGCGTTTGTTGGTCCGACCCTGATTGCATACTCTATGGAACACCCGGGCGCTGCTGCTCGTCTGTTCAAAGAGTTCGCGAAAGCGAATGCAAAATTTGAGGTCAAAGCCGCTGCCTTTGAAGGTGAGCTGATCCCGGCGTCTCAGATCGACCGCCTGGCAACTCTGCCGACCTACGAAGAAGCAATTGCACGC  >  NC_000913/4180155‑4180409

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