Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,718,111 G→C 14.5% 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,1110GC14.5% 129.2 / 15.1 49L256F (TTG→TTCrsmJribosomal RNA small subunit methyltransferase J
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base C (3/4);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

GCGCTGCTCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCA  >  NC_002947/1717965‑1718260
                                                                                                                                                  |                                                                                                                                                      
gcgcTGCTCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAg                                                                                                                                                     <  2:913792/150‑1 (MQ=255)
  gctgctCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGc                                                                                                                                                   >  2:543220/1‑150 (MQ=255)
   ctgctCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGcc                                                                                                                                                  <  1:268183/150‑1 (MQ=255)
    tgctCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCt                                                                                                                                                 <  1:873542/150‑1 (MQ=255)
       tCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGAc                                                                                                                                              >  1:650555/1‑150 (MQ=255)
       tCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGAc                                                                                                                                              >  2:1244355/1‑150 (MQ=255)
          aaGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttcttcg                                                                                                                                           <  1:543220/150‑7 (MQ=255)
          aaGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttcttcg                                                                                                                                           <  2:1294409/150‑7 (MQ=255)
             ccGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGccc                                                                                                                                        <  1:680024/150‑1 (MQ=255)
                 tctGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTAt                                                                                                                                    <  2:650555/150‑1 (MQ=255)
                 tctGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttcttcg                                                                                                                                           <  2:1186743/143‑7 (MQ=255)
                 tctGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttcttcg                                                                                                                                           >  1:1186743/1‑137 (MQ=255)
                       gcTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGcc                                                                                                                              >  2:868969/1‑150 (MQ=255)
                         tGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGcccc                                                                                                                                       >  2:509859/1‑139 (MQ=255)
                         tGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGcccc                                                                                                                                       <  1:509859/139‑1 (MQ=255)
                              agccaCCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGc                                                                                                                       >  2:622786/1‑150 (MQ=255)
                                     gggTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTg                                                                                                                >  1:1089099/1‑150 (MQ=255)
                                           tAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAg                                                                                                         >  2:1235427/1‑150 (MQ=255)
                                                      ccACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTg                                                                                                                                           <  1:535649/106‑1 (MQ=255)
                                                      ccACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTg                                                                                                                                           >  2:535649/1‑106 (MQ=255)
                                                      ccACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTgg                                                                                                                                >  1:543599/1‑117 (MQ=255)
                                                      ccACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTgg                                                                                                                                <  2:543599/117‑1 (MQ=255)
                                                      ccACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGcc                                                                                              >  1:385182/1‑150 (MQ=255)
                                                                cACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCg                                                                                    >  1:67019/1‑150 (MQ=255)
                                                                   cGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTg                                                                                                                                 >  1:1365315/1‑103 (MQ=255)
                                                                   cGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTg                                                                                                                                 <  2:1365315/103‑1 (MQ=255)
                                                                                gCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCAt                                                                    >  1:415163/1‑150 (MQ=255)
                                                                                  cGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttct                                                                                                                                              >  2:224502/1‑72 (MQ=255)
                                                                                  cGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTCAAGGCCttct                                                                                                                                              <  1:224502/75‑4 (MQ=255)
                                                                                      gCCGAGCCATAGCCGTGAGGGTAAGTCGAGCCGGTGTGACATCTACCCGAAGAAGGCGTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGAGTGTAGGAGCGGATATAGCGGCGAACACCGACAACGTCGGTGCCATacacgg                                                              >  2:101569/1‑148 (MQ=255)
                                                                                               tAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGtctc                                                     >  1:698317/1‑150 (MQ=255)
                                                                                               tAGCCTTGAGGGTAAGTCGAGCAGGTTGGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGtctc                                                     >  2:1194923/1‑150 (MQ=255)
                                                                                                                  aGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGCCGTCTCCTGCGCGTTAAACCTGCTc                                  <  1:1235427/150‑1 (MQ=255)
                                                                                                                     cGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTa                               <  1:622786/150‑1 (MQ=255)
                                                                                                                       gTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACa                             <  1:868969/150‑1 (MQ=255)
                                                                                                                       gTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCACGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGt                                             <  2:1273207/134‑1 (MQ=255)
                                                                                                                       gTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCACGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGt                                             >  1:1273207/1‑134 (MQ=255)
                                                                                                                         aTGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAAt                           <  2:67019/150‑1 (MQ=255)
                                                                                                                            aCATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGa                         >  2:321199/1‑149 (MQ=255)
                                                                                                                            aCATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGa                         <  1:321199/149‑1 (MQ=255)
                                                                                                                            aCATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGAc                        <  2:698317/150‑1 (MQ=255)
                                                                                                                            aCATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGAc                        <  2:895375/150‑1 (MQ=255)
                                                                                                                              aTCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCt                      <  2:385182/150‑1 (MQ=255)
                                                                                                                                   cccGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGc                                    >  1:707437/1‑131 (MQ=255)
                                                                                                                                   cccGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGc                                    <  2:707437/131‑1 (MQ=255)
                                                                                                                                   cccGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTAt                 >  1:821402/1‑150 (MQ=255)
                                                                                                                                       aagaagGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCACCCCCGGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCt             <  1:1244355/150‑1 (MQ=255)
                                                                                                                                         gaagGCCTTTACGTCCTGACCTGCCCCTATCTGGCAGGCTAGCCCGAATGAAACGACAGGTCATGGTTGGTACTCCTAAAGCGGCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGc           <  1:101569/150‑1 (MQ=255)
                                                                                                                                                 tGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTc   <  1:1194923/150‑1 (MQ=255)
                                                                                                                                                  gAAGGCCTGACCTGCCCCTATCTGGCCGGCCAGCACGATTGTAGGGGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCa  >  1:275011/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GCGCTGCTCGAAGCCGCTCTGGCGCTGGCCAGCCACCGGGTGGTAGTGAAGCGCCCACGCAAGGCACCGATCATCGACGGGCCGAAGCCGAGCCATAGCCTTGAGGGTAAGTCGAGCCGGTATGACATCTACCCGAAGAAGGCCTTGAAGGCCTGACCTGCCCCTATCTGGCCGGCTAGCACGATTGTAGGAGCGGCTTTAGCCGCGAACACCGACAACGCCGGTGCCATCCACTGCGTCGTCTCCTGCGCGTTAAACCTGCTCCTACAATGACCTGTTATTTCTGCCAGTAGTCA  >  NC_002947/1717965‑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: