Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 296,504 T→G N198T (AAC→ACC)  intF ← CP4‑6 prophage; putative phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913296,5040TG85.7% 10.9 / ‑2.4 7N198T (AAC→ACC) intFCP4‑6 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (6/0);  total (6/1)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

TTGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCAAACTAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTCC  >  NC_000913/296443‑296568
                                                             |                                                                
ttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGa                                                               <  1:1763952/65‑1 (MQ=255)
    cTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCaa                                                           >  1:3145256/1‑65 (MQ=255)
            aCTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCGAAAAAGccc                                                   >  1:258161/1‑64 (MQ=255)
            aCTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGAGGCAAATAGccc                                                   >  1:3964243/1‑64 (MQ=255)
            aCTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGGGGGGAAACTAGCCa                                                   >  1:351086/1‑65 (MQ=255)
              tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTGGGGCAAACTAGCCaaa                                                 >  1:689256/1‑64 (MQ=255)
              tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGGAAATAGCCAAt                                                 >  1:2658794/1‑65 (MQ=255)
              tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAATACCCAAt                                                 >  1:976972/1‑65 (MQ=255)
              tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGAGCGCCAACTCGCCAAt                                                 >  1:4196246/1‑65 (MQ=255)
                 gCCACGCTGCTCTGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTTAACAAACTAGCCAAtcgg                                              >  1:112018/1‑62 (MQ=255)
                          ctATGTAAGCGGGGTTTAGCTCCGATAACGGGAGGTTGAGAAAATTAGCCAATGGTCCAGCCGaa                                     <  1:2133525/65‑1 (MQ=255)
                            aTGTAATCTGGGGTTCGCTCCGATACGGGCGGGTTGAGCAAACTAGCCAATGGTCCAGCCGAAGt                                   <  1:635566/65‑1 (MQ=255)
                              gTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGGGCAAACCAGCCACTGGGCCACCCGCAGTCg                                 >  1:3160801/1‑65 (MQ=255)
                              gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGAGGAAAATAGCCAAAGGGCCAGCCGCAGGCg                                 >  1:4550305/1‑65 (MQ=255)
                              gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGAGCAAACTAGCCCATGGGCCCGCCGCAGACg                                 >  1:3314564/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAACTAGCCCATGGTCCAGCCGGAGGCggg                               >  1:2806973/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAACTAGCCAAAGGGCCAGCCGGAGGCggg                               >  1:583989/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAATACCCCACGGTCAGGTGGAAGTCggg                               >  1:3433212/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAATAAGCAAAGGGGCAGGCCAAGGCggg                               >  1:682974/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAACTAGCCAATGGGCCCGCCGGAGTCggg                               >  1:884939/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAACTAGCCAATGGGCCAGCCGAAGTCggg                               >  1:4358089/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAACAGCCCACGGGGCCGCCCCAGTCggg                               >  1:1191314/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGAGCAAACTAGCCAATGGGCCAGCCGCAGGCggg                               >  1:696947/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGAGCAAACTAGCCAATGGGCCAGCCGAAGGCggg                               >  1:631568/1‑65 (MQ=255)
                                aaTCTGGGGTTAGCTCCGATAACGGCAGGGGGGGCAAACCAGCCCAAGGTGCCGCCGCAGGCggg                               >  1:2187182/1‑65 (MQ=255)
                                     ggggTTAGCTCCGATAACGGCAGGGTGAGGAAAATAGCGAAAGGGCCAGCCGAAAGCGGGGCTTg                          >  1:960479/1‑65 (MQ=255)
                                        gTTAGCTCCGATAACGGCAGGGTGAGCAAACTAGCCCATGGTCCAGCCGAAGTCGGGGCGTGGcc                       >  1:1640081/1‑65 (MQ=255)
                                         ttAGCTCCGATAACGGCAGGGTGGGGAAACTAGCCAAAGGGCCAGCCCGAGGGGGGCCTTTGccg                      >  1:23983/1‑64 (MQ=255)
                                         ttAGCTCCGATAACGGCAGGGTGGGGAAACCATCCCACGGTTCTGCCCAACTCGGGGCGTGGcgg                      >  1:619719/1‑63 (MQ=255)
                                            gCTCCGATAACGGCGGTTTGAGCAAATTAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCtt                   <  1:1339100/65‑1 (MQ=255)
                                            gCTCCGATAACGGCAGGGTGAGGAAAATAACCCATGGGCCAGCCGGAGGCGGGGCCTTGcttctt                   >  1:475332/1‑65 (MQ=255)
                                              tCCGATAACGGCAGGGTGGGGAAAATAGCCAATGGGCCAGCCCCAGGCGGGGCCTGGGCTCtttt                 >  1:3453580/1‑65 (MQ=255)
                                                           gattGAGCAATCTAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCt    <  1:3576739/63‑1 (MQ=255)
                                                             ttGAGCAAACTAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTcc  >  1:4410896/1‑65 (MQ=255)
                                                             |                                                                
TTGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCAAACTAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTCC  >  NC_000913/296443‑296568

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: