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%	71.4 / NA	23	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 (10/13);  total (10/13)

CCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTG  >  NC_000913/1987271‑1987734                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            ccGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCAcc                                                                                                                                                                                                                       >  1:346896/1‑251 (MQ=255)              gcgcGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTaaa                                                                                                                                                                                                           <  1:266217/250‑1 (MQ=255)                            aaGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAg                                                                                                                                                                                            <  1:298555/251‑1 (MQ=255)                                 cagcagCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCCGTGCGCAACTTGCGGGTCTGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGcgcg                                                                                                                                                                                        >  2:192323/1‑250 (MQ=255)                                                                     gatgaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTg                                                                                                                                                   <  1:347709/251‑1 (MQ=255)                                                                        gaAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTgcgg                                                                                                                                                <  1:367653/251‑1 (MQ=255)                                                                                  gTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAAGAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATg                                                                                                                                      <  1:278714/251‑1 (MQ=255)                                                                                    gCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATgg                                                                                                                                     <  2:346896/250‑1 (MQ=255)                                                                                                     gtgCCCTGTGCAGCCCGAGCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCTTTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCGCAGCTGTTTGAAGCAGAGGTTTTGTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTGTTGGATTGGCCGCATGTTGACCGGGTCTGCGGCGGCGGGATGGTGAGTAAAGTCGTTAttc                                                                                                                   <  2:435028/251‑1 (MQ=255)                                                                                                           tgtgCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTTAAGTCGATATTCTTCAGc                                                                                                             >  2:107306/1‑251 (MQ=255)                                                                                                                              gTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCAc                                                                                          >  1:308034/1‑251 (MQ=255)                                                                                                                                                    aGCTGCATTGATCGTCATCGGCGCTAATACATTACGCCGATGCAGGGTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTCAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAGCGGGCGCGAtttc                                                                    >  2:216432/1‑249 (MQ=255)                                                                                                                                                        gCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCg                                                                <  2:414969/251‑1 (MQ=255)                                                                                                                                                              atcatcatcGGCGCTAATGCATAGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACACTGTGGGGAAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGTTTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCCTCAGCTGACTGGACATTTTTTAACGCACTTCAACGGGCGCGATTACATCGAGAAc                                                           >  1:72864/1‑250 (MQ=255)                                                                                                                                                                catcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGata                                                                                                                      <  1:408124/189‑1 (MQ=255)                                                                                                                                                                catcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGata                                                                                                                      >  2:408124/1‑189 (MQ=255)                                                                                                                                                                catcatcGCCGCTAATGCATTGCGCCGATGCAGGTTTTGACCAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTACACTTGAGTTCACTAAGCGCGATGGCTTTTTGGATTGCCCGCATGGCGACCGGGTCTGCGGCGGCGGGATGGTTAGTAACGTAGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATc                                                        >  2:419864/1‑251 (MQ=255)                                                                                                                                                                 atcatcGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGTGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATcc                                                       <  1:253675/251‑1 (MQ=255)                                                                                                                                                                                        ccGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTg                                >  1:338182/1‑251 (MQ=255)                                                                                                                                                                                                    ttGAGAAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGGGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATGTGTTTGATAACCAACTGGCAGGCAACca                     <  2:35486/250‑1 (MQ=255)                                                                                                                                                                                                      gagaAACCGCTGCCTCATCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACcaca                   <  1:275648/250‑1 (MQ=255)                                                                                                                                                                                                       agaAACCGCTGCCTCAGCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCACCCACATc                 <  2:363346/251‑1 (MQ=255)                                                                                                                                                                                                                      aTCTGTTTGAAGCAGAGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTg  >  1:219121/1‑251 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            CCGTCAATTTCCTGCGCGATCTGGAAAAAGAACAGCAGCATGATGGTCTGTTGCTGCAAACCATTCTTGATGAAGTGCGCAGTGCGAAACTTGCGGGTATGTGCCCTGTGCAGACCGACCAACATGTTCTGAATGTCGTGTCACACCAGCTGCATTGATCATCATCGGCGCTAATGCATTGCGCCGATGAAGGTTTTGAGAAACCGCTGCCTCATCTGTTTGAAGCAGCGGTTTTTTTAATGGGATTCACCCTGTGGGGTAAACTTGAGTTCAATAAGCGCGATGGCTTTTTGGATTGCCCGCATGGTGACCGGGTCTGCGGCGGCGGGATGGTTAGTAAAGTCGATATTCTTCAGCTGACTGGACATTTTTTCACGAACTTCAACGGGCGCGATTACATCGAGAACATCCAGAATTTGTTTGATAACCAACTGGCAAGCAACCACATCAGAAACCAATTCCTG  >  NC_000913/1987271‑1987734

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