Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009132,774,8510A→C100.0% 1.9 3K512Q (AAG→CAG) yfjWCP4‑57 prophage; putative inner membrane protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (1/0):  total (3/0)

GATGATATTATTAATTACAACAATAGCCCTTCCGTACTACTCAAGACTGATTTTGCGTTCTATAAAAAAACATACCAGAAGCTTGATAAAATATATGATGACATAAGCAATGGAAAATTATCTTCACTGAGAGCAACAGGTATCTCCCAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAA  >  NC_000913/2774563‑2775002
                                                                                                                                                                                                                                                                                                |                                                                                                                                                       
gatgatATTATTAATTACAACAATAGCCCTTCCGTACTACTCAAGACTGATTTGGCGTTCTATAAAAAAACATACCAGAAGCTTGATAAAATATATGATGACATAAGCAATGGAAAATTATCTTCACTGAGAGCAACAGGTATCTCCCAGTTTAGTATCAACGGGACACATTTGTCGTTAAGGCCAAAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAACAATCAAAATATATGGCACAGCCGATCTTGTTTTTGTTTATAAACAGATAGTGAAc                                                                                                                                              >  1:508712/1‑300 (MQ=255)
                   acaaTAGCCCTTCCGTACTACTCAAGACTGATTTTGCGTTCTATAAAACAACATACCAGAAGCTTGATCACATCTATGATGACATAAGCACTGGACCCTTATCTTCACTGAGCGCAACAGGTATCTCCCCGTTTAGTATCAACGGGCAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTCATAAAAGCAAACAACTCAAAATATCTGGCACAGCAGATCTTGTTTTTGTTGATAACCAGTTAATGAAccctgtcccttctccacca                                                                                                                            >  1:15380/1‑282 (MQ=255)
                                                                                                                                           gTATCTCCCAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAGATTTTTGTTTATTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTAATTCATTCTTACTTTAAGAATGATGtaaa  >  2:765057/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                       
GATGATATTATTAATTACAACAATAGCCCTTCCGTACTACTCAAGACTGATTTTGCGTTCTATAAAAAAACATACCAGAAGCTTGATAAAATATATGATGACATAAGCAATGGAAAATTATCTTCACTGAGAGCAACAGGTATCTCCCAGTTTAGTATCAACGGGAAACATTTGTCGTTAAGGCCAGAATCGGAGATTATTATTTCTGAAGGGAGTTTGTATGGACTAGTTAATAAAAGCAAAAAAATCAAAATATATGGCACAGCAGATCTTGTTTTTGTTGATAACAAGATAATGAACCTTCGTAAAATAACTTATTTGCAATCTAAGCTAGAAATTTTTGGTTCTTCTATTATGGATATATTAAAGTATATATTTGGTTTAGGTCTGCTAGCAATTTCTATAAAATTCATTCATTCTTACTTTAAGAATGATGTGAA  >  NC_000913/2774563‑2775002

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.