Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,437,930	(A)6→5	coding (224/939 nt)	ytfI →	protein YtfI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,437,925	0	A	.	100.0%	78.6 / NA	19	coding (219/939 nt)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (8/11);  total (8/11)

GATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAA  >  NC_000913/4437848‑4438001                                                                              |                                                                             gATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCaaa                                                                            >  2:895062/1‑80 (MQ=255)          ttAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGaaaa                                                                   <  2:895045/79‑2 (MQ=255)                 tGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACg                                                           <  2:895046/80‑1 (MQ=255)                  gAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGc                                                          >  2:895069/1‑80 (MQ=255)                    gggAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCAt                                                       <  2:895047/81‑1 (MQ=255)                      gAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATg                                                      >  2:895070/1‑80 (MQ=255)                           cTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATaa                                                 <  2:895048/80‑1 (MQ=255)                            tttCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGt                                               >  2:895066/1‑81 (MQ=255)                                  tGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTa                                         >  2:895068/1‑81 (MQ=255)                                            aGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTc                                 <  1:895049/79‑1 (MQ=255)                                            aGATTGATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAg                               >  2:895067/1‑81 (MQ=255)                                                 gATAAAAGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCtt                           >  1:895064/1‑80 (MQ=255)                                                     aaaGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAg                       <  1:895050/80‑1 (MQ=255)                                                      aaGTGGTTATCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGAt                     <  2:895051/81‑1 (MQ=255)                                                          ggTTATCAGGTGCTTTCATC‑AAAAATGAAGCTGTTAACGCCATGGATGTGTTTCTTATCAGTTTCAGCCTTAGGGATAACg                 <  2:895052/81‑1 (MQ=255)                                                               tCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGcc            <  2:895054/81‑1 (MQ=255)                                                               tCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGcc            <  2:895053/81‑1 (MQ=255)                                                               tCAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGcc            >  1:895065/1‑81 (MQ=255)                                                                cAGGTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGcc            >  1:895071/1‑80 (MQ=255)                                                                   gTGCTTTCATC‑AAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAat        <  2:895055/81‑1 (MQ=255)                                                                         tcatc‑aAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACaa  >  1:895074/1‑81 (MQ=255)                                                                              |                                                                             GATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAA  >  NC_000913/4437848‑4438001

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