Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,609,555 Δ1 bp 18.2% coding (1126/1452 nt) uxaB ← altronate oxidoreductase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,609,5550G.18.2% 118.2 / 26.5 44coding (1126/1452 nt)uxaBaltronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (20/16);  new base . (5/3);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

GTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAC  >  NC_000913/1609458‑1609654
                                                                                                 |                                                                                                    
ntttCGCCGTTACGCTCACAACGATAGAACGCAATCAATGCCGCTACTGCGAAAGCAAGACGCAACGGAAGAGTGCCGTTTGCCTTCTGCCCTGCcagca                                                                                                    >  2:1142509/2‑100 (MQ=255)
gTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCcagca                                                                                                    >  1:836816/1‑100 (MQ=255)
   tCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCt                                                                                                 >  2:730333/1‑100 (MQ=255)
     gCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAACGTAGGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTCCCCTCCCAGCCGCtgt                                                                                              >  2:779499/1‑100 (MQ=255)
      ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtgtg                                                                                             <  2:28403/100‑1 (MQ=255)
       cGTTACGCTCACCATGATAGAAAGCAATCAATGAAGCTAATGCGAAACTCAGACGAGCCGGAAGTGAGCCGTTTGCATTCCGCTCTGCCCGCAGCtgtcg                                                                                            >  2:866541/1‑98 (MQ=255)
       cGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTgg                                                                                            <  1:1502432/100‑1 (MQ=255)
            cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGa                                                                                       >  2:1506786/1‑100 (MQ=255)
            cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAggg                                                                                       <  1:1056748/100‑2 (MQ=255)
            cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAggg                                                                                       <  1:597941/100‑2 (MQ=255)
            cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCACAGCTGCTGGCAggg                                                                                       <  1:514161/100‑2 (MQ=255)
                  ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCggg                                                                                 >  1:1003450/1‑100 (MQ=255)
                  ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCggg                                                                                 <  1:1848369/100‑1 (MQ=255)
                     cGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGc                                                                              >  2:1737855/1‑100 (MQ=255)
                      gATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCg                                                                             <  2:2275616/100‑1 (MQ=255)
                        tAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCgg                                                                            <  1:984363/99‑1 (MQ=255)
                          gAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAAc                                                                         >  2:147920/1‑100 (MQ=255)
                             cGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTg                                                                      >  1:400037/1‑100 (MQ=255)
                              gCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTgg                                                                     <  2:1521484/100‑1 (MQ=255)
                               caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGt                                                                    >  1:244683/1‑100 (MQ=255)
                                 atcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCa                                                                  <  1:304308/100‑1 (MQ=255)
                                              aTGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcgc                                                     <  1:1120320/100‑1 (MQ=255)
                                                gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGa                                                   <  2:1003439/100‑1 (MQ=255)
                                                gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGa                                                   >  1:1902779/1‑100 (MQ=255)
                                                 cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                    >  2:2051335/1‑98 (MQ=255)
                                                   aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                >  1:2387506/1‑100 (MQ=255)
                                                   aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                >  2:2130983/1‑100 (MQ=255)
                                                     aGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAc                                              >  1:996868/1‑100 (MQ=255)
                                                      gTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                             <  2:2015791/100‑1 (MQ=255)
                                                             gcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTg                                      <  2:2178079/100‑1 (MQ=255)
                                                              cgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGa                                     >  1:1088199/1‑100 (MQ=255)
                                                               gcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGTATGCGGGTGCGGAACTTGGTCATACCGTGGAGCGCGATAGACAGCAGCTGAt                                    <  2:1094441/100‑1 (MQ=255)
                                                                  ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGATCGCGATAGACAGCAGCTGATGCt                                 <  2:1350103/100‑1 (MQ=255)
                                                                  ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCt                                 <  1:454800/100‑1 (MQ=255)
                                                                     aGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                              >  1:1745534/1‑100 (MQ=255)
                                                                               ttGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGtt                    >  2:429640/1‑100 (MQ=255)
                                                                                    ttCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGa               <  1:1813522/100‑1 (MQ=255)
                                                                                      gcatcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGTTCATACCGTTGAGCGCGATAGACAGAAGCTGATGCTTACTGTACGGGTTACGGaaa             >  2:2284855/5‑100 (MQ=255)
                                                                                      ctgccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACAGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaaa             >  2:145690/1‑100 (MQ=255)
                                                                                       catcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:1760799/4‑100 (MQ=255)
                                                                                       catcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:560002/4‑100 (MQ=255)
                                                                                        gccctgccAGCAGCTGTGGCAGGATGCTGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg           >  2:175001/1‑100 (MQ=255)
                                                                                        gccctgccAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg           >  2:1573720/1‑100 (MQ=255)
                                                                                        atcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGATACGGaaa             >  1:437576/3‑98 (MQ=255)
                                                                                          cctgccACAGCTGCTGGCAGGATGCGGGTGAGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTGCGGAAACGcc         >  2:1574729/1‑100 (MQ=255)
                                                                                              ccAGCAGCTGTGGCAGGATGCGGGTCCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATTTATGGGTTACGGAAACGCCCGGt     <  2:157918/100‑1 (MQ=255)
                                                                                                 gcagcTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGCGAc  >  2:443565/1‑100 (MQ=255)
                                                                                                 |                                                                                                    
GTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAC  >  NC_000913/1609458‑1609654

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: