Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,149,502 A→T 7.0% N237Y (AAC→TAC)  PP_2765 → sulfonate monooxygenase MsuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,149,5020AT 7.0% 157.7 / 2.8 57N237Y (AAC→TAC) PP_2765sulfonate monooxygenase MsuD
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base T (2/2);  total (29/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 = 8.07e-01

CTCGCGGGCCGCGCGCGACATGGCCTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGAAGTGAAGAACGCCGGTGCGGCCAGCCCGGAAGGCGAGGGCAA  >  NC_002947/3149360‑3149647
                                                                                                                                              |                                                                                                                                                 
cTCGCGGGCCGCGCGCGACATGGCCTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCtt                                                                                                                                             <  4:297608/149‑1 (MQ=255)
cTCGCGGGCCGCGCGCGACATGGCCTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCtt                                                                                                                                             <  2:94053/149‑1 (MQ=255)
   gcgGGCCGCGCGCGACATGGCCTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTtc                                                                                                                                           <  8:92348/148‑1 (MQ=255)
         cgcgcgcgACATGGACTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGACGGGGTAAACCATTTCATCATCg                                                                                                                                    >  7:169751/1‑149 (MQ=255)
          gcgcgcgACATGGCCTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGGGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGc                                                                                                                                   <  6:218294/149‑1 (MQ=255)
            gcgcgACATGGCCTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGccc                                                                                                                                 >  4:80380/1‑149 (MQ=255)
                     ggCCTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAATGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACAcc                                                                                                                         >  6:249814/1‑148 (MQ=255)
                     ggCCTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACAcc                                                                                                                         >  1:114190/1‑148 (MQ=255)
                      gCCTAGATGGTCTCCGACTGGTACTTCACCAATGCCAACAGTGTCGAGGGCATCAAGACCCAGCTCGACTATATCCGCACCAAGGCACCGGCAAACGTGCATTCGGTGAAGGTCGGGGTAAACTCTTTCATATTCTCCCGCTACACCGa                                                                                                                       <  8:169751/149‑1 (MQ=255)
                           gAAGGTCTCCGACTTGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAgg                                                                                                                     >  5:97934/1‑146 (MQ=255)
                           gAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGaa                                                                                                                   <  6:151202/148‑1 (MQ=255)
                             aGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAagag                                                                                                                <  7:95900/149‑1 (MQ=255)
                                          gTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCAGCGCGGTGc                                                                                                   >  7:132007/1‑149 (MQ=255)
                                           tACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCt                                                                                                  <  2:218712/149‑1 (MQ=255)
                                                   cAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCCAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCaga                                                                                            >  4:248575/1‑147 (MQ=255)
                                                   cAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAgaga                                                                                          <  1:292422/149‑1 (MQ=255)
                                                   cAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTAGGGGTAAACGCATTCATCATAGCCCGCGACACCGAGGAAGAGGCCCGAGCGGTTCTGGCa                                                                                              >  2:133030/1‑145 (MQ=255)
                                                      tGGCAACAGTGTCGAGGGAATCAAGACCCAGGTCGACGATATCCGCACAAAGTCAGCGGCCAACGGGCATTCGGTGAAGCTCGCGGTAAACGCTTTAATCATCGCCCGCGACACCGAGGAAGCAGCCCGCGCGGGGCTGGCAGAGatca                                                                                       >  4:63559/1‑149 (MQ=255)
                                                          aaCAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTATACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGc                                                                                   >  6:302290/1‑149 (MQ=255)
                                                          aaCAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTATACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGc                                                                                   >  6:297143/1‑149 (MQ=255)
                                                             aGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCaa                                                                                <  2:114190/149‑1 (MQ=255)
                                                                 tCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGcc                                                                            >  3:134724/1‑149 (MQ=255)
                                                                           cAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAgg                                                                  >  6:271852/1‑149 (MQ=255)
                                                                            aaGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACGCCGAGGAAGAGGCCCGCGCGGTGCTGACAGAGATCATCGCCAAGGCCGACCCGGAGGc                                                                 >  2:29171/1‑149 (MQ=255)
                                                                            aaGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGc                                                                 >  6:20198/1‑149 (MQ=255)
                                                                             aGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCg                                                                <  8:132007/149‑1 (MQ=255)
                                                                                 ccAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTc                                                             >  1:251568/1‑148 (MQ=255)
                                                                                   aGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTATACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAAc                                                          <  5:297143/149‑1 (MQ=255)
                                                                                   aGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTATACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAAc                                                          <  5:302290/149‑1 (MQ=255)
                                                                                   aGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCAACGCCAAGGCCGACCCGGAGGCGGTCAAc                                                          <  1:208157/149‑1 (MQ=255)
                                                                                      tCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGc                                                       >  6:143611/1‑149 (MQ=255)
                                                                                      tCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGATAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACggg                                                       >  8:346517/1‑148 (MQ=255)
                                                                                      tCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGATAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGc                                                       >  2:104546/1‑149 (MQ=255)
                                                                                        gacgaTATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCtt                                                     <  5:271852/149‑1 (MQ=255)
                                                                                        gacgaTATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCtt                                                     >  5:363480/1‑149 (MQ=255)
                                                                                         tgcttatCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGTGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCTACCCGGAGGCGGTCAACGGCttt                                                    <  1:133030/145‑1 (MQ=255)
                                                                                         acgaTATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCttt                                                    >  2:289523/1‑149 (MQ=255)
                                                                                         acgaTATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCttt                                                    <  6:97934/149‑1 (MQ=255)
                                                                                                 cGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCtt                                                     >  6:161850/1‑140 (MQ=255)
                                                                                                 cGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCtt                                                     <  5:161850/140‑1 (MQ=255)
                                                                                                  gCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCagtgaa                                           <  6:363480/149‑1 (MQ=255)
                                                                                                  gCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCagtg                                             >  2:336229/1‑147 (MQ=255)
                                                                                                     ccAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCagtgaagtg                                        <  1:289523/149‑1 (MQ=255)
                                                                                                     ccAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCagtgaagtg                                        <  4:134724/149‑1 (MQ=255)
                                                                                                      cgaGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCagtgaagtga                                       <  4:194148/147‑1 (MQ=255)
                                                                                                            aGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGa                                                                                                                       <  6:220958/63‑1 (MQ=255)
                                                                                                            aGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGa                                                                                                                       >  5:220958/1‑63 (MQ=255)
                                                                                                               ggCCAACGGGCATTCGGGGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGAAGTGAAGAACGCCg                              <  3:248575/149‑1 (MQ=255)
                                                                                                                      gggCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCagtgaagt                                         <  4:292781/131‑1 (MQ=255)
                                                                                                                      gggCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCagtgaagt                                         >  3:292781/1‑131 (MQ=255)
                                                                                                                             cGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGAAGTGAAGAACGCCGGTGCGGCCAGccc                 <  3:63559/148‑1 (MQ=255)
                                                                                                                               gTGAAGGTCGGGGTAAACGCTTTCATCCTCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGAAGTGCAGAACGCCGGTGCGGCCAGCCCgg               >  4:226027/1‑148 (MQ=255)
                                                                                                                                      tCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCAGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGAAGTGAAGAACGCCGGTGCGGCCAGCCCGGAAGGCGAg       >  6:167594/1‑149 (MQ=255)
                                                                                                                                       cGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGAAGTGAAGAACGCCGGTGCGGCCAGCCCGGAAGGCGAgg      <  1:336229/149‑1 (MQ=255)
                                                                                                                                        ggggTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGAAGTGAAGAACGCCGGTGCGGCCAGCCCGGAAGGCGAggg     <  5:20198/149‑1 (MQ=255)
                                                                                                                                          ggTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGAAGTGAAGAACGCCGGTGCGGCCAGCCCGGAAGGCGAGGGc    >  1:318682/1‑148 (MQ=255)
                                                                                                                                           gTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGATGTGAAGAACGCCGGTGCGGCCAGCCCGGAAGGCGAGGGCaa  >  6:366338/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
CTCGCGGGCCGCGCGCGACATGGCCTCGAAGGTCTCCGACTGGTACTTCACCAATGGCAACAGTGTCGAGGGCATCAAGACCCAGGTCGACGATATCCGCACCAAGGCAGCGGCCAACGGGCATTCGGTGAAGGTCGGGGTAAACGCTTTCATCATCGCCCGCGACACCGAGGAAGAGGCCCGCGCGGTGCTGGCAGAGATCATCGCCAAGGCCGACCCGGAGGCGGTCAACGGCTTTGGCAGTGAAGTGAAGAACGCCGGTGCGGCCAGCCCGGAAGGCGAGGGCAA  >  NC_002947/3149360‑3149647

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