Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,888,321 T→C 7.0% F109S (TTC→TCC)  PP_3431 → ThiJ/PfpI family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,888,3210TC 7.0% 164.2 / 2.4 58F109S (TTC→TCC) PP_3431ThiJ/PfpI family protein
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base C (2/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

AGATGGCCCCGGACCTGGACTGGCAGTGGTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAACTGCCTGCTGCTCAAAAACCGCAAGGTGAC  >  NC_002947/3888177‑3888465
                                                                                                                                                |                                                                                                                                                
aGATGGCCCCGGACCTGGACTGGCAGTGGTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCa                                                                                                                                              <  4:173290/149‑1 (MQ=255)
       cccGGACCTGGACTGGCAGTGGTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAAcc                                                                                                                                       <  6:367814/149‑1 (MQ=255)
         cGGACCTGGACTGGCAGTGGTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCAt                                                                                                                                     <  3:66146/149‑1 (MQ=255)
                 gACTGGCAGTGGTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTg                                                                                                                               <  3:17917/147‑1 (MQ=255)
                      gCAGTGGTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGAcc                                                                                                                        >  8:181792/1‑149 (MQ=255)
                         gTGGTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTgg                                                                                                                     <  5:141809/149‑1 (MQ=255)
                            gTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGAc                                                                                                                                           >  5:168660/1‑124 (MQ=255)
                            gTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGAc                                                                                                                                           <  6:168660/124‑1 (MQ=255)
                            gTACGACGACAAGCGCTTCATGACGCATCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCAc                                                                                                                  <  5:3597/149‑1 (MQ=255)
                                gacgacAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGGGGGCATGGGGTGATGTATGACTACGCAGACAACCATGCTCTGCAGGACAGGGCAGGct                                                                                                               >  2:347477/1‑147 (MQ=255)
                                      aaGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTa                                                                                                        <  7:55336/149‑1 (MQ=255)
                                      aaGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTa                                                                                                        <  3:53192/149‑1 (MQ=255)
                                       aGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTAc                                                                                                       >  4:43254/1‑149 (MQ=255)
                                            ttCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGaaa                                                                                                   >  4:39069/1‑148 (MQ=255)
                                              cATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAg                                                                                                >  5:233490/1‑149 (MQ=255)
                                                 gACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTg                                                                                             >  7:178206/1‑149 (MQ=255)
                                                 gACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTCCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTg                                                                                             >  8:3399/1‑149 (MQ=255)
                                                 gACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTCCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTg                                                                                             >  8:8792/1‑149 (MQ=255)
                                                            gCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTc                                                                                                                                                 <  7:154502/86‑1 (MQ=255)
                                                            gCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTc                                                                                                                                                 >  8:154502/1‑86 (MQ=255)
                                                            gCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCgca                                                                                  <  3:39069/149‑1 (MQ=255)
                                                               gTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCgcagca                                                                               <  3:148133/149‑1 (MQ=255)
                                                                tACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGcagcag                                                                              >  8:96079/1‑149 (MQ=255)
                                                                 aCCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGt                                                                             >  5:133998/1‑149 (MQ=255)
                                                                   cTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCt                                                                           >  7:21784/1‑149 (MQ=255)
                                                                      ggCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTCCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGt                                                                         <  7:3399/148‑1 (MQ=255)
                                                                      ggCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTCCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGt                                                                         <  7:8792/148‑1 (MQ=255)
                                                                         cccAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGt                                                                             >  3:235843/1‑141 (MQ=255)
                                                                         cccAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACg                                                                     <  7:181792/149‑1 (MQ=255)
                                                                         cccAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACg                                                                     <  3:43254/149‑1 (MQ=255)
                                                                            aGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGTGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGt                                                                                              <  7:210180/121‑1 (MQ=255)
                                                                            aGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTAATACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGt                                                                                              >  8:210180/1‑121 (MQ=255)
                                                                               gCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAg                                                               <  8:178206/149‑1 (MQ=255)
                                                                                       aGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCctg                                                       <  6:233490/149‑1 (MQ=255)
                                                                                          cTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCctgctg                                                    <  6:133998/149‑1 (MQ=255)
                                                                                           tGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCctgctg                                                    >  3:126442/1‑148 (MQ=255)
                                                                                             aCGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGa                                                   <  8:185192/147‑1 (MQ=255)
                                                                                             aCGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAAc                                                 <  7:96079/149‑1 (MQ=255)
                                                                                               gAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATTTATGACTTCGCAGACAACCATGCTCTGCAGGAGCTGGCACGCAAAATCTACGAAAAAGGTCGCAGTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACAt                                               <  5:301938/149‑1 (MQ=255)
                                                                                                   aCAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCt                                                                           >  8:24209/1‑117 (MQ=255)
                                                                                                   aCAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCt                                                                           <  7:24209/117‑1 (MQ=255)
                                                                                                    cAGCGTCATCTATTACACCGGTGGGCATTTGGTGATGTATGACTTCGCATACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGTTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAg                                           <  1:347477/148‑1 (MQ=255)
                                                                                                       cGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGa                                       <  4:126442/149‑1 (MQ=255)
                                                                                                       cGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGa                                       >  4:48895/1‑149 (MQ=255)
                                                                                                        gTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGa                                       <  1:165817/148‑1 (MQ=255)
                                                                                                            tCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTaa                                   >  2:373389/1‑148 (MQ=255)
                                                                                                            tCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTaa                                   >  4:250537/1‑148 (MQ=255)
                                                                                                            tCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTaa                                   >  1:355059/1‑148 (MQ=255)
                                                                                                                  acacCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAActgc                            >  8:188708/1‑149 (MQ=255)
                                                                                                                    acCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAActgcct                          >  3:187946/1‑149 (MQ=255)
                                                                                                                     ccGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAActgcctg                         >  1:141081/1‑149 (MQ=255)
                                                                                                                      cGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTAGGCCTGCTGAACATCAACCTGAGTAACAActgcctgc                        >  3:338653/1‑149 (MQ=255)
                                                                                                                                 gggTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAACTGCCTGCTGCTCAAAAAc             >  2:276786/1‑149 (MQ=255)
                                                                                                                                   gTGAAGTATGACTTCTCAGACAACCAAGATCATCAGGACGTGGCACTAAAAATCTACGAAAAAGGTGTCATTGTCGAAGCACTCATTCACGGAGTAGTCGGCATGCTGAACATAAAGAAGAGTACCAAATGACAGCTGCTCAAAAAccc           >  4:297847/1‑148 (MQ=255)
                                                                                                                                    tGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAACTGCCTGCTACTCAAAAACCGc          <  1:373389/149‑1 (MQ=255)
                                                                                                                                         tATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAACTGCCTGCTGCTCAAAAACCGCAAGGt     >  4:346587/1‑149 (MQ=255)
                                                                                                                                            gACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAACTGCCTGCTGCTCAAAAACCGCAAGGTGAc  >  2:34786/1‑149 (MQ=255)
                                                                                                                                            gACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACCGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAACTGCCTGCTGCTCAAAAACCGCAAGGTgcc  >  6:83225/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AGATGGCCCCGGACCTGGACTGGCAGTGGTACGACGACAAGCGCTTCATGACGCGTCTTGGCAGTACCTTGGCCCCAGGGCAGGTGAAGGCTGACGAATACAGCGTCATCTATTACACCGGTGGGCATGGGGTGATGTATGACTTCGCAGACAACCATGCTCTGCAGGACCTGGCACGCAAAATCTACGAAAAAGGTGGCATTGTCGCAGCAGTCTGTCACGGCGTAGTCGGCCTGCTGAACATCAAGCTGAGTAACAACTGCCTGCTGCTCAAAAACCGCAAGGTGAC  >  NC_002947/3888177‑3888465

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