Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,217,609 T→C 8.5% S228P (TCG→CCG)  moaA → cyclic pyranopterin monophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,217,6090TC 8.5% 124.5 / 2.9 47S228P (TCG→CCG) moaAcyclic pyranopterin monophosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (23/20);  new base C (2/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

GGTGCGCTTTGCCATCGAACGCGAGCTGGACATCACCTTCATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCACCGTCGAAGGCCGCCTGCTGCTGT  >  NC_002947/5217461‑5217756
                                                                                                                                                    |                                                                                                                                                    
ggTGCGCTTTGCCATCGAACGCGAGCTGGACATCACCTTCATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAAt                                                                                                                                                      >  7:134735/1‑149 (MQ=255)
      cTTTGCCATCGAACGCGAGCTGGACATCACCTTCATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCg                                                                                                                                                >  2:261025/1‑149 (MQ=255)
       tttGCCATCGAACGCGAGCTGGACATCACCTTCATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGAcc                                                                                                                                                 >  1:53745/1‑147 (MQ=255)
         tGCCATCGAACGCGAGCTGGACATCACCTTCATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAg                                                                                                                                              <  8:51902/148‑1 (MQ=255)
                            gACATCACCTTCATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTAGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCg                                                                                                                          <  1:214845/149‑1 (MQ=255)
                             aCATCACCTTCATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGc                                                                                                                         <  2:53745/149‑1 (MQ=255)
                                   ccTTCATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTgg                                                                                                                   <  4:157409/149‑1 (MQ=255)
                                       cATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGcc                                                                                                               >  3:267740/1‑149 (MQ=255)
                                         tCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTg                                                                                                             >  5:80357/1‑149 (MQ=255)
                                           gAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGc                                                                                                           >  1:128560/1‑149 (MQ=255)
                                           gAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGGCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGc                                                                                                           <  8:101909/149‑1 (MQ=255)
                                               aaaTGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGaa                                                                                                       <  7:23446/149‑1 (MQ=255)
                                               aaaTGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGaa                                                                                                       <  5:252925/149‑1 (MQ=255)
                                                     cGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCcgcg                                                                                                 <  8:134735/149‑1 (MQ=255)
                                                            ggTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACAccc                                                                                          >  8:32720/1‑149 (MQ=255)
                                                                    gCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGgc                                                                                                                  >  6:240411/1‑117 (MQ=255)
                                                                    gCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGgc                                                                                                                  <  5:240411/117‑1 (MQ=255)
                                                                            gAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAACCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCcgc                                                                                                  >  6:179571/1‑125 (MQ=255)
                                                                            gAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAACCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCcgc                                                                                                  >  6:184041/1‑125 (MQ=255)
                                                                            gAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAACCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCcgc                                                                                                  <  5:179571/125‑1 (MQ=255)
                                                                            gAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAACCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCcgc                                                                                                  <  5:184041/125‑1 (MQ=255)
                                                                              aCGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTccc                                                                        <  6:45748/149‑1 (MQ=255)
                                                                              aCGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTccc                                                                        <  4:150999/149‑1 (MQ=255)
                                                                                gTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑ca                                                                     <  3:217800/149‑1 (MQ=255)
                                                                                         cGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGACCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCca                                                             >  3:234/1‑148 (MQ=255)
                                                                                           ttCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACaa                                                          >  5:184360/1‑149 (MQ=255)
                                                                                           ttCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACaa                                                          <  7:32720/149‑1 (MQ=255)
                                                                                                cTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCAGTTCCC‑CACACAGCCACAACTTCt                                                     >  4:187191/1‑149 (MQ=255)
                                                                                                   cAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTgcg                                                  >  6:138169/1‑149 (MQ=255)
                                                                                                   cAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTgcg                                                  >  3:207195/1‑149 (MQ=255)
                                                                                                       gATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCAc                                              <  4:207195/149‑1 (MQ=255)
                                                                                                          gAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCtt                                                                                <  6:321678/113‑1 (MQ=255)
                                                                                                          gAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCtt                                                                                >  5:321678/1‑113 (MQ=255)
                                                                                                                  ggCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACcgcg                                   >  3:224368/1‑149 (MQ=255)
                                                                                                                       gCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTgcgc                              >  1:82480/1‑149 (MQ=255)
                                                                                                                         cTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCt                            >  5:102547/1‑149 (MQ=255)
                                                                                                                         cTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCAGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCt                            >  7:125179/1‑149 (MQ=255)
                                                                                                                            gCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCa                          <  2:235843/148‑1 (MQ=255)
                                                                                                                             ccGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCAcc                        <  4:234/149‑1 (MQ=255)
                                                                                                                                   aGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCACCGTCGaa                  <  4:267740/149‑1 (MQ=255)
                                                                                                                                     ttCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGCACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCACCGTCGAAgg                >  6:301836/1‑149 (MQ=255)
                                                                                                                                      tCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCACCGTCGAAGgc               >  5:320112/1‑149 (MQ=255)
                                                                                                                                            tGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTAATGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTACC‑CACACAGCCAAAACTTCTGCGCCACCTGCAACCGCGTGCTCCTCACCGTCCAACGCCGCctg         >  4:164484/1‑149 (MQ=255)
                                                                                                                                              gTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCCACACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCACCGTCGAAGGCCGCctgc        <  1:261025/149‑1 (MQ=255)
                                                                                                                                                 gAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCACCGTCGAAGGCCGCctgctgct    >  1:212361/1‑149 (MQ=255)
                                                                                                                                                  aaTCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCACCGTCGAAGGCCGCctgctgctg   >  8:286956/1‑149 (MQ=255)
                                                                                                                                                   aTCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACAGCCGGGTCGGCTTCATTTCCC‑CACAGAGCCACAACTTCTGCGCCTCCTGCAACCGCGTGCGCCTCACCGTCGAAGGCCGCCTGCTGCtgt  <  3:164484/149‑1 (MQ=255)
                                                                                                                                                   aTCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCACCGTCGAAGGCCGCCTGCTGCtgt  >  2:82073/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTGCGCTTTGCCATCGAACGCGAGCTGGACATCACCTTCATCGAGGAAATGCCGCTGGGGGTGATAAGCGAGCATGAACGTGGCGAATCGTTCTGCTCCAGCGATGAAGTACGGGCACGCCTGGCCGAGCAGTTCACGCTGGTCGAATCGACCGAGTCGTCCATGGGGCCGGCCCGCTACTGGCGCCTGGCCGAAGCCGCGAACACCCGGGTCGGCTTCATTTCCC‑CACACAGCCACAACTTCTGCGCCACCTGCAACCGCGTGCGCCTCACCGTCGAAGGCCGCCTGCTGCTGT  >  NC_002947/5217461‑5217756

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