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

GCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCAC  >  NC_000913/1609463‑1609659
                                                                                                   |                                                                                                  
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:1848369/100‑1 (MQ=255)
             ccacGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCggg                                                                                      >  1:1003450/1‑100 (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                                                        >  1:1902779/1‑100 (MQ=255)
                                           gCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGa                                                        <  2:1003439/100‑1 (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:560002/4‑100 (MQ=255)
                                                                                  catcctgccACAGCTGCTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAAc                 >  1:1760799/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)
                                                                                                 tgtgGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCAc  <  1:1552548/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GCCGTTACGCTCACCACGATAGAACGCAATCAATGCCGCTAATGCGAAAGTAAGGCGCGCCGGAAGTGTGCCGTTTGCCTTCTGCCCTGCCAGCAGCTGTGGCAGGATGCGGGTGCGGAACTTGGTCATACCGTTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCAC  >  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: