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

GTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAC  >  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)
     gCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAACGTAGGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTCCCCTCCCAGCCGCtg‑t                                                                                              >  2:779499/1‑100 (MQ=255)
      ccGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCtg‑tg                                                                                             <  2:28403/100‑1 (MQ=255)
       cGTTACGCTCACCATGATAGAAAGCAATCAATGAAGCTAATGCGAAACTCAGACGAGCCGGAAGTGAGCCGTTTGCATTCCGCTCTGCCCGCAGCtg‑tcg                                                                                            >  2:866541/1‑98 (MQ=255)
       cGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑Tgg                                                                                            <  1:1502432/100‑1 (MQ=255)
            cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGa                                                                                       >  2:1506786/1‑100 (MQ=255)
            cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  1:1056748/100‑2 (MQ=255)
            cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  1:597941/100‑2 (MQ=255)
            cGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCA‑CAGCTGCTGGCAggg                                                                                       <  1:514161/100‑2 (MQ=255)
                  ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                 >  1:1003450/1‑100 (MQ=255)
                  ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                 <  1:1848369/100‑1 (MQ=255)
                     cGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGc                                                                              >  2:1737855/1‑100 (MQ=255)
                      gATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCg                                                                             <  2:2275616/100‑1 (MQ=255)
                        tAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCgg                                                                            <  1:984363/99‑1 (MQ=255)
                          gAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAAc                                                                         >  2:147920/1‑100 (MQ=255)
                             cGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTg                                                                      >  1:400037/1‑100 (MQ=255)
                              gCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTgg                                                                     <  2:1521484/100‑1 (MQ=255)
                               caatcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGt                                                                    >  1:244683/1‑100 (MQ=255)
                                 atcaatGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCa                                                                  <  1:304308/100‑1 (MQ=255)
                                              aTGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAgcgc                                                     <  1:1120320/100‑1 (MQ=255)
                                                gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGa                                                   <  2:1003439/100‑1 (MQ=255)
                                                gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGa                                                   >  1:1902779/1‑100 (MQ=255)
                                                 cGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGcgcg                                                    >  2:2051335/1‑98 (MQ=255)
                                                   aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                >  1:2387506/1‑100 (MQ=255)
                                                   aaaGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAg                                                >  2:2130983/1‑100 (MQ=255)
                                                     aGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAc                                              >  1:996868/1‑100 (MQ=255)
                                                      gTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGAca                                             <  2:2015791/100‑1 (MQ=255)
                                                             gcgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTg                                      <  2:2178079/100‑1 (MQ=255)
                                                              cgcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGa                                     >  1:1088199/1‑100 (MQ=255)
                                                               gcCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGTATGCGGGTGCGGAACTTGGTCATACCGTGGAGCGCGATAGACAGCAGCTGAt                                    <  2:1094441/100‑1 (MQ=255)
                                                                  ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGATCGCGATAGACAGCAGCTGATGCt                                 <  2:1350103/100‑1 (MQ=255)
                                                                  ggAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCt                                 <  1:454800/100‑1 (MQ=255)
                                                                     aGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTaa                              >  1:1745534/1‑100 (MQ=255)
                                                                               ttGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGtt                    >  2:429640/1‑100 (MQ=255)
                                                                                    ttCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGa               <  1:1813522/100‑1 (MQ=255)
                                                                                      gcatcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGTTCATACCGTTGAGCGCGATAGACAGAAGCTGATGCTTACTGTACGGGTTACGGaaa             >  2:2284855/5‑100 (MQ=255)
                                                                                      ctgccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACAGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGaaa             >  2:145690/1‑100 (MQ=255)
                                                                                       catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:1760799/4‑100 (MQ=255)
                                                                                       catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc            >  1:560002/4‑100 (MQ=255)
                                                                                        gccctgccAGCAGCTG‑TGGCAGGATGCTGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg           >  2:175001/1‑100 (MQ=255)
                                                                                        gccctgccAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACg           >  2:1573720/1‑100 (MQ=255)
                                                                                        atcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGATACGGaaa             >  1:437576/3‑98 (MQ=255)
                                                                                          cctgccA‑CAGCTGCTGGCAGGATGCGGGTGAGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTGCGGAAACGcc         >  2:1574729/1‑100 (MQ=255)
                                                                                              ccAGCAGCTG‑TGGCAGGATGCGGGTCCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATTTATGGGTTACGGAAACGCCCGGt     <  2:157918/100‑1 (MQ=255)
                                                                                                 gcagcTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGCGAc  >  2:443565/1‑100 (MQ=255)
                                                                                                 |                                                                                                    
GTTTCGCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGAC  >  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: ‑