Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 267,012 A→G intergenic (‑459/+172) yafY ← / ← yafZ CP4‑6 prophage; inner membrane lipoprotein YafY/CP4‑6 prophage; conserved protein YafZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913267,0120AG100.0% 89.1 / NA 26intergenic (‑459/+172)yafY/yafZCP4‑6 prophage; inner membrane lipoprotein YafY/CP4‑6 prophage; conserved protein YafZ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/15);  total (11/15)

CAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCG  >  NC_000913/266935‑267087
                                                                             |                                                                           
cagccGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCtt                                                                          >  2:56285/1‑81 (MQ=255)
cagccGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCtt                                                                          <  1:56235/81‑1 (MQ=255)
 agccGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCttt                                                                         <  1:56236/81‑1 (MQ=255)
    cGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTa                                                                        >  1:56268/1‑79 (MQ=255)
     ggTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAag                                                                     <  2:56237/81‑1 (MQ=255)
               gcgTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCg                                                            >  1:56279/1‑80 (MQ=255)
                 gTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGa                                                          >  2:56280/1‑80 (MQ=255)
                       ccTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGt                                                    <  1:56238/80‑1 (MQ=255)
                       ccTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGt                                                    <  2:56239/80‑1 (MQ=255)
                          gtgtCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATg                                                 >  1:56277/1‑80 (MQ=255)
                                  agaCGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGAcaccac                                         <  2:56240/80‑1 (MQ=255)
                                      gTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGaaa                                    >  1:56273/1‑81 (MQ=255)
                                       ttCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAAc                                   <  1:56241/81‑1 (MQ=255)
                                         ctcCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTa                                 >  2:56281/1‑81 (MQ=255)
                                             gTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTAtttt                             >  2:56283/1‑81 (MQ=255)
                                               gAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTAttttt                            >  2:56288/1‑80 (MQ=255)
                                                aaaGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTc                           <  1:56242/80‑1 (MQ=255)
                                                     tAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCt                     >  1:56286/1‑81 (MQ=255)
                                                       actgactgAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGc                   <  1:56244/81‑1 (MQ=255)
                                                       actgactgAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGc                   <  1:56243/81‑1 (MQ=255)
                                                           actgaAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCaa                 <  2:56245/79‑1 (MQ=255)
                                                               aaaaTGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGc             >  2:56276/1‑79 (MQ=255)
                                                               aaaaTGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCa           <  1:56246/81‑1 (MQ=255)
                                                                   tgatgTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAgg       <  2:56247/81‑1 (MQ=255)
                                                                     atgTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCa     <  2:56248/81‑1 (MQ=255)
                                                                        tgATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCg  <  2:56249/81‑1 (MQ=255)
                                                                             |                                                                           
CAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCG  >  NC_000913/266935‑267087

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