Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	267,042	A→G	intergenic (‑489/+142)	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,042	0	A	G	100.0%	118.2 / NA	35	intergenic (‑489/+142)	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 (19/16);  total (19/16)

CAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAG  >  NC_000913/266917‑267170                                                                                                                              |                                                                                                                                 cAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTa                                                                                                                    >  1:160719/1‑140 (MQ=255)      acagcctgacagccAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTAttttt                                                                                                               >  1:524330/1‑140 (MQ=255)          cctgacagccAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCAtt                                                                                                           >  1:826567/1‑140 (MQ=255)                    agccGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTTAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAg                                                                                                 >  1:901897/1‑140 (MQ=255)                        ggTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCag                                                                                             <  2:815862/140‑1 (MQ=255)                                       cTGCCTGTATCATGAGTCGTTCTCCGTGAAAGTTACCTGGCTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTa                                                                              <  2:482941/140‑1 (MQ=18)                                        tGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTAc                                                                             >  2:844592/1‑140 (MQ=255)                                                   tGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCt                                                                  <  2:637681/140‑1 (MQ=255)                                                    gagaCGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTc                                                                 <  2:123395/140‑1 (MQ=255)                                                          ttCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGgtgag                                                           >  1:203130/1‑140 (MQ=255)                                                              ccGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGcc                                                                        <  1:137297/123‑1 (MQ=255)                                                              ccGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGcc                                                                        >  2:137297/1‑123 (MQ=255)                                                              ccGTGAAAGTTAACTCACTGAACATGATGAGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGtgagtgag                                                       <  1:492945/140‑1 (MQ=255)                                                                 tGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCGCCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGAc                                                    >  1:209453/1‑140 (MQ=255)                                                                           ctgactgaAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACg                                          >  1:97026/1‑140 (MQ=255)                                                                             gactgaAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGt                                        <  2:203130/140‑1 (MQ=255)                                                                                  aaaaTGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGAt                                   <  2:143546/140‑1 (MQ=255)                                                                                   aaaTGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCGCCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATg                                  <  2:209453/140‑1 (MQ=255)                                                                                   aaaTGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATg                                  <  1:844592/140‑1 (MQ=255)                                                                                     atgatgTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCg                                <  2:160719/140‑1 (MQ=255)                                                                                       gatgTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGa                              <  2:697448/140‑1 (MQ=255)                                                                                        atgTGATTACTTTAAATAGGGGCTCGGACAGGGAAAGGGTACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGAc                             <  2:444621/140‑1 (MQ=37)                                                                                         tgtgATTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTTGGTGAGTAAGGCCCATCAGCACGGTGCGATGCGGACa                            >  1:381247/1‑140 (MQ=25)                                                                                             aTTACTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGACAGCACGGTGCGATGCGGACAggtg                        >  1:465999/1‑140 (MQ=255)                                                                                                 cTTTAAAGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGgtgt                    <  1:878490/140‑1 (MQ=255)                                                                                                     aaaGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGa                >  2:1315/1‑140 (MQ=255)                                                                                                      aaGAGGGGGTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGat               <  2:902674/140‑1 (MQ=255)                                                                                                           ggggTCGGACAGGGTATGGGCACCACAGAAACTGTTGTTCATTTCTGCAAAAGCCAGAGGGATTTGTGGTTACAGGCCGTCCCGCGGTTATTAAGTACCTGCAGCACGGGGCGAGGCGGAGATGTGCGGGCCGATATTaa          >  2:639134/1‑140 (MQ=255)                                                                                                           ggggTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTaa          >  1:221200/1‑140 (MQ=255)                                                                                                           ggggTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTaa          >  2:987030/1‑140 (MQ=255)                                                                                                            gggTCGGACAGGGTATGGGCACCACGGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAAc         >  2:695359/1‑140 (MQ=255)                                                                                                            gggTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTGAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAAc         >  2:327360/1‑140 (MQ=255)                                                                                                            gggTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAAc         >  1:589126/1‑140 (MQ=255)                                                                                                            gggTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAAc         >  2:895338/1‑140 (MQ=255)                                                                                                            gggTCGGACAGGGTATGGGCACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAAc         >  2:936665/1‑140 (MQ=255)                                                                                                                   aCAGGGTATGGGCACCTCAGAATCTATTTTTCAGTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAagcag  <  2:465999/140‑1 (MQ=255)                                                                                                                              |                                                                                                                                 CAGTGACAGCCTGACAGCCAGCCGGTCATGGCGGCGTTCTGCCTGTGTCATGAGACGTTCTCCGTGAAAGTTAACTGACTGAAAATGATGTGATTACTTTAAAGAGGGGGTCGGACAGGGTATGGACACCACAGAAACTATTTTTCATTTCTGCAAAAGCCAGAGGCAGCGGGGGTTACAGGCCTTCCCTCGGTGAGTGAGGACCGTCAGCACGGTGCGATGCGGACAGGTGGTGTCCGATATTAACAAAGCAG  >  NC_000913/266917‑267170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑