Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,529 T→C H152R (CAC→CGC)  yagA ← CP4‑6 prophage; integrase core domain‑containing protein YagA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,5290TC95.2% 64.8 / NA 22H152R (CAC→CGC) yagACP4‑6 prophage; integrase core domain‑containing protein YagA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (8/12);  total (9/13)
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

TGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGC  >  NC_000913/281407‑281663
                                                                                                                          |                                                                                                                                      
tGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGGGGATGGCAGCGTGCACGGCCAAAAGGAAAGAGGCCCTTAAAATCCATc                                                                                                                       >  2:195360/1‑140 (MQ=255)
          ctgctgcACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGc                                                                                                             >  1:475732/1‑140 (MQ=255)
                      ctcGCGCCGTTCATCTGTACAGTGCCCCATGCACAGGCAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGcgt                                                                                                 <  1:195360/140‑1 (MQ=255)
                      ctcGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGcgt                                                                                                 <  2:69375/140‑1 (MQ=255)
                       tcccgcCGTTCATCGGTATAGTGCGCCAGGCAGAGCGAGAAACGGGAGTGGTCGTCCAGCAGAGTGAGCGCATGGCAGCGTCCCCCGCCAAAAGGAAAGCGGCCCTTAAAATCCATATGCCAGAGGCGGTTCGGCGcgtc                                                                                                <  2:583276/137‑1 (MQ=255)
                               ttCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGa                                                                                        >  1:114919/1‑140 (MQ=255)
                                     ggTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGc                                                                                  <  2:475732/140‑1 (MQ=255)
                                          aGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGGGTTCGAACCGGCCCGTg                                                                             >  2:435843/1‑140 (MQ=255)
                                            tGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATGCATCTGCCGGAGGCGGGTCGGCGCGGCGTGATCGAACCGGCCCGTggg                                                                           >  2:185164/1‑139 (MQ=255)
                                            tGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGc                                                                           >  2:376022/1‑140 (MQ=255)
                                              cgcCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCTTCTGCCAGGGGCGGGTAGCCGCGCCACGTTTGAACAGGCCCGGTGCgg                                                                         >  2:169514/1‑140 (MQ=255)
                                              cgcCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCgg                                                                         >  1:43483/1‑140 (MQ=255)
                                              cgcCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCCTGGCgg                                                                         >  2:20629/1‑140 (MQ=255)
                                              cgcCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGGCGTGTTCGAACCCGCCCGGGGCgg                                                                         >  2:375835/1‑140 (MQ=255)
                                                ccAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGa                                                                       <  2:114919/140‑1 (MQ=255)
                                                             aaaaCGGGAGTGGTCGTCCAGCAGGGTGAGCAGATCGCATCGTCCACCACCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAg                                                          <  2:242652/140‑1 (MQ=255)
                                                                       tGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGgcag                                                <  2:566187/140‑1 (MQ=255)
                                                                                         ggCGGATGACAGCGTCGTCCGCAAAAAGGATGGTGGCCCATAAAATCCTTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCAt                              <  2:325933/139‑1 (MQ=255)
                                                                                          gCGGATGGCGGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATc                             <  2:483637/140‑1 (MQ=255)
                                                                                          gCGGATGGCAGCGTCTACCGCCAAACGGACAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATc                             <  1:185164/140‑1 (MQ=255)
                                                                                          gCGGATGGCAGCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATc                             <  2:68811/140‑1 (MQ=255)
                                                                                                    gCGTCCACCGCCAAAAGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGa                   <  2:140257/140‑1 (MQ=255)
                                                                                                                aaaaGGAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAgg       <  1:376022/140‑1 (MQ=255)
                                                                                                                     gAAAGCGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGc  <  1:375835/140‑1 (MQ=255)
                                                                                                                          |                                                                                                                                      
TGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGC  >  NC_000913/281407‑281663

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