Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 835,596 Δ1 bp coding (349/963 nt) ybiB → putative family 3 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913835,5930G.97.5% 156.7 / ‑2.5 40coding (346/963 nt)ybiBputative family 3 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (20/19);  minor base C (1/0);  total (21/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 = 1.00e+00

GCTGACGCCGCCAGCAGGCAAGCCGATGCCGATTGTCATCCCCAGTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTCGACGAACATCAACCGGTGTTTATGCCAGTCGGCGC  >  NC_000913/835463‑835726
                                                                                                                                  |                                                                                                                                     
gCTGACGCCGCCAGCAGGCAAGCCGATGCCGATTGTCATCCCCAGTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCg                                                                                                                              <  1:544535/139‑1 (MQ=255)
           cagcagGCAAGCCGATGCCGATTGTCATCCCCAGTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAAcc                                                                                                                   <  2:29200/139‑1 (MQ=255)
                  cAAGCCGATGCCGATTGTCATCCCCAGTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGc                                                                                                            <  2:448869/139‑1 (MQ=255)
                         atgccgatTGTCATCCCCAGTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGctgactga                                                                                                     >  1:1437782/1‑139 (MQ=255)
                                     aTCCCCAGTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCCCGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGa                                                                                         >  2:218057/1‑139 (MQ=255)
                                        cccAGTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCAAAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATATTTGAAtt                                                                                      >  2:958125/1‑139 (MQ=255)
                                           aGTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGAt                                                                                   <  1:218057/139‑1 (MQ=255)
                                            gTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATg                                                                                  >  2:588315/1‑139 (MQ=255)
                                             ttACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATgg                                                                                 >  1:714351/1‑139 (MQ=255)
                                                       gcgcgTAAATACGCCAACCTGACCCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGcc                                                                       >  1:1298185/1‑139 (MQ=255)
                                                       gcgcgTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAAc                                                                                                  >  1:751642/1‑112 (MQ=255)
                                                       gcgcgTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAAc                                                                                                  <  2:751642/112‑1 (MQ=255)
                                                       gcgcgTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGcc                                                                       >  2:1383292/1‑139 (MQ=255)
                                                        cgcgTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATg                                                                                  >  1:1598229/1‑127 (MQ=255)
                                                        cgcgTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATg                                                                                  <  2:1598229/127‑1 (MQ=255)
                                                        cgcgTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCa                                                                      >  2:463281/1‑139 (MQ=255)
                                                              aaCAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTg                                                                <  1:588315/139‑1 (MQ=255)
                                                                   gCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTg                                                                                                             >  1:1270086/1‑89 (MQ=255)
                                                                   gCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTg                                                                                                             <  2:1270086/89‑1 (MQ=255)
                                                                   gCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGct                                                                                                           <  1:516707/91‑1 (MQ=255)
                                                                   gCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGct                                                                                                           >  2:516707/1‑91 (MQ=255)
                                                                   gCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCAcgg                                                           <  2:1356526/139‑1 (MQ=255)
                                                                    ccAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACgggtt                                                                                                                                   >  2:761236/1‑65 (MQ=255)
                                                                    ccAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACgggtt                                                                                                                                   <  1:761236/67‑3 (MQ=255)
                                                                    ccAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCAcggc                                                          >  2:1121442/1‑139 (MQ=255)
                                                                      aaCCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCAcggcgg                                                        <  1:1121442/139‑1 (MQ=255)
                                                                                  ttgttgGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGgcaggcgcaggc                                            >  2:326711/1‑139 (MQ=255)
                                                                                            ttCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTcgac                                  >  2:76093/1‑139 (MQ=255)
                                                                                             tctcCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTcgacg                                 <  1:1383292/139‑1 (MQ=255)
                                                                                                        aaCTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTCGACGAACATCAACCg                      <  1:958125/139‑1 (MQ=255)
                                                                                                        aaCTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTCGACGAACATCAACCg                      <  1:326711/139‑1 (MQ=255)
                                                                                                         aCTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGaa                                                                                        >  1:844773/1‑72 (MQ=255)
                                                                                                         aCTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGaa                                                                                        <  2:844773/72‑1 (MQ=255)
                                                                                                         aCTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGgcaggcg                                                 >  1:78832/1‑111 (MQ=255)
                                                                                                         aCTCGGTTTTCCTGTGGTGGTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGgcaggcg                                                 <  2:78832/111‑1 (MQ=255)
                                                                                                                       tggtggTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTAt                                                                             >  2:1237220/1‑69 (MQ=255)
                                                                                                                       tggtggTTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTAt                                                                             <  1:1237220/69‑1 (MQ=255)
                                                                                                                            gtTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTCGACGAACATCAACCg                      <  2:1480466/119‑1 (MQ=255)
                                                                                                                            gtTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTCGACGAACATCAACCg                      >  1:1480466/1‑119 (MQ=255)
                                                                                                                            gtTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTCGACGAACATCAACCGGTGTTTATGCCAGTCgg     <  2:1432012/136‑1 (MQ=255)
                                                                                                                            gtTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTCGACGAACATCAACCGGTGTTTATGCCAGTCgg     >  1:1432012/1‑136 (MQ=255)
                                                                                                                            gtTCACGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTCGACGAACATCAACCGGTGTTTATGCCAGTCGgcgc  >  1:757381/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                     
GCTGACGCCGCCAGCAGGCAAGCCGATGCCGATTGTCATCCCCAGTTACAACGGCGCGCGTAAACAGGCCAACCTGACGCCGTTGTTGGCGATTCTCCTGCACAAACTCGGTTTTCCTGTGGTGGTTCACGGGGTTAGCGAAGATCCAACCCGCGTGCTGACTGAAACCATTTTTGAATTGATGGGTATTACGCCAACGTTGCACGGCGGGCAGGCGCAGGCGAAGCTCGACGAACATCAACCGGTGTTTATGCCAGTCGGCGC  >  NC_000913/835463‑835726

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: