Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,113,851	G→A	H208Y (CAT→TAT)	fpr ←	flavodoxin/ferredoxin‑NADP(+) reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,851	0	G	A	100.0%	99.0 / NA	28	H208Y (CAT→TAT)	fpr	flavodoxin/ferredoxin‑NADP(+) reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/11);  total (17/11)

TGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTAGT  >  NC_000913/4113613‑4114095                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     tGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCt                                                                                                                                                                                                                                            <  2:126942/249‑1 (MQ=255)        ggTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCAt                                                                                                                                                                                                                                    <  1:342083/250‑1 (MQ=255)           catcatcaCTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGc                                                                                                                                                                                                                                >  2:321777/1‑251 (MQ=255)                                        gCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACt                                                                                                                                                                                                   <  1:271031/251‑1 (MQ=255)                                                atCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAAtt                                                                                                                                                                                            >  1:97503/1‑250 (MQ=255)                                                              aaTAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTAtccgt                                                                                                                                                                             >  1:122366/1‑251 (MQ=255)                                                              aaTAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTAtccgt                                                                                                                                                                             <  1:39520/251‑1 (MQ=255)                                                                   gATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATCAGCAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGt                                                                                                                                                                        >  2:28794/1‑251 (MQ=255)                                                                      gAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTg                                                                                                                                                                       >  1:257762/1‑249 (MQ=255)                                                                                     cATCGCGCCGTTGATCGTTACGTAAACGCTTAGCAGTACTCCTCCGCTGGCATATGGCCCAGTCGGCGACGTAAATGTTTCGACATCTGACGGGTCTCTTTCAGCGACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCAACTAGCTGGTTTCTTTATTCATCGGAAGGACCAGCGTGCTTACCCATTCCCCACTTTCCATTAATGCCGTGAtccgtcccggggctgacccccttgccgt                                                                                                                                                      >  2:192958/1‑230 (MQ=255)                                                                                            ccGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGa                                                                                                                                               <  1:54341/251‑1 (MQ=255)                                                                                                  aTCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGacc                                                                                                                                         >  2:291484/1‑251 (MQ=255)                                                                                                   tCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTc                                                                                                                                                   >  2:77677/1‑240 (MQ=255)                                                                                                   tCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTc                                                                                                                                                   <  1:77677/240‑1 (MQ=255)                                                                                                                gCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATgcg                                                                                                                           >  1:162666/1‑251 (MQ=255)                                                                                                                               tCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAGCTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGt                                                                                                             >  2:269542/1‑250 (MQ=255)                                                                                                                                                              aaTGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCaaa                                                                             <  1:366157/251‑1 (MQ=255)                                                                                                                                                                      gccATCTGCCGGGTCTCATTCAGCAACTGGTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATAACATAGCTGGTTTCTTTATTCATCGGCAGGCCACTCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGGCCACCGCCTGAACGCGCCGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCCGTGGCAAATAACTTa                                                                      >  2:397377/3‑250 (MQ=255)                                                                                                                                                                             gCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTcggc                                                               <  2:86492/250‑1 (MQ=255)                                                                                                                                                                               cGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCggcggc                                                            >  1:327448/1‑251 (MQ=255)                                                                                                                                                                                           aGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTgcggc                                                >  1:52831/1‑251 (MQ=255)                                                                                                                                                                                               aCTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGc                                            <  2:42587/251‑1 (MQ=255)                                                                                                                                                                                                               cACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGAtttt                            >  1:337164/1‑251 (MQ=255)                                                                                                                                                                                                                               ccGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCGAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATc            <  2:25990/251‑1 (MQ=255)                                                                                                                                                                                                                                   acaGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTa        <  2:134762/251‑1 (MQ=255)                                                                                                                                                                                                                                      gCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCt     >  2:355215/1‑251 (MQ=255)                                                                                                                                                                                                                                       cATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTACTGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCt     >  1:33958/1‑250 (MQ=255)                                                                                                                                                                                                                                         tCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTAGt  >  1:7430/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     TGTCTGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTAGT  >  NC_000913/4113613‑4114095

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