Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,987,499	C→A	intergenic (+71/+8)	ftnB → / ← yecJ	putative ferritin‑like protein/DUF2766 domain‑containing protein YecJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,987,499	0	C	A	100.0%	56.2 / NA	21	intergenic (+71/+8)	ftnB/yecJ	putative ferritin‑like protein/DUF2766 domain‑containing protein YecJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (9/12);  total (9/12)

AGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCA  >  NC_000913/1987304‑1987707                                                                                                                                                                                                    |                                                                                                                                                                                                                 agcagcATGATGGTCTGTTGCTGCAAACCAGTCTTGATGAAGTGCGCAGTGCGAAACGTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATTTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCAGTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGGTTGAAGCAGAGGGTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGAt                                                                                                                                                           <  1:287855/251‑1 (MQ=255)                   tgctgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGGGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGccc                                                                                                                                          <  1:455810/250‑1 (MQ=255)                     ctgcAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGGGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACAAGCTGCCTTGATAATCATAGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGGTTTTTTAATGGGATTCACCCTGTGGGGAAAACTTGAGTTCAATAAGCGCGATGGCGTTTTGGATTGCACGc                                                                                                                                        >  1:345053/1‑250 (MQ=255)                                             cgcAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAAAGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGcggcggcgg                                                                                                               <  1:2097/251‑1 (MQ=255)                                                      aaaCTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTCTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGtt                                                                                                        >  1:498551/1‑249 (MQ=255)                                                                        ccTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACCCCAGCTGCATTGATCATCATCGGCGCTAACGCATTGCGCCGATGACGGTTTTGAGAACCCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCg                                                                                                                                         >  1:251933/1‑198 (MQ=255)                                                                                gaccgaccAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTg                                                                            <  1:458065/251‑1 (MQ=255)                                                                                                        cGTGTCACACCAGCTGCATGGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGGTTTGAGAAACCGCTGCCTCAGCTGTTTGAAGCAGATGTTTTTTTAATGTGATTCACCCTGTGGGGTAAACTTGAGTTCAAGAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCTGGATGGTTATTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCaaag                                                    <  1:462028/251‑3 (MQ=255)                                                                                                         gtgtCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCAc                                                                                                                                                                                             >  1:387364/1‑113 (MQ=255)                                                                                                                               catcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAACCCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAACCTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCTCCTGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATc                             >  1:8542/1‑251 (MQ=255)                                                                                                                                          cTAATGCATTGCGCCGATGAAGGTTTTGAGCAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTg                                                                                                                                             >  1:197987/1‑128 (MQ=255)                                                                                                                                          cTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTATTGGGATTCACCCTGTGGGGTAAACTTg                                                                                                                                                                            >  1:277178/1‑97 (MQ=255)                                                                                                                                          cTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCCCCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGAGATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTAc                                         <  1:493447/228‑1 (MQ=255)                                                                                                                                                                           ccGCTNCATCANCNNTTTGAAGCAGAGGTTTTTTNAATGGGATTCACCCTNTGGGNTAAACNTGNgn                                                                                                                                                                         >  2:191652/1‑66 (MQ=255)                                                                                                                                                                           ccGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGAGGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACAGTTTTTTACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGAt              <  1:191652/222‑1 (MQ=255)                                                                                                                                                                            cGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTNGGANGGCCNGCNTGGTGACCGGGTCTGNGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGNNANNNTNTCACGAAc                                                           <  2:362437/176‑1 (MQ=255)                                                                                                                                                                            cGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGAGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAAc                                                           >  1:362437/1‑176 (MQ=255)                                                                                                                                                                                  cTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGAGTCACCCTGTGGGGTAAACTTGAGTTCAAGAAGCGCGATGGCTGTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGGAAAGTCGATATGCTTCAGCTGACTGTACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAAtttgtttn                <  1:43300/213‑2 (MQ=255)                                                                                                                                                                                  cNCANCTGTTNGAAGNAGAGGTTTTTTTAATNGGATTCACCCTGTGGGGTTNNCNTGAGTNCAATa                                                                                                                                                                   <  2:186399/66‑1 (MQ=255)                                                                                                                                                                                   tCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGAGTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGgca  <  1:215478/226‑1 (MQ=255)                                                                                                                                                                                             gAAGCAGAGGTTTTTTTAAGGGGATTCACCCGGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTg                                                                                                                                                  <  1:245658/72‑1 (MQ=255)                                                                                                                                                                                                    |                                                                                                                                                                                                                 AGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCA  >  NC_000913/1987304‑1987707

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