Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,316,498 G→A H241Y (CAC→TAC)  phnM ← RPnTP hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,316,4980GA93.3% 37.3 / NA 15H241Y (CAC→TAC) phnMRPnTP hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (7/7);  minor base C (1/0);  total (8/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

GGCGACGCGAAATGCCGCATCGAGCAGGCTGGCGGGGTAGTAGTCGGAAGAGAGGATATCCAGCAGGCCAAGCTGCGCCAGTTCACTGGCCGCCACGTTGCCGGAGTGCGAGCCGCCGCGCACAATATTCGGCGCGCCCATCAGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTGAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGG  >  NC_000913/4316274‑4316701
                                                                                                                                                                                                                                |                                                                                                                                                                                                           
ngCGACGCGAAATGCCGCATCGAGCAGGCTGGCGGGGTAGTAGTCGGAAGAGAGGATATCCAGCAGGCCAAGCTGCGCCAGTTCACTGGCCGCCACGTTGCCGGAGTGCGAGCCGCCGCGCACAATATTCGGCGCGCCCATCAGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCgtgggtgg                                                                                                                                                                                   >  1:37373/2‑251 (MQ=255)
       ngAAATGCCGCATCGAGCAGGCTGGCGGGGTAGTAGTCGGAAGAGAGGATATCCAGCAGGCCAAGCTGCGCCAGTTCACTGGCCGCCACGTTGCCGGAGTGCGAGCCGCCGCGCACAATATTCGGCGCGCCCATCAGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATc                                                                                                                                                                            >  1:214538/2‑251 (MQ=255)
                        cAGGCTGGCGGGGTAGTAGTCGGAAGAGAGGATATCCAGCAGGCCAAGCTGCGCCAGTTCACTGGCCGCCACGTTTCCGGAGTGCGAGCCGCCGCGCACAATATTCGGCGCGCCCATCAGCACGTTCATGCCATGCTTTCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCTATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCa                                                                                                                                                            <  1:19516/250‑1 (MQ=255)
                         aGGCTGGCGGGGTAGTAGTCGGAAGAGAGGATATCCAGCAGGCCAAGCTGCGCCAGTTCACTGGCCGCCACGTTGCCGGAGTGCGAGCCGCCGCGCACAATATTCGGCGCGCCCATCAGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAAt                                                                                                                                                          <  1:122604/251‑1 (MQ=255)
                                                                                                                 cgccgcGCACAATATTCGGCGCGCCCATCAGCACGTTCATGGCATGCTTTCGCTAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTTGCAAGCTGGTAAGATTCAGCAACGTGTGCGTGGGTGGCGTCATCGTGGCTGGCAAGGGCAATTTTTCGCGCGTGGCGCATGGCGGCTATGGTTTCGCGATTCGGCTGCTACGAGCTTGCGGCGAGCGCCAGTTGCTCTTCTTCGTActgc                                                                  <  1:138988/251‑1 (MQ=255)
                                                                                                                  gccgcGCACAATATTCGGCGCGCCCATCAGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGATGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCg                                                                                                           <  1:13354/209‑1 (MQ=255)
                                                                                                                          naaTNTTCGGCGCGCCCATCAGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGt                                                                                                                                                                          <  2:214538/137‑1 (MQ=255)
                                                                                                                                    cgnnnNNATCNGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGACGATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAAt                                                                                                                                                          <  2:211225/144‑1 (MQ=255)
                                                                                                                                      ngcCCATCAGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAgg                                             >  1:90877/2‑251 (MQ=255)
                                                                                                                                        cccATCAGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGCAATTCGGCGATCACGCTGCCACGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTgcggc                                                                                              >  1:133874/1‑200 (MQ=255)
                                                                                                                                               gCACGTTCCTGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCATAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTg                                    >  1:353667/1‑251 (MQ=255)
                                                                                                                                                  ngTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGccc                                 >  1:285492/2‑251 (MQ=255)
                                                                                                                                                                         gCTTCNGCCGCNTNNAACGCGGAGGGAAATTAGGNGATCACGCTGCCAAGNTGGTNAGAttncgn                                                                                                                                                                                                    >  2:137752/1‑61 (MQ=255)
                                                                                                                                                                            tCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTAAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGGGACCAGCGTGCGGCGAGCGCCAGTTGTTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGATGGAGTATTTGCCCTGATAATATTCGCGATACTTCTcg         <  1:137752/249‑1 (MQ=255)
                                                                                                                                                                                           gtggtgGGAAATTCGGCGATCACGCTGCCAAGCTGGTCAGCTTCAGCAACGTGGGCGTGGCCGCTGTCATCGTGGCTGGCACGCGCAATTTTTTGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGc                                                                                                               >  1:195717/1‑132 (MQ=255)
                                                                                                                                                                                                  gAAATTCGGCGATCACGCTGCCAAGCTGGTAAGATTCAGCGACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCACTTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGAGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCCTCTGCGCATCAGTGAGGGCGTATTTGCCCTGATAATATTCGCGCTACTTCTCGCGGTTgg  >  1:227786/1‑234 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                           
GGCGACGCGAAATGCCGCATCGAGCAGGCTGGCGGGGTAGTAGTCGGAAGAGAGGATATCCAGCAGGCCAAGCTGCGCCAGTTCACTGGCCGCCACGTTGCCGGAGTGCGAGCCGCCGCGCACAATATTCGGCGCGCCCATCAGCACGTTCATGCCATGCTTGCGCGAGGCTTCCGCCGCTTCGAACGTGGTGGGAAATTCGGCGATCACGCTGCCAAGCTGGTGAGATTCAGCAACGTGGGCGTGGGTGGCGTCATCGTGGCTGGCAAGCGCAATTTTTCGCGCGCGGCACAGGGCGGCGATGGATTCGCGATTCGGCTGCGACCAGCGTGCGGCGAGCGCCAGTTGCTCTTCTTCGTACTGCTGCATCTGCGCATCAGTGAGGGAGTATTTGCCCTGATAATATTCGCGATACTTCTCGCGGTTGG  >  NC_000913/4316274‑4316701

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