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,5040TG80.0% 15.7 / 1.3 10N198T (AAC→ACC) intFCP4‑6 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (8/0);  total (9/1)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-02

ATTTTGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCAAACTAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTCCACG  >  NC_000913/296440‑296571
                                                                |                                                                   
aTTTTGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGa                                                                  <  1:1153997/68‑1 (MQ=255)
  tttGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGc                                                                >  1:3989730/1‑68 (MQ=255)
                 tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAACGAGACAATggg                                                 >  1:3414574/1‑67 (MQ=255)
                 tCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAACCAGCCAATggg                                                 >  1:2941929/1‑67 (MQ=255)
                    gCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGGAACTAGACAAAGGTCCa                                              >  1:1620588/1‑68 (MQ=255)
                    gCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAACTAGCCAATGGGccc                                              >  1:2702691/1‑67 (MQ=255)
                      cACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAATAGCCAAAGGGGCAGc                                            >  1:2189810/1‑68 (MQ=255)
                         gctgctATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGCAAAAAGCCAAAGGGCCCGccca                                         >  1:3127572/1‑66 (MQ=255)
                         gctgctATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAATAACCAATGGGCCAGCCGa                                         >  1:2064454/1‑68 (MQ=255)
                              tATGTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAACTAGCCCAAGGGGCAGCCGCAgagg                                    >  1:3427622/1‑65 (MQ=255)
                                 gaaaTCGGGGTTTACCCCCGATATCGGCAGGTTGACCAACCTAGCCAATGGTCCAGCCGAAGTCGGGc                                 <  1:552793/66‑1 (MQ=255)
                                 gTATTCTGGGTTTGGCCCCAAAAACGGCGGGTGGAGAAAATTAGCCAATGGTCCAGCCGAAGTCGGGc                                 <  1:2781311/68‑1 (MQ=255)
                                 gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAACTAGCCAATGGGTCCGCCGGAGGCgggg                                 >  1:1061931/1‑67 (MQ=255)
                                 gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAATAACCAAAGGTCCAGCCGAAGGCGGGc                                 >  1:1558729/1‑68 (MQ=255)
                                 gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAACTAACCAATGGTCCCGGCGGAGGCgggg                                 >  1:809730/1‑67 (MQ=255)
                                 gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGAGGAAACTAGCCAATGGGCCAGCCGAAGGCGGGc                                 >  1:306735/1‑68 (MQ=255)
                                 gTAATCTGGGGTTAGCTCCGATAACGGCAGGGTGAGCAAAATAGCCAATGGGCCAGCCGAAGTCGGGc                                 >  1:2860778/1‑68 (MQ=255)
                                   aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAAAAGCCAATGGGCCAGCCGGAGGCGGGGCt                               >  1:796012/1‑68 (MQ=255)
                                   aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCCAAACAGCCCATGGTGCAGCCGCAGTCGGGGCt                               >  1:546895/1‑68 (MQ=255)
                                   aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAACTAGCCCACGGTCCCCCCCCACTCGGGCCt                               >  1:1835633/1‑68 (MQ=255)
                                   aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAACAAACCAATGGGCCAGGCCCAGGCGGGcgt                               >  1:3140381/1‑66 (MQ=255)
                                   aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAATAACCAATGCGGCAGCCGAAGGCGGGGCt                               >  1:3595917/1‑68 (MQ=255)
                                   aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGGGCAAAACAGCCCAAGGGCCAGCCGCAGGCggggcc                               >  1:1993622/1‑65 (MQ=255)
                                   aaTCTGGGGTTAGCTCCGATAACGGCAGGGTGAGCAAACCAGCCAATGGGCCACCCGAAGGCGGGCCt                               >  1:2342944/1‑68 (MQ=255)
                                   aaTCTGGGGTTAGCTCCGATAACGGCAGGGGGAGGAAAATAGCCCATCGGCCAGCCCCAGGCGGGCCt                               >  1:550134/1‑68 (MQ=255)
                                      cTGGGGTTAGCTCCGATAACGGCAGGGTGGGGAAAATAGCCAATGGGCCAGCCGAAGGCGGGGCtttg                            >  1:2335003/1‑66 (MQ=255)
                                          ggTTAGCTCCGATAACGGCAGGGTGGGCAAAATAGCCAAAGGTGCAGCCGAAGACGGGCCTTTGGctc                        >  1:1242541/1‑68 (MQ=255)
                                            ttAGCTCCGATAACGGCAGGGTGGGCAAAATAGCCAATGGGCCCGCCGAAGACGGGGCCTGGCCCCtt                      >  1:3875946/1‑68 (MQ=255)
                                            ttAGCTCCGATAACGGCAGGGGGAGCGAACTAACCAATGGTCCCGCCCCAGTCGGGCCCTGGCctctc                      >  1:686099/1‑67 (MQ=255)
                                               gCTCCGATAACGGCAGGGTGGGGAAACTAGCCAATGGGCCAGCCGGAGGCGGGGCTTGGGCTCttttt                   >  1:2279777/1‑68 (MQ=255)
                                               gCTCCGATAACGGCAGGGTGGGCAAAATAGCCAAAGGGCCCGCCGAAGACGGGCCCTGGCCCCttttt                   >  1:3336851/1‑68 (MQ=255)
                                               gCTCCGATAACGGCAGGGTGGGCAAAAAAGCCCAAAGTCGAGCCGCAGTAGGGCCCTGGCGGCttttt                   >  1:564230/1‑68 (MQ=255)
                                                cTCCGATAACGGCAGGGTGAGGAAAATAGCCAATGGTCCCGCCCAGGTGCGGCCTTTGGCCCtttttt                  >  1:1149532/1‑68 (MQ=255)
                                                           cgcggTTTAACAATCTACCCAATGGCCCAGCCAAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCt       <  1:1800535/65‑1 (MQ=255)
                                                               gTTGAGCAAACTAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTCCAc   >  1:3454492/1‑68 (MQ=255)
                                                                ttGAGCAAACTAGCCAATGGGCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTCCCCg  >  1:880341/1‑68 (MQ=255)
                                                                |                                                                   
ATTTTGCCTTTCTGTACTCAGCCACGCTGCTATGTAATCTGGGGTTAGCTCCGATAACGGCAGGTTGAGCAAACTAGCCAATGGTCCAGCCGAAGTCGGGCCTTGGCCTCTTTTTTTACTTTCGCCTCCACG  >  NC_000913/296440‑296571

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