Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,843 C→G 12.0% intergenic (‑279/‑43) nicC ← / → nicX 6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,447,8430CG12.0% 136.7 / 9.1 50intergenic (‑279/‑43)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (19/24);  new base G (3/3);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

ATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCAC  >  NC_002947/4447695‑4447984
                                                                                                                                                      |                                                                                                                                             
aTCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACAGAGATGCAGCACCAGGCAc                                                                                                                                               >  3:290597/1‑149 (MQ=255)
             ttACCGTAAATACTGGTCAATATGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                   <  1:136805/148‑2 (MQ=255)
                 cGTAAATACTTGTCAATGTGTGGATAGGTAAAGTTTTAACAGATTACTTGCGCTCTCTGGCAGGATGCGGTTTCATCTGCACAATTATGGTGTACAAACTATTATTTTGCACCGAGATTCAGCACCAGGCACTGCCATGGTGCGCCCCt                                                                                                                              <  2:133086/149‑1 (MQ=255)
                 cGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCt                                                                                                                              >  6:150459/1‑149 (MQ=255)
                  gTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTg                                                                                                                             >  7:246665/1‑149 (MQ=255)
                   tAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGt                                                                                                                            <  6:59030/149‑1 (MQ=255)
                     aaTACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCa                                                                                                                          <  4:248475/149‑1 (MQ=255)
                         cTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGc                                                                                                                      >  2:47480/1‑149 (MQ=255)
                          tGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                   >  1:232979/1‑134 (MQ=255)
                          tGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                   <  2:232979/135‑2 (MQ=255)
                            gTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCAc                                                                                                                   >  8:131092/1‑149 (MQ=255)
                               aaTGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCAcc                                                                                                                  >  2:78935/1‑147 (MQ=255)
                                aTGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAg                                                                                                               <  5:65919/149‑1 (MQ=255)
                                     tgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGaa                                                                                                          >  4:84856/1‑149 (MQ=255)
                                       gATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAAc                                                                                                         <  3:84856/148‑1 (MQ=255)
                                            gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACAAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                   >  3:109199/1‑149 (MQ=255)
                                                cGTTTGAACAGTTGTCTTGCGCTCGCTGGGAGGTTGCGTTTTCATCTGCACATGTAGGGTGGACACACTATTATTTTGCAGCGAGATGGAGCACCAGGCACTGCCATGGTGCGCCCCTGTCATCGCCACCCAGGAGAACCCCGATGcgg                                                                                               <  4:109199/149‑3 (MQ=255)
                                                     gAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAg                                                                                          >  2:296748/1‑149 (MQ=255)
                                                      aaCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCAGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAg                                                                                          <  7:242688/148‑1 (MQ=255)
                                                          gATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATg                                                                                     <  1:78935/149‑1 (MQ=255)
                                                                ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTAATGCACCGAGATGAAGCACAAGGCCATGCAATGGTGCCCCCCTGGAAACGCCCCCCAGGAAAACCCCGATGCCGGTGGGCATTGcaca                                                                                 >  1:87996/1‑148 (MQ=255)
                                                                  gcgcTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACt                                                                              <  2:274143/148‑1 (MQ=255)
                                                                        gctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTTTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCa                                                                        <  5:290619/148‑1 (MQ=255)
                                                                        gctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAg                                                                       <  5:150459/149‑1 (MQ=255)
                                                                              aGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGAcc                                                                          >  4:11900/1‑140 (MQ=255)
                                                                              aGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGAcc                                                                          <  3:11900/140‑1 (MQ=255)
                                                                                  tGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGc                                                             <  2:87996/149‑1 (MQ=255)
                                                                                    cGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCAc                                                           <  1:296748/149‑1 (MQ=255)
                                                                                     gTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACg                                                          <  8:246665/149‑1 (MQ=255)
                                                                                        ttCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTa                                                        <  4:290597/148‑1 (MQ=255)
                                                                                         tCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTAc                                                       <  1:47480/148‑1 (MQ=255)
                                                                                                   aaGTAGGGTGTACACACTATTATTTTGCACCGAGCTGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGtc                                             >  5:18789/1‑148 (MQ=255)
                                                                                                     gTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGc                                                             <  4:8877/130‑1 (MQ=255)
                                                                                                     gTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGc                                                             >  3:8877/1‑130 (MQ=255)
                                                                                                     gTAGGGTGTAAACACTATTATATTGCACCGAGATGCAGCTCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAATCTGtctgg                                          <  6:226542/149‑3 (MQ=255)
                                                                                                         ggTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                       <  4:69336/148‑1 (MQ=255)
                                                                                                          gtgtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCg                                                                                               >  7:187447/1‑91 (MQ=255)
                                                                                                          gtgtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCg                                                                                               <  8:187447/91‑1 (MQ=255)
                                                                                                                acacTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGaaa                               >  1:283024/1‑149 (MQ=255)
                                                                                                                     attattTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACg                             >  8:202244/1‑146 (MQ=255)
                                                                                                                                aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTg                                      >  7:238429/1‑127 (MQ=255)
                                                                                                                                aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTg                                      <  8:238429/127‑1 (MQ=255)
                                                                                                                                 ccGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGt               <  7:131092/149‑1 (MQ=255)
                                                                                                                                        gcagcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGTTc           >  4:49890/1‑146 (MQ=255)
                                                                                                                                        gcagcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTc           >  6:75050/1‑146 (MQ=255)
                                                                                                                                         cagcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCaa         >  2:156837/1‑147 (MQ=255)
                                                                                                                                          cgcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCaa         >  7:274012/2‑146 (MQ=255)
                                                                                                                                          cgcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAaga       <  8:274012/147‑1 (MQ=255)
                                                                                                                                           gcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGc     >  8:35523/1‑149 (MQ=255)
                                                                                                                                             aCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCa   <  7:202244/149‑1 (MQ=255)
                                                                                                                                               cAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCac  >  4:264873/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
ATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCAC  >  NC_002947/4447695‑4447984

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: