Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,413,535 G→A E170K (GAA→AAA)  yjfC → putative acid‑‑amine ligase YjfC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,413,5350GA94.1% 40.2 / ‑6.2 17E170K (GAA→AAA) yjfCputative acid‑‑amine ligase YjfC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (8/8);  minor base C (1/0);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GAAAGCTGGCGCGCCCGCGATCCTTCGCTGTATGGACGGATGGATTTTGCCTGGTGTGGCAATGCGCCGGTGAAGCTGCTGGAGTACAACGCCGATACGCCAACTTCATTGTACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGGAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACTTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTCGCGGTTTATCTACATTGAAGATCTCGGTTTGGGCGTCGGCGGCGTACTGACCGATCTTGATGATAATGTCATCCAGCGTGCATTTAAGCTGTATCCGCTGGAGTGGATGATGCGTGACGATAACGGTCCGCTGCTGCGCAAGCGTCGC  >  NC_000913/4413292‑4413777
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
gccacaatccgcgcCCGCGATCCTTCGCTGTATGGACGGATGGATTTTGCCTGGTGTGGCAATGCGCCGGTGAAGCTGCTGGAGTACAACGCCGATACGCCAACTTCATTGTACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGtt                                                                                                                                                                                                                                             >  1:317466/10‑251 (MQ=255)
        gcgcgcCCGCGATCCTTCGCTGTATGGACGGATGGATTTTGCCTGGTGTGGCAATGCGCCGGTGAAGCTGCTGGAGTACAACGCCGATACGCCAACTTCATTGTACGAGTCGGCTTATTTACAGTGGCTGTGGCTGGAGGATGCCCGGTGCAGCGGCATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTAttttt                                                                                                                                                                                                                                     >  1:441290/1‑251 (MQ=255)
                     ccTTCGCTGTATGGACGGATGGATTTTGCCTGGTGTGGCAATGCGCCGGTGAAGCTGCTGGAGTACAACGCCGATACGCCAACTTCATTGTACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGTTCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTATTTGTGCTGCTGTCAgg                                                                                                                                                                                                                         <  1:119249/250‑1 (MQ=255)
                                                         ggCAATGCGCCGGTGAAGCTGCTGGAGTACAACGCCGATACGCCAACTTCATTGTACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACTTgn                                                                                                                                                                                     <  1:47824/250‑2 (MQ=255)
                                                                    ggTGAAGCTGCTGGAGTACAACGCCGATACGCCAACTTCATTGTACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGCTCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTCCTTGCAGGACTGCGccc                                                                                                                                                                         >  1:231711/1‑251 (MQ=255)
                                                                          gctgctGGCGTACCCCGCCGATACGCNACCTTCNNTGNACGCGTCGGCTTATTCCCANCGNNNGNNNCNTGACGNTGNNCGGCGCCGCNGCCTTATTCCGCGTGATGCTGTTCTGTAACCA‑GCTACTCAGGAACGCCTGGTTNCACGCGTTTNNGGA‑CCTTACAGTTCGCAACCGNTTTNNTNNtnctgncgn                                                                                                                                                                                                                             >  2:95854/1‑189 (MQ=255)
                                                                                           ccGATACGCCAACTTCATTGTACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACTTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTcgcgn                                                                                                                                                  <  1:266533/251‑2 (MQ=255)
                                                                                                          cgttGTACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAAACGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACTTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTCGCGGTTTAt                                                                                                                                             >  1:411396/3‑241 (MQ=255)
                                                                                                                                 cAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTAGCGTGCTGTACTTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTCGCGGTTTATCTACATTGAAGATCTCGGTTTGGGCGTcggcgg                                                                                                            <  1:444725/251‑1 (MQ=255)
                                                                                                                                           ggctgGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGATCCGTA‑CAATGATATCCAAGAAAGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTATTATTTTTGTTGCTGTccggcc                                                                                                                                                                                                                       >  1:197068/1‑129 (MQ=255)
                                                                                                                                                           gcgcAGCGGCATTATTCCGCGTGATGCCGATAAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTTTACTTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTCGCGGTTTATCTACATTGAAGATCTCGGTTTGGTCGTCAGCGGCGTACTGACCGATCATGATGATAATg                                                                                  <  1:372237/251‑1 (MQ=255)
                                                                                                                                                             gcAGCGGCATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCCCGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGa                                                                                                                                                                                                       >  1:257177/1‑132 (MQ=255)
                                                                                                                                                                   ncATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACTTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTCGCGGTTTATCTACATTGAAGATCTCGGTTTGGGCGTCGGCGGCGTACTGACCGa                                                                                                >  1:90075/2‑229 (MQ=255)
                                                                                                                                                                                                          cAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGAAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACTTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTCGCGGTTTATCTACATTGAAGATCTCGGTTTGGGCGTCGGCGGCGTACTGACCGATCTTGATGATAATGTCATCCAGCGTGCATTTAAGCTGTGTCCGCTGGAGTGGATGATGc                                     <  1:146547/250‑1 (MQ=255)
                                                                                                                                                                                                                  ccTGATTTCACGA‑TTTAGCGAGCTTTACAGTCGGAACCCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACTTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTCGCGGTTTATCTACATTGAAGATCTCGGTTTGGGCGTCGGCGGCGTACTGAcc                                                                                                  >  1:341358/1‑181 (MQ=255)
                                                                                                                                                                                                                                   aGCGAGCTTTACAGTCGGAAACCGTTTTATGTGTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACTTGCAGGACTgcgc                                                                                                                                                                           <  1:167784/92‑1 (MQ=255)
                                                                                                                                                                                                                                             aCAGTCGGAAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACGTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTCGCGGTTTATCTACATTGAAGATCTCGGTTTGGGCGTCGGCGGCGTACTGACCGATCTTGATGATAATGTCATCCAGCGTGCATTTAAGCTGTAGCCGCTGGAGTGGATGATGCGTGACGATAACGGTCCGCTGCGGCGCAAGCGTcgc  <  1:95854/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
GAAAGCTGGCGCGCCCGCGATCCTTCGCTGTATGGACGGATGGATTTTGCCTGGTGTGGCAATGCGCCGGTGAAGCTGCTGGAGTACAACGCCGATACGCCAACTTCATTGTACGAGTCGGCTTATTTCCAGTGGCTGTGGCTGGAGGATGCCCGGCGCAGCGGCATTATTCCGCGTGATGCCGATCAGTA‑CAATGCTATTCAGGAACGCCTGATTTCACGC‑TTTAGCGAGCTTTACAGTCGGGAACCGTTTTATTTTTGCTGCTGTCAGGACACCGATGAAGACAGGAGTACCGTGCTGTACTTGCAGGACTGCGCCCAGCAGGCAGGGCAGGAGTCGCGGTTTATCTACATTGAAGATCTCGGTTTGGGCGTCGGCGGCGTACTGACCGATCTTGATGATAATGTCATCCAGCGTGCATTTAAGCTGTATCCGCTGGAGTGGATGATGCGTGACGATAACGGTCCGCTGCTGCGCAAGCGTCGC  >  NC_000913/4413292‑4413777

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