Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 292,606 C→G 7.6% intergenic (+22/‑38) ssuA → / → ssuD aliphatic sulfonatesABC transporter substrate‑binding protein/alkanesulfonate monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947292,6060CG 7.6% 150.5 / 2.0 49intergenic (+22/‑38)ssuA/ssuDaliphatic sulfonatesABC transporter substrate‑binding protein/alkanesulfonate monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base G (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCAG  >  NC_002947/292459‑292754
                                                                                                                                                   |                                                                                                                                                    
ccAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGcgcc                                                                                                                                                     >  8:267460/1‑149 (MQ=255)
ccAGAGGTGGTGGCCGCGCAAAAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCATTGATTGACTGACAGTGACGGGGcgcc                                                                                                                                                     >  8:196164/1‑149 (MQ=255)
       tggtggCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCgg                                                                                                                                              <  8:55448/149‑1 (MQ=255)
          tggCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGGCGCCCCGGcac                                                                                                                                           <  2:197832/149‑3 (MQ=255)
            gCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTaga                                                                                                                                         <  8:12479/149‑1 (MQ=255)
                      agaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCAccc                                                                                                                                  >  7:212471/1‑146 (MQ=255)
                      agaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCAccc                                                                                                                                  <  8:212471/146‑1 (MQ=255)
                       gaagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTga                                                                                                                              <  1:235711/149‑1 (MQ=255)
                        aagaTCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                             <  7:267460/149‑1 (MQ=255)
                                  acacCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                             <  7:216758/139‑1 (MQ=255)
                                  acacCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                             >  8:216758/1‑139 (MQ=255)
                                   cacCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTc                                                                                                                   >  1:227719/1‑148 (MQ=255)
                                      cTTCCAGGCGCTCAAGCTGATAACCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAAt                                                                                                               <  6:267092/149‑1 (MQ=255)
                                            ggCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCt                                                                                                         >  2:243856/1‑149 (MQ=255)
                                                   aaGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCTTTGAGGAGACGACTCCAATGAGCCTTAACAt                                                                                                   <  1:31147/148‑1 (MQ=255)
                                                    aGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACAttt                                                                                                 >  7:43188/1‑149 (MQ=255)
                                                       tGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCt                                                                                              >  3:154931/1‑149 (MQ=255)
                                                          tACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGt                                                                                           <  4:154931/149‑1 (MQ=255)
                                                            cccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTaa                                                                                                      >  8:110262/1‑136 (MQ=255)
                                                            cccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTaa                                                                                                      <  7:110262/136‑1 (MQ=255)
                                                             ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtcc                                                                                        >  2:178391/1‑149 (MQ=255)
                                                             ccAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTCATTGACTGCCAGTGACGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTtc                                                                                         >  7:200464/1‑148 (MQ=255)
                                                              cAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                             <  2:124827/111‑1 (MQ=255)
                                                              cAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTgag                                                                                                                             >  1:124827/1‑111 (MQ=255)
                                                                         gCATCAAGGACGTTATGTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTTCCGGTGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACg                                                                            <  7:196164/149‑1 (MQ=255)
                                                                         gCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACg                                                                            <  5:56241/149‑1 (MQ=255)
                                                                         gCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACg                                                                            >  6:56241/1‑149 (MQ=255)
                                                                           aTCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGc                                                                                                           <  6:13855/116‑1 (MQ=255)
                                                                           aTCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGt                                                                          >  8:123721/1‑149 (MQ=255)
                                                                            tCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGTCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGt                                                                          >  2:158550/1‑148 (MQ=255)
                                                                            tCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTg                                                                         >  8:56389/1‑149 (MQ=255)
                                                                              aaGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGa                                                                        <  1:178391/148‑1 (MQ=255)
                                                                                         cTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCt                                                            >  4:234494/1‑149 (MQ=255)
                                                                                              cacCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCa                                                       <  2:227719/149‑1 (MQ=255)
                                                                                                 ccccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACtt                                                    <  8:200464/149‑1 (MQ=255)
                                                                                                  cccTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTc                                                   >  1:155720/1‑149 (MQ=255)
                                                                                                       ccAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGaa                                               <  2:155720/148‑1 (MQ=255)
                                                                                                         aaaGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGgc                                            >  4:200936/1‑149 (MQ=255)
                                                                                                                   gcgcgcCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCgg                                  <  1:243856/149‑1 (MQ=255)
                                                                                                                    cgcgcCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGt                                 <  7:56389/149‑1 (MQ=255)
                                                                                                                         cTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGAcc                            >  6:281876/1‑149 (MQ=255)
                                                                                                                              ttgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGgcgc                                          >  4:248623/1‑130 (MQ=255)
                                                                                                                              ttgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGgcgc                                          <  3:248623/130‑1 (MQ=255)
                                                                                                                               tgaCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGt                       >  2:264775/1‑148 (MQ=255)
                                                                                                                                 aCTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTAc                    <  3:234494/149‑1 (MQ=255)
                                                                                                                                   tGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCt                  >  5:253650/1‑149 (MQ=255)
                                                                                                                                     ccAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTg                 >  6:245166/1‑148 (MQ=255)
                                                                                                                                      tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCccacccac                                                                             >  3:6569/5‑87 (MQ=255)
                                                                                                                                      tctagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCccacccac                                                                             <  4:6569/83‑1 (MQ=255)
                                                                                                                                        tagCCGGGGCGGCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTgcag              >  2:104110/3‑148 (MQ=255)
                                                                                                                                               ggCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATcgc      <  6:253650/149‑1 (MQ=255)
                                                                                                                                                  gccgccCCGTCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCa   <  3:200936/149‑1 (MQ=255)
                                                                                                                                                  gccgccCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATcgcg     >  5:24035/1‑147 (MQ=255)
                                                                                                                                                   ccgccCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCAg  >  2:249603/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCAGAGGTGGTGGCCGCGCAACAGAAGATCGCCGACACCTTCCAGGCGCTCAAGCTGATACCCAAGCCGCTGAGCATCAAGGACGTGATCTGGACACCCCCTGCCAAAGTCGCTAGCGCGCCTTGATTGACTGCCAGTGCCGGGGCGCCGCCCCGGCTAGAGCCACCCTTGAGGAGACGACTCCAATGAGCCTTAACATTTTCTGGTTCCTCCCCACCCACGGTGACGGCAAGTACCTGGGCACTTCCGAAGGCGCACGTGCGGTCGACCACGGTTACCTGCAGCAGATCGCGCAG  >  NC_002947/292459‑292754

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: ‑