Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	278,518	Δ1 bp	coding (285/1047 nt)	afuC ←	CP4‑6 prophage; ABC transporter ATP‑binding protein AfuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	278,516	0	C	.	100.0%	117.3 / NA	28	coding (287/1047 nt)	afuC	CP4‑6 prophage; ABC transporter ATP‑binding protein AfuC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (11/17);  total (11/17)

TCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGAT  >  NC_000913/278436‑278588                                                                                 |                                                                         tCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATtctctc                                                                          <  2:58503/81‑1 (MQ=255) tCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATtctctc                                                                          <  2:58504/81‑1 (MQ=255)   gCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTcc                                                                         <  2:58505/80‑1 (MQ=255)      aaCGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTccagcg                                                                     <  2:58506/81‑5 (MQ=255)       aCGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGa                                                                   <  1:58507/81‑1 (MQ=255)             ctTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATg                                                             >  1:58544/1‑81 (MQ=255)               ttGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCg                                                           >  1:58555/1‑81 (MQ=255)               ttGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCg                                                           >  1:58556/1‑81 (MQ=255)                  aCGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCggg                                                         <  2:58508/80‑1 (MQ=255)                   cgcgGGCTTTCGGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGaa                                                       <  1:58509/81‑1 (MQ=255)                      gggCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAg                                                    >  1:58548/1‑81 (MQ=255)                      gggCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAg                                                    <  1:58510/81‑1 (MQ=255)                           ttCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCa                                                >  1:58551/1‑80 (MQ=255)                            tCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCAt                                               <  1:58511/80‑1 (MQ=255)                                ctctGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGAc                                          >  1:58552/1‑81 (MQ=255)                                     cgcgcgGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGaaa                                     <  1:58513/81‑1 (MQ=255)                                     cgcgcgGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGa                                       <  1:58512/79‑1 (MQ=255)                                      gcgcgGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGaaa                                     >  2:58566/1‑80 (MQ=255)                                        gcgGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAacac                                  <  2:58514/81‑1 (MQ=255)                                          ggTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCa                                <  2:58515/81‑1 (MQ=255)                                          ggTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCa                                >  2:58559/1‑81 (MQ=255)                                          ggTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCa                                <  1:58516/81‑1 (MQ=255)                                           gTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCAt                               <  2:58517/81‑1 (MQ=255)                                                  gAGCATTTTCAGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGata                        <  1:58518/81‑1 (MQ=255)                                                            aGGCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCgctgct              <  2:58519/81‑1 (MQ=255)                                                             ggCCATAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGctgctg             <  2:58520/81‑1 (MQ=255)                                                                  tAACCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATag        <  2:58521/81‑1 (MQ=255)                                                                    aCCGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATagag      <  1:58522/81‑1 (MQ=255)                                                                     ccGACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGc     >  1:58557/1‑81 (MQ=255)                                                                       gACATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGa   >  2:58554/1‑81 (MQ=255)                                                                        aCATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGAt  <  1:58523/81‑1 (MQ=255)                                                                        aCATTCTCT‑CCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGAt  >  2:58563/1‑81 (MQ=255)                                                                                 |                                                                         TCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGAT  >  NC_000913/278436‑278588

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