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_000913	267,012	0	A	G	100.0%	118.0 / NA	35	intergenic (‑459/+172)	yafY/yafZ	CP4‑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/24);  total (11/24)

CCAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGG  >  NC_000913/266934‑267082                                                                               |                                                                       ccagccGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCt                                                                      <  2:61955/81‑1 (MQ=255) ccagccGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCt                                                                      >  2:62002/1‑81 (MQ=255)      cGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTaaa                                                                 <  1:61956/81‑1 (MQ=255)      cGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTa                                                                   >  2:62007/1‑79 (MQ=255)        gTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTaaa                                                                 >  1:61999/1‑79 (MQ=255)               cggcgTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTc                                                        <  1:61957/81‑1 (MQ=255)                    ttCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACa                                                   <  1:61958/81‑1 (MQ=255)                         ccTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGt                                               <  1:61959/80‑1 (MQ=255)                         ccTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGt                                               <  2:61960/80‑1 (MQ=255)                         ccTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTa                                              >  1:62005/1‑81 (MQ=255)                          cTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTAt                                             <  1:61961/81‑1 (MQ=255)                                 aTGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGacac                                       <  2:61962/80‑1 (MQ=255)                                  tGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGAcacca                                     <  1:61964/81‑1 (MQ=255)                                  tGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGAcacca                                     <  2:61963/81‑1 (MQ=255)                                    agaCGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGAcaccac                                    <  2:61966/80‑1 (MQ=255)                                    agaCGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACcaca                                   <  1:61965/81‑1 (MQ=255)                                     gaCGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGAcacca                                     >  2:62008/1‑78 (MQ=255)                                     gaCGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGAcacca                                     >  2:62014/1‑78 (MQ=255)                                       cGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGaa                                <  1:61967/81‑1 (MQ=255)                                       cGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGaa                                <  1:61968/81‑1 (MQ=255)                                         ttCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAAc                              <  1:61969/81‑1 (MQ=255)                                           ctcCGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACt                             <  1:61970/80‑1 (MQ=255)                                              cGTGAAAGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTAttt                         <  1:61971/81‑1 (MQ=255)                                                  aaaGTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTc                      <  2:61972/80‑1 (MQ=255)                                                     gTTAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCAttt                  >  2:62012/1‑81 (MQ=255)                                                      ttAACTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTc                 <  2:61973/81‑1 (MQ=255)                                                        aaCTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCt                <  2:61975/80‑1 (MQ=255)                                                        aaCTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCt                <  2:61974/80‑1 (MQ=255)                                                        aaCTGACTGAAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTg               >  2:62009/1‑81 (MQ=255)                                                             actgaAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCaa            >  1:62006/1‑79 (MQ=255)                                                              ctgaAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAg         <  2:61976/81‑1 (MQ=255)                                                              ctgaAAATGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAg         >  1:62022/1‑81 (MQ=255)                                                                 aaaaTGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGc        >  1:62019/1‑79 (MQ=255)                                                                  aaaTGATGTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCag     <  1:61977/81‑1 (MQ=255)                                                                     tgatgTGATTGCTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAgg  <  2:61978/81‑1 (MQ=255)                                                                               |                                                                       CCAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGG  >  NC_000913/266934‑267082

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