Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,446,852 T→C 8.0% Q238R (CAG→CGG)  nicC ← 6‑hydroxynicotinate 3‑monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,446,8520TC 8.0% 141.9 / 2.3 50Q238R (CAG→CGG) nicC6‑hydroxynicotinate 3‑monooxygenase
Reads supporting (aligned to +/- strand):  ref base T (20/26);  new base C (2/2);  total (22/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.97e-01

CACGGCTCCACAACGGCAACGGGTTGCGGTTGCGCAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCGAACACATCGGCATGCTGTGCCAGG  >  NC_002947/4446706‑4446992
                                                                                                                                                  |                                                                                                                                            
cACGGCTCCACAACGGCAACGGGTTGCGGTTGCGCAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTg                                                                                                                                             >  7:51391/1‑148 (MQ=255)
  cGGCTCAAAAACGGAAACGGGTTGGGGTTGCGAAGGGGCCTCCCCAAGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGaa                                                                                                                                          <  1:105483/149‑1 (MQ=255)
         acaACGGCAACGGGTTGCGGTTGCGCAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAg                                                                                                                                   >  6:47161/1‑149 (MQ=255)
         acaACGGCAACGGGTTGCGGTTGCGCAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAg                                                                                                                                   <  5:322522/149‑1 (MQ=255)
           aacggcaacggGTTGCGGTTGCGCAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCa                                                                                                                                    <  4:81655/146‑1 (MQ=255)
           aacggcaacggGTTGCGGTTGCGCAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCa                                                                                                                                    >  3:81655/1‑146 (MQ=255)
           aacggcaacggGTTGCGGTTGCGCAGCCGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGc                                                                                                                                 <  3:292264/149‑1 (MQ=255)
              ggcaacggGTTGCGGTTGCGCAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTc                                                                                                                              <  7:2551/149‑1 (MQ=255)
                  acggGTTGCGGTTGCGCAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCAtgt                                                                                                                          <  1:138121/149‑1 (MQ=255)
                                cgcAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGa                                                                                                             <  5:71699/148‑1 (MQ=255)
                                 gcAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGAc                                                                                                            >  7:126825/1‑148 (MQ=255)
                                      ggCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGt                                                                                                       >  8:16796/1‑148 (MQ=255)
                                      ggCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTa                                                                                                      <  8:126825/149‑1 (MQ=255)
                                         cACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATa                                                                                                   <  7:222073/149‑1 (MQ=255)
                                         cACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATa                                                                                                   <  5:47161/149‑1 (MQ=255)
                                          aCTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTc                                                                                                                              >  5:129666/1‑121 (MQ=255)
                                          aCTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTc                                                                                                                              <  6:129666/121‑1 (MQ=255)
                                          aCTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTaa                                                                                                     >  1:25245/1‑146 (MQ=255)
                                          aCTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTaa                                                                                                     <  2:25245/146‑1 (MQ=255)
                                          aCTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATAc                                                                                                  >  1:54302/1‑149 (MQ=255)
                                           cTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCATGCCGGTGACGAAGTAATACt                                                                                                 >  4:298101/1‑149 (MQ=255)
                                              ggTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACtc                                                                                                <  2:301187/147‑1 (MQ=255)
                                                 gATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATa                                                                                                   <  3:5427/141‑1 (MQ=255)
                                                 gATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATa                                                                                                   >  4:5427/1‑141 (MQ=255)
                                                    cGAGTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTTGCTGCTGTCGACGATCGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCgcgc                                                                                         <  5:138842/148‑1 (MQ=255)
                                                      aTTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTg                                                                                      >  5:264402/1‑149 (MQ=255)
                                                                   gATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCCGGAAGTcc                                                                                                                                      >  8:268580/1‑88 (MQ=255)
                                                                   gATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCCGGAAGTcc                                                                                                                                      >  8:272884/1‑88 (MQ=255)
                                                                   gATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCCGGAAGTcc                                                                                                                                      <  7:272884/88‑1 (MQ=255)
                                                                   gATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCCGGAAGTcc                                                                                                                                      <  7:268580/88‑1 (MQ=255)
                                                                    aTCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCATGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTgtagta                                                                        <  3:298101/149‑1 (MQ=255)
                                                                       aGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGa                                                                      <  3:201141/148‑1 (MQ=255)
                                                                              gCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCAtg                                                              >  7:195487/1‑149 (MQ=255)
                                                                                aCCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCggtggt                                                                              <  7:61525/131‑1 (MQ=255)
                                                                                aCCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCggtggt                                                                              >  8:61525/1‑131 (MQ=255)
                                                                                         tggtAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTc                                                   <  7:16796/149‑1 (MQ=255)
                                                                                                   cGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCAt                                         >  1:40832/1‑149 (MQ=255)
                                                                                                     aaGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCAttt                                       <  2:40832/149‑1 (MQ=255)
                                                                                                           gcgcgcATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCAtt                                        >  6:105478/1‑142 (MQ=255)
                                                                                                           gcgcgcATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCAtt                                        <  5:105478/142‑1 (MQ=255)
                                                                                                               gcATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGc                             >  7:232660/1‑149 (MQ=255)
                                                                                                                     cttcCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCGaa                        <  5:137232/148‑1 (MQ=255)
                                                                                                                     cttcCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCGAac                       >  4:94989/1‑149 (MQ=255)
                                                                                                                             ctgctgTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCGAACACATCGGc               <  2:54302/149‑1 (MQ=255)
                                                                                                                                  gTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCAATTGACGCAGGGCTCGAACACATCGGCATGCt          >  6:81021/1‑149 (MQ=255)
                                                                                                                                   tCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCGAACACATCGGCATGCtg         >  7:296425/1‑149 (MQ=255)
                                                                                                                                     gacgaACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCGAACACATCGGCATGCtgt        <  6:264402/148‑1 (MQ=255)
                                                                                                                                        gaACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCg                          >  7:273135/1‑127 (MQ=255)
                                                                                                                                        gaACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCg                          <  8:273135/127‑1 (MQ=255)
                                                                                                                                          aCGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCGAACACATCGGCATGCTGTGCCAgg  <  8:232660/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CACGGCTCCACAACGGCAACGGGTTGCGGTTGCGCAGCGGCCACTTGGTGATCGATTCGGTAGCGTCGATCAGCTTCTGCACCGTGGGGTGGTAGCCTTCGAAGGCTGCGCGCATTTCTTCCTGGCTGCTGTCGACGAACGCACCCTGGAAGTCCCAGGCTTCATGTGGCACGCCGGTGACGAAGTAATACTCGTCGCGCTTGCCGGTGGTGTAGTAGACCATCATGTGGCGGTCTTCGGACCACCATTTGACGCAGGGCTCGAACACATCGGCATGCTGTGCCAGG  >  NC_002947/4446706‑4446992

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