Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,114,419	G→A	A18A (GCC→GCT)	fpr ←	ferredoxin‑NADP reductase; flavodoxin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,114,419	0	G	A	100.0%	64.5 / NA	23	A18A (GCC→GCT)	fpr	ferredoxin‑NADP reductase; flavodoxin reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/13);  total (10/13)

TTGGTAAATTGCCCGGCGGTAAACGGAAGCACGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCT  >  NC_000913/4114360‑4114483                                                            |                                                                 ttGGTTATTTGCCCGGCGGTATACGTTAGCACGGGGTCGTGAACGGTGAGACTAAACTGAGCGTCGGt                                                          <  2:207691/68‑1 (MQ=255) ttGGTAAATTGCCCGGCGGTAAACGGAAGCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGt                                                          >  2:653492/1‑68 (MQ=255) ttGGTAAATTGCCCGGCGGTAAACGGAAGCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGt                                                          <  2:379180/68‑1 (MQ=255) ttGGTAAATTGCCCGGCGGTAAACGGAAGCACGGGGGCGTGAACGGCGAGACTAAACAGAGCGTCGGt                                                          <  2:374596/68‑1 (MQ=255)      aaaTTGCCCGGCGGTAAACGGAAGCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGTCCAGt                                                     <  1:505227/68‑1 (MQ=255)       aaTTGCCCGGCGGTAAACGGAAGCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGTCCAGtt                                                    >  1:213752/1‑68 (MQ=255)           gCCCGGCGGTAAACGGAAGCACGGGGGCGTGAACGGTGAGACTACACAGAGCGTCGGTCCAGTTCCGc                                                <  1:170626/68‑1 (MQ=255)                       aCGGAAGCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACt                                    >  2:90508/1‑68 (MQ=255)                       aCGGAAGCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACt                                    <  1:620036/68‑1 (MQ=255)                          gAAGCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGTCCAGTTCTGCGCTTTAGTGACTTTg                                 >  1:263941/1‑68 (MQ=255)                            aGCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACt                                    <  1:2889/63‑1 (MQ=255)                            aGCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACt                                    >  2:2889/1‑63 (MQ=255)                             gCACGGGGGCGTGAACGGTGAGACTAAACAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCt                              <  2:293806/68‑1 (MQ=255)                                    ggCGTGAACGGTGAGACTAAACAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTAccc                       >  2:211367/1‑68 (MQ=255)                                           aCGGTGAGACTAAACAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCTATCAGc                >  2:489221/1‑68 (MQ=255)                                             ggTGAGACTAAACAGAGCGTCGGTCCAGTTCTGCACTGTAGTGGCTTTGCCTGTTGCCCAATCAGccg              >  1:710152/1‑67 (MQ=255)                                                 agaCTAAACAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGtt          <  2:313622/68‑1 (MQ=255)                                                 agaCTAAACAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGtt          <  1:211367/68‑1 (MQ=255)                                                 agaCTAAACAGAGCGGCGGTCCAGTACTGCACTTTAGGGACTTTGCCTGTTACCCAATCAGCCATGtt          <  2:719383/68‑1 (MQ=255)                                                      ataCAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTtct     <  2:436413/66‑1 (MQ=255)                                                      aaaCTGAGCGTCGGCCCAGTTCTGCACATTAGTGACTTTGCCTGTAACCCAATCAGCCATGTTTTtct     <  2:123111/68‑1 (MQ=255)                                                         cAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCt  >  1:457879/1‑68 (MQ=255)                                                         cAGAGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTATTTCTCCt  >  2:62944/1‑68 (MQ=255)                                                            |                                                                 TTGGTAAATTGCCCGGCGGTAAACGGAAGCACGGGGGCGTGAACGGTGAGACTAAACAGGGCGTCGGTCCAGTTCTGCACTTTAGTGACTTTGCCTGTTACCCAATCAGCCATGTTTTTCTCCT  >  NC_000913/4114360‑4114483

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