Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,174 C→T L270L (CTG→CTA yagA ← integrase core domain‑containing protein YagA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,1740CT90.9% 22.7 / ‑4.3 11L270L (CTG→CTAyagAintegrase core domain‑containing protein YagA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (4/6);  minor base G (0/1);  total (4/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 = 7.48e-01

TCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGAGTGCCCCACCCGGATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGC  >  NC_000913/280933‑281406
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        
tCAGCCCGCCCCGTTCCCCCCTGAACCCCTGCCCGGCGCTCAGACTTCCCCCTTCCCCGCTCAGCTTTCCGCTGATATCCCCTTCCCTGACCATCCCCCCTTCACCGTATTCCGGGGGCGTTGTGTGGCGGCTGTACGGCCGCGCTGACGGCTGATACCGCGAGCGCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGTAGTTCAc                                                                                                                                                                                                                                   <  1:81202/249‑1 (MQ=255)
        ccccGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTCCCCCTTTCCCGCTCAGCTTTCCGCTGATATCCCCTTTCCTGCCCATCCCCCCTTCATCGTATTCCGGGGGCTTTGTGTTCCCGCTGTCCTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGGAGTTCACCGCTGTCTGc                                                                                                                                                                                                                         <  1:45544/250‑1 (MQ=255)
                   cctgacccccTTGCCTGCGCTCAGACTTCCCCCTTCCCCGCTCAGCTTTCCGCTGATATCCCCTTTCCTTCCCATCCCCCCTTCATCGTATTCGGGGGGCGTTGTGTTGCCGCTGTCCTGCCGCGCTGACGGCTGATCCCGCGGGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGTAGTTCACCGCTGTCTGCGAACCATTTTc                                                                                                                                                                                                              <  1:9452/243‑1 (MQ=255)
                                                                                                     ccATCGTTTTCGGGGGGGTTTGTGTTGCCGCTGTACGGCCGGGCGGCCGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGTAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTTCGTCTGCGGATGATAAGGCCGGGAGTGc                                                                                                                              <  1:53320/248‑1 (MQ=255)
                                                                                                          gTATTCCGGGGGCGTTGTGTTGCCGCTGTACTTCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGTAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGAGTGCCCCAccc                                                                                                                       >  1:56512/1‑251 (MQ=255)
                                                                                                                                ccGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGCCCGTCCGCCAGTGGTCGAAGGCGCGCTGTAGTCCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGCGTGCCCCACCCGGATACCCAGGCGCATCAGcc                                                                                                  >  1:106348/1‑250 (MQ=255)
                                                                                                                                                                                  atCCAGCGCCTCGTGCGGGCGTTCAAGGTGATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGTAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATTATAAGGCCGGGAGTGCCCCACCCGGATACCCAGGCTCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGt                                               <  1:111989/251‑1 (MQ=255)
                                                                                                                                                                                                   ggCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGTAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGNGGANNNNNAgnnnnnnnnnn                                                                                                                               >  2:101759/1‑144 (MQ=255)
                                                                                                                                                                                                      gTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGTAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGAGTg                                                                                                                               <  1:81449/151‑1 (MQ=255)
                                                                                                                                                                                                                              aGTGGTCGACGGCGCGATGTAGTTCACCGCTGCCTGCGAACAATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAACACGCTCCAGCTTCCCCTGCGTCTGCGGATGCTAAGGCCGGGAGTGCCCCACCCGGATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGATCAAACACg   >  1:139815/1‑251 (MQ=255)
                                                                                                                                                                                                                               gTGGTCGAAGGCGCGCTGTAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGAGTGCCCCACCCGGATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGTAGGCCGTAACGCTCAAACACgn  <  1:101759/251‑2 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        
TCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGTTCACCGCTGTCTGCGAACCATTTTCCCTGCAGCACTTCCGCCTTCAGGCTGCGGTGAAAACGCTCCAGCTTCCCCTGCGTCTGCGGATGATAAGGCCGGGAGTGCCCCACCCGGATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGC  >  NC_000913/280933‑281406

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