Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 658,184 C→T intergenic (‑246/‑763) crcB ← / → tatE F(‑) channel protein/twin arginine protein translocation system ‑ TatE protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913658,1840CT94.1% 42.6 / NA 17intergenic (‑246/‑763)crcB/tatEF(‑) channel protein/twin arginine protein translocation system ‑ TatE protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/6);  minor base G (0/1);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCG  >  NC_000913/657940‑658345
                                                                                                                                                                                                                                                    |                                                                                                                                                                 
gCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATcn                                                                                                                                                              <  1:150644/250‑2 (MQ=255)
                                                        nttGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTCCATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAg                                                                                                       >  1:53114/2‑249 (MQ=255)
                                                                                 ggaggTGGCAATGCTGGTTGCCGCAGNACAGTTNNCTNTTACCTCTGTGTGGGAAAANAANNNTNNNANTTGTGNCTNNTTGATGGCGNAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGNGCGATGCTCNNGATGCAGATTTATCGGTTAAANCAGNNCNncngct                                                                                                                                          >  2:56891/1‑189 (MQ=255)
                                                                                 ggaggTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACAGCTGTGTGGGAAAAGAACGCGGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCgg                                                                                                                        <  1:56891/207‑1 (MQ=255)
                                                                                                     ccGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAACACGACGCATCGATGTTTGCCCTGCCGGAAGCCTTGCTGGCGCGCGATGATCATGATGCAGCTTTATCGGTTAAATCAGCACAGCTGCTGGAAGGc                                                                                                                                   >  1:142029/1‑176 (MQ=255)
                                                                                                           gACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCgg                                                  <  1:11368/251‑1 (MQ=255)
                                                                                                                                     gAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGc                         >  1:83633/1‑250 (MQ=255)
                                                                                                                                                  agaTTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCNg                                                                                                                                                    <  1:198539/114‑1 (MQ=255)
                                                                                                                                                   naTTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGGAAACGTAACATGATGACGACAATCCTGACTATTCATGTTCCTTCAACGCCGGGGCGAGCAAGGAATATGCTGGAGGCACCACAAGCAGGaa           >  1:204798/2‑250 (MQ=255)
                                                                                                                                                           ccTCGTTGATGGCGCAGGCGGCGGAAAACGCCGCATCGCTGTTTGCCCTGCCGGCAGCATTGCTGGCGCGCGATGATCATGATGCCGATTTCTCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGGAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACAtcgtcg  >  1:98454/1‑251 (MQ=255)
                                                                                                                                                                           ggcggcGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTca                   >  1:14995/1‑218 (MQ=255)
                                                                                                                                                                                cggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCgaga                                                                                                          >  1:202456/1‑126 (MQ=255)
                                                                                                                                                                                 ggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCAGCACAGNNGCTGGan                                                                                                                                      >  2:13799/1‑96 (MQ=255)
                                                                                                                                                                                 ggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATCAGCACAGCTGCTGGAAgg                                                                                                                                    <  1:13799/99‑1 (MQ=255)
                                                                                                                                                                                                gCTGTTGGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTGATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAAGTCCTCGGACTTTACGGGGAGAAAGTAAGCGTAACAGGATCACGACAATTCTGACGATTCATGTTCCTTCAACGCCgg                                                  <  1:146154/166‑1 (MQ=255)
                                                                                                                                                                                                         ccTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATTTATCGGTTAAATAAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGgcg                                              >  1:235140/1‑161 (MQ=255)
                                                                                                                                                                                                                        gCTGGCGCGCGATGATCATGATGCAGATGTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCGCGTACTTTACGGCGGGAAAGTAAACGGAACTGGTTGACGACATGTCTGACGATTCATGTTCCt                                                            <  1:84907/132‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                 
GCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGTAACATGATGACGACAATTCTGACGATTCATGTTCCTTCAACGCCGGGGCGCGCATGGAATATGCTGGTGGCACTTCAGGCAGGAAACATCGTCG  >  NC_000913/657940‑658345

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