Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,609,561:1 +C 19.5% coding (1120/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,5611.C19.5% 112.9 / 32.9 41coding (1120/1452 nt)uxaBaltronate oxidoreductase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base . (18/15);  new base C (5/3);  total (23/18)
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.83e-01

GCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCAC  >  NC_000913/1609463‑1609659
                                                                                                   |                                                                                                  
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:1848369/100‑1 (MQ=255)
             ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCggg                                                                                      >  1:1003450/1‑100 (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                                                        >  1:1902779/1‑100 (MQ=255)
                                           gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGa                                                        <  2:1003439/100‑1 (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:560002/4‑100 (MQ=255)
                                                                                  catcctgccA‑CAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc                 >  1:1760799/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)
                                                                                                 tg‑tgGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCAc  <  1:1552548/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTG‑TGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCAC  >  NC_000913/1609463‑1609659

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