Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,848,300 2 bp→TT coding (376‑377/552 nt) ydjA ← putative oxidoreductase YdjA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,848,3000CT92.9% 31.4 / ‑5.9 14G126K (GGG→AAG) ydjAputative oxidoreductase YdjA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/5);  minor base G (1/0);  total (9/5)
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
*NC_0009131,848,3010CT100.0% 39.7 / NA 15G126K (GGG→AAG) ydjAputative oxidoreductase YdjA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/5);  total (10/5)

GCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAAT  >  NC_000913/1848115‑1848540
                                                                                                                                                                                         ||                                                                                                                                                                                                                                               
gcgACAATTATCAGAAATAAGTTACANACGGCGNNGGNTCCGGGACGTTAATCGACGNAGNNNCNNNCNGCTGCNGCNNACCGAGGTANAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACNGGACTTTCANNTAATGCGCCACTGCGCCAGATNCCGNNANNCNTCTnggn                                                                                                                                                                                                                                           >  2:83910/1‑192 (MQ=255)
                                       ccGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCcaccac                                                                                                                                           <  1:83910/250‑1 (MQ=255)
                                         gggACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTc                                                                                                                                                          <  1:24550/233‑1 (MQ=255)
                                                                                          gaAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCCCTGCGCCAGCTGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAAc                                                                                                                  >  1:316018/1‑224 (MQ=255)
                                                                                                   naaTTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCcgcg                                                                                                                                                                               >  1:66695/2‑154 (MQ=255)
                                                                                                          aTCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTATTACCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGCTAGCTTTGTCATCACt                                                                            >  1:203526/1‑246 (MQ=255)
                                                                                                             cTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCCAACTTCTGGGCCACTGCTGCCCTTTGCATCGCCATGTCCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGCTAGCTTTGTCATCACTACCGgc                                                                      >  1:191865/1‑249 (MQ=255)
                                                                                                               gCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGAATCCGGCAGCCCTTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTgcgc                                                                                                                       >  1:315812/1‑198 (MQ=255)
                                                                                                               gCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACGGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGgcagc                                                                   <  1:199783/250‑1 (MQ=255)
                                                                                                                                 cTTCACGCACTACCGGACTTTCAGTTAAGGCGCCACTGCGCCAGAGGCCGCAAAACTTCTGGGCATCGGCTGCCATTTGCATCGCCATGACCGCGCATCCGTCAGACATTTCCAGGTCCCAGCGCGGGACGTTATGATCCTCTTCGCATGGCGCCTCCACCGTGATAATGAGCGGTGCGCGGAACGgc                                                                                                               <  1:172284/188‑1 (MQ=255)
                                                                                                                                   tCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGc                                                                                                                                                                                                                           <  1:20633/78‑1 (MQ=255)
                                                                                                                                   tCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGc                                                                                                                                                                                                                           >  2:20633/1‑78 (MQ=255)
                                                                                                                                                             tgcgcacatGCGCCAGAGGCCGCCAAACGTCTGGGCAATTGCCGCCATTTACATCGCCACAACCGCGCATCCGGAAGACATTTCCTATTCCCAGAACGGGAATTTATGATTCTCTCCGCCTTTCACCACCACCGCGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGATTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCACTGAAACGCTCGCGcccagctc                    >  1:47972/9‑246 (MQ=255)
                                                                                                                                                              gcgcCACTGCGCCAGATGCCGCCAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCgg                                                                                                                            >  1:327542/1‑146 (MQ=255)
                                                                                                                                                                               gccgccAAACTTCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAAt  >  1:161692/1‑251 (MQ=255)
                                                                                                                                                                                         ||                                                                                                                                                                                                                                               
GCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGAGAAAACCGACAATTTTATCCTGCTCACGGCAACCGAATGCTTCACGCACTACCGGACTTTCAGTTAATGCGCCACTGCGCCAGATGCCGCCAAACCCCTGGGCAACTGCTGCCATTTGCATCGCCATGACCGCGCATCCGGCAGACATTTCCTGTTCCCAGCGCGGGACTTTATGATTCTCTTCGCATTTCGCCACCACCGTGATGATGAGCGGTGCGCGGAACGGCGCATTACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAAT  >  NC_000913/1848115‑1848540

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