Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,376,438 G→A S41N (AGC→AAC)  ecnA → entericidin A lipoprotein, antidote to entericidin B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,376,4380GA88.9% 16.9 / ‑5.9 9S41N (AGC→AAC) ecnAentericidin A lipoprotein, antidote to entericidin B
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (2/6);  total (2/7)
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

GTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCATCTCTCGCGCTGCCAGCTAATTTTTCTTCTCTTCCGAAAAATCATCAGATTCCCATCATTTTTGGCGATGTTGTCTATTATTAATTTGCTATAGGCAAACATAAATAACATTACCTAAAAGGAAGACGTTATGGTGAAGAAGACAATTGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGC  >  NC_000913/4376270‑4376677
                                                                                                                                                                        |                                                                                                                                                                                                                                               
gTTGTGGTGCGGCCTGCAGGCTGGACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTTTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGATGTCATCAAACAGCTCGGCAACCCCATCGCTCGCGCTGCCAACTAATTTCTCGTCTCTTCCGAAAAATCTTCAGATTGCCAGCATTTTTGGCGATGTTGTCtatt                                                                                                                                                                                  <  1:214405/232‑1 (MQ=255)
gTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCATCTCTCGCGCTGCCAACTAATTTTTCTTCTCTTCCGAAAAATCATCAGATTCCCATCATTTTTGGCGAtgtn                                                                                                                                                                                         <  1:155434/225‑2 (MQ=255)
                                    tCAGGGCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCATCTCTCGCGCTGCCAACTAATTTTTCTTCTCTTCCGAAAAATCATCAGATTCCCATCAGTTTTGGCGATTTTGTCTATTAGTAAGTTGCTATATGCAAACATAAATAACATTCCCTAAAAGGAAGACGTTATg                                                                                                                            <  1:182041/250‑1 (MQ=255)
                                                                                         ntgctCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCATCTCTCGCGCTGCCAACTAATTTTTCTTCTCTTCCGAAAAATCATCAGATTCCCATCATTTTTGGCGATGTTGTCTATTATTAATTTGCTATAGGCAAACATAAATAACATTACCTAAAAGGAAGACGTTATGGTGAAGAAGACAATTGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTAtt                                                                       >  1:156763/2‑250 (MQ=255)
                                                                                                  ggcggTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCGCGGCAACTCTATCTCTCGCGCTGCCAGCTAATTTTTCTTCTCTTCCTAAAAATCATCAGATTCCCATCATTTTTGGCGAtgttgt                                                                                                                                                                                       <  1:188518/125‑1 (MQ=255)
                                                                                                  ggCTGNNACNNNGCNNGCGGTTTCGGCGAAGACATCAAANATCTCGGCAACTCCATCTCTCGCGCTGCCNANTNCTTTNNCTTNTCTTCCGAAANATCATCAGATTCCCATCATTTTTGGCGAtgttgt                                                                                                                                                                                       >  2:188518/1‑129 (MQ=255)
                                                                                                                          ggCGAAGACATCAAACATCTCGGTAACTCCATCTCTCGCGCTGCCAACTAATTTTTCTGCTCTTCCGAAAAATCATCAGATTCCCATCATTTTTGGCGATGTTGTCTATTATTAATTTGCTATAGGCAAACATAAATAACATGACCTAAAAGGAAGACGTTATGGTGAAGAAGACAAGTGCAGCGATCTTGTCTGTTCTGGCGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTgn                                               <  1:53302/241‑2 (MQ=255)
                                                                                                                                                              cgcgCTGCCAACTAATTTTTCTTCTCTTCCGAAAAATCATCAGATTCCCATCATTTTTGGCGATGTTGTCTATTATTAATTTGCTATAGGCAAACATAAATAACATTACCTAAAAGGAAGACGTTATGGTGAAGAAGACAAGTGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGgcgc  <  1:83916/250‑1 (MQ=255)
                                                                                                                                                                    gCCAACTAATTTTTCTTCTCTTCCGAAAAATCATCAGATTCCCATCATTTTTGGCGATGTTGTCTATTATTAATTTGCTATAGTCAAACATAAATAACATTACCTAATAGTAATACGTTATTGGGAAGAAGACAATTGCAGCGTTCTTTTCTGTTCTGGTGCTTTCAACAGTATTAActgc                                                                 <  1:110930/181‑1 (MQ=255)
                                                                                                                                                                        |                                                                                                                                                                                                                                               
GTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCGGCAACTCCATCTCTCGCGCTGCCAGCTAATTTTTCTTCTCTTCCGAAAAATCATCAGATTCCCATCATTTTTGGCGATGTTGTCTATTATTAATTTGCTATAGGCAAACATAAATAACATTACCTAAAAGGAAGACGTTATGGTGAAGAAGACAATTGCAGCGATCTTTTCTGTTCTGGTGCTTTCAACAGTATTAACTGCCTGCAACACCACGCGTGGCGTTGGTGAAGACATTTCTGATGGCGGTAACGCGATTTCTGGCGC  >  NC_000913/4376270‑4376677

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