Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,718,111 G→C 11.6% L256F (TTG→TTC rsmJ → ribosomal RNA small subunit methyltransferase J

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,718,1110GC11.6% 122.9 / 7.1 43L256F (TTG→TTCrsmJribosomal RNA small subunit methyltransferase J
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base C (2/3);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCA  >  NC_002947/1717975‑1718260
                                                                                                                                        |                                                                                                                                                     
aaGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTg                                                                                                                                          >  1:215436/1‑150 (MQ=255)
 aGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGc                                                                                                                                         >  2:125694/1‑150 (MQ=255)
 aGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttcttcgg                                                                                                                                         <  1:715294/150‑8 (MQ=255)
                       caCCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACg                                                                                                                   >  1:423587/1‑150 (MQ=255)
                       caCCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCAACGCGTTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACg                                                                                                                   >  1:337086/1‑150 (MQ=255)
                        aCCGGGTGCTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGa                                                                                                                  >  2:596046/1‑150 (MQ=255)
                          cGGGTGCTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGAtt                                                                                                                <  1:596046/150‑1 (MQ=255)
                                        gcgcCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGc                                                                                                                          >  1:148468/1‑126 (MQ=255)
                                        gcgcCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGc                                                                                                                          <  2:148468/126‑1 (MQ=255)
                                        gcgcCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCttt                                                                                                  >  2:336296/1‑150 (MQ=255)
                                          ttaccctcAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAg                                                                                                >  2:411521/8‑150 (MQ=255)
                                                     gCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAATGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTTTAGTAGCGGCTTTAGCCGCGAACAcc                                                                                     <  2:269936/150‑1 (MQ=255)
                                                            tcatcGACGGGCCGAAGCCGAGCCATAGCCAGGAGGGTAAGTCGAGCCGGTATGACATCTAACCGAAGAAGGCGTTGAAGGCCTGACCTGCCCCTATCTGGCCCGCTAGCACGATTGTAGCAGCGGCATTAGCCGCGAACACCGACagcc                                                                              >  2:96703/1‑147 (MQ=255)
                                                               tcGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCg                                                                           >  2:790767/1‑150 (MQ=255)
                                                                 gACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGt                                                                         >  1:316312/1‑150 (MQ=255)
                                                                     ggCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCa                                                                     <  2:215436/150‑1 (MQ=255)
                                                                      gCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCAt                                                                    >  2:379318/1‑150 (MQ=255)
                                                                         gAAGCCTAGCCATAGCCTTGAGGGTAAGTCGAGTCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttc                                                                                                                                              <  1:764625/73‑3 (MQ=255)
                                                                         gAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGa                                                                                   <  2:266250/132‑1 (MQ=255)
                                                                         gAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGa                                                                                   >  1:266250/1‑132 (MQ=255)
                                                                         gAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttc                                                                                                                                              >  2:764625/1‑71 (MQ=255)
                                                                                                  aaGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAAc                                        <  2:337086/150‑1 (MQ=255)
                                                                                                  aaGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAAc                                        <  1:411521/150‑1 (MQ=255)
                                                                                                    gTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCt                                      <  1:125694/150‑1 (MQ=255)
                                                                                                      cGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGc                                    <  1:292812/150‑1 (MQ=255)
                                                                                                                tGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATg                          <  2:423587/150‑1 (MQ=255)
                                                                                                                 gACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCa                                                                 >  2:639485/1‑110 (MQ=255)
                                                                                                                 gACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCa                                                                 <  1:639485/110‑1 (MQ=255)
                                                                                                                         cccGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTAt                 <  1:790767/150‑1 (MQ=255)
                                                                                                                          ccGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGTCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTAtt                <  2:541569/150‑1 (MQ=255)
                                                                                                                            gaagaaGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTc              <  1:379318/150‑1 (MQ=255)
                                                                                                                               tcagGCCTTCAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGt                                             <  2:174135/114‑1 (MQ=255)
                                                                                                                               tcagGCCTTCAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGt                                             >  1:174135/3‑116 (MQ=255)
                                                                                                                               gaagGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGc           >  2:9027/1‑150 (MQ=255)
                                                                                                                                  ggCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGtt                   <  2:372469/139‑1 (MQ=255)
                                                                                                                                  ggCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGtt                   >  1:372469/1‑139 (MQ=255)
                                                                                                                                  ggCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCag        >  2:318712/1‑150 (MQ=255)
                                                                                                                                  ggCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCag        >  2:57479/1‑150 (MQ=255)
                                                                                                                                   gCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTcc                                 >  1:701266/1‑124 (MQ=255)
                                                                                                                                   gCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTcc                                 <  2:701266/124‑1 (MQ=255)
                                                                                                                                     cTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTg                     >  1:714911/1‑134 (MQ=255)
                                                                                                                                     cTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTg                     <  2:714911/134‑1 (MQ=255)
                                                                                                                                     cTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACAATCGAGACCGTTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCagtag     >  1:297441/1‑150 (MQ=255)
                                                                                                                                        gAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCa  >  1:686209/1‑150 (MQ=255)
                                                                                                                                        gAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCCTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCa  <  2:316312/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                     
AAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCA  >  NC_002947/1717975‑1718260

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