Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	266,547	G→T	R3S (CGT→AGT)	yafY ←	CP4‑6 prophage; inner membrane lipoprotein YafY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	266,547	0	G	T	100.0%	117.2 / NA	33	R3S (CGT→AGT)	yafY	CP4‑6 prophage; inner membrane lipoprotein YafY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/17);  total (16/17)

TCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCA  >  NC_000913/266474‑266623                                                                          |                                                                             tcatcaCTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCatt                                                                       <  1:56140/81‑1 (MQ=255) tcatcaCTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCat                                                                        >  1:56185/1‑80 (MQ=255)  catcaCTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCatta                                                                      >  2:56178/1‑81 (MQ=255)     caCTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTAttt                                                                   >  1:56182/1‑81 (MQ=255)      aCTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTAttt                                                                   <  2:56141/80‑1 (MQ=255)       cTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTAttttt                                                                 <  1:56142/81‑1 (MQ=255)        ttCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTAtttttt                                                                >  2:56179/1‑81 (MQ=255)          cTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTc                                                              >  2:56187/1‑81 (MQ=255)             tCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTc                                                              <  2:56143/78‑1 (MQ=255)             tCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTc                                                              <  1:56144/78‑1 (MQ=255)               aTCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGaa                                                         <  1:56145/81‑1 (MQ=255)               aTCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGa                                                          >  2:56177/1‑80 (MQ=255)                    aTGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAg                                                    <  2:56146/81‑1 (MQ=255)                    aTGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAg                                                    <  2:56147/81‑1 (MQ=255)                             gaAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCAc                                           <  1:56148/81‑1 (MQ=255)                                  agagTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTgg                                      >  2:56195/1‑81 (MQ=255)                                  agagTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTgg                                      <  2:56149/81‑1 (MQ=255)                                            cAACCAGCGCCAGCATAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACTAACCACTTTGGCAAAGACATa                            <  2:56150/81‑1 (MQ=14)                                             aaCCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATa                            >  1:56183/1‑80 (MQ=255)                                              aCCAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATaaa                          <  2:56151/81‑1 (MQ=255)                                                cAGCGCCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATg                        >  1:56196/1‑81 (MQ=255)                                                    gcCAGCAGAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCa                    >  1:56193/1‑81 (MQ=255)                                                     ccAGCAGAGGCAGTGTTCTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCAc                   <  1:56152/81‑1 (MQ=255)                                                       agcagAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACg                  >  2:56184/1‑80 (MQ=255)                                                       agcagAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGa                 <  2:56153/81‑1 (MQ=255)                                                         cagAGGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGaa                >  2:56188/1‑80 (MQ=255)                                                           gagGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGaaaa              >  1:56186/1‑80 (MQ=255)                                                            agGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATa            <  2:56154/81‑1 (MQ=255)                                                            agGCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATa            >  1:56191/1‑81 (MQ=255)                                                             ggCAGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATAc           >  2:56205/1‑81 (MQ=255)                                                                aGTGTTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCa        <  1:56155/81‑1 (MQ=255)                                                                 gtgtTTTACTTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCAc       <  1:56156/81‑1 (MQ=255)                                                                      ttACTTTTCATTATTTTTTTCCTGACATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCa  >  2:56194/1‑81 (MQ=255)                                                                          |                                                                             TCATCACTTCTGTCATCGCATGCGATGAGAAACAGAGTGGTGGCAACCAGCGCCAGCAGAGGCAGTGTTTTACGTTTCATTATTTTTTTCCTGAAATCAGACGAACCACTTTGGCAAAGACATAAATGCCCACGAAAATACCCACCGGCA  >  NC_000913/266474‑266623

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: ‑