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%	110.4 / NA	32	H208Y (CAT→TAT)	fpr	flavodoxin/ferredoxin‑NADP(+) reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/13);  total (19/13)

TGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTA  >  NC_000913/4113617‑4114093                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   tgtgGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCCGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGAGCCGTTTATCGATAAGCAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTc                                                                                                                                                                                                                                    >  1:226794/1‑251 (MQ=255)    ggTCATCATCACTCATGTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATc                                                                                                                                                                                                                                 <  2:229423/251‑1 (MQ=255)                 cATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAAt                                                                                                                                                                                                                      <  1:212385/249‑1 (MQ=255)                  aTTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGt                                                                                                                                                                                                                   <  2:78214/251‑1 (MQ=255)                                                      gatCAATCAGATGAGTGCGGATTTGCTCCTCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCCATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTAt                                                                                                                                                                               >  2:28613/1‑251 (MQ=255)                                                                aTGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTg                                                                                                                                                                     >  1:139243/1‑251 (MQ=255)                                                                              gCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTGCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTg                                                                                                                                                        <  1:277720/250‑1 (MQ=255)                                                                               cTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGcc                                                                                                                                                      >  1:74593/1‑251 (MQ=255)                                                                                tCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCAGTATCCCTCCAGTGAGCGACCCCGCttc                                                                                                                                                       >  2:155658/1‑247 (MQ=255)                                                                                               ccGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTGCCAGTTCCCCACTTTCAAGGAATGCCGGTATCCGGCCGGTGAGCGACCCCGCTGCCGTGTCCCGACGGacca                                                                                                                                      <  2:132795/250‑1 (MQ=255)                                                                                                    aaGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTc                                                                                                                                 >  1:55891/1‑251 (MQ=255)                                                                                                                aCCAGTAATTCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCCACCGGTGAGCGACCCCCCTGCCGTTTCCCGACTGACCACCGTCTGAAAgcgccgt                                                                                                                     >  2:226409/1‑248 (MQ=255)                                                                                                                 ccAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGGTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCGGTGGATTGCCGCACAGCATCACATAGCTGGGTTCTTTATTCATCGGCAGGCCAATCGTGCGTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGGTGCCCGACTGTCCACCGTCTGAATGCGCAGtt                                                                                                                    <  2:82228/251‑1 (MQ=255)                                                                                                                     tAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTcc                                                                                                                >  2:165590/1‑251 (MQ=255)                                                                                                                          cTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGCCCACCGTCTGCATGCGCCGTTTTCCTTc                                                                                                             >  1:81579/1‑249 (MQ=255)                                                                                                                          cTCCGCTGGCATAGGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTGTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCTCCGTCTGAATGCGCAGTTTGCCTTc                                                                                                             <  2:81579/249‑1 (MQ=255)                                                                                                                          atCCGCTGTAATATGGCCCGGTCGGCGACGTAACTGTTTAGTCATCTGCCGGGTCTCTTTCAGCAACTGCTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGGGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCAGCTGCCGTTTCCCGACTGCCCACCGTCTGACTGCGCAGTTTTCCTTCGt                                                                                                           >  2:98466/2‑251 (MQ=255)                                                                                                                                                cggagacgtctaTGTGTCGTCATATGCCGGCTCTCTTTCAGCAACTGATGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACCTAGCTGGATTCTTTATTCATCGGCCGGCCAATCGTGCTTTCCAGTTCCTCACCTTCAATTAATGCCGGTATCCGTCCGGGGAGCGACCCCGCTGCCGCTTCTCGACTGACCACCGTCTGCATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGctt                                                                                     >  2:248588/12‑249 (MQ=255)                                                                                                                                                cGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCAt                                                                                     >  1:178546/1‑251 (MQ=255)                                                                                                                                                       gTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGc                                                                              >  2:143621/1‑251 (MQ=255)                                                                                                                                                              tttCGTCATCTGCCTGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCACCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGTGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGTCAAATAg                                                                        <  1:164649/250‑1 (MQ=255)                                                                                                                                                                   tCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTa                                                                    >  1:164760/1‑249 (MQ=255)                                                                                                                                                                     aTCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTc                                                                >  1:177772/1‑251 (MQ=255)                                                                                                                                                                                                         cgcACCATCTGTGGATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGAtt                            >  1:211228/1‑251 (MQ=255)                                                                                                                                                                                                                 cTGTGGATTGCAGCACAGCATCACATACCTGGTTTCTTTATTCCTCGGCAGGCCAATCGGGCCTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCAACCCCGCTGCCGTTTCCCGACTGACCAACGCCTGAATGCGCAGTTTTCCTTCGAAGCGTTTTTCCAGTTCCGGCTTAAGTGGCAAATACATTAAGACGGCGGCATCACGTGCGGCGTGCACCAGGACCAGATTTTTGAAgg                    >  2:175887/1‑250 (MQ=255)                                                                                                                                                                                                                     ggATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCCGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGGGCACCAGGCCCAGATTTTTGAAGCGATc                >  2:144040/1‑251 (MQ=255)                                                                                                                                                                                                                      gATTGCCGCACAGCATCTCATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCGCTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCt               <  2:67458/251‑1 (MQ=255)                                                                                                                                                                                                                      gATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCt               >  2:199820/1‑251 (MQ=255)                                                                                                                                                                                                                      gATTGCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCt               <  1:196919/251‑1 (MQ=255)                                                                                                                                                                                                                          gCCGCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTGCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAAt           <  1:172598/251‑1 (MQ=255)                                                                                                                                                                                                                             gCACAGCATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCtt        <  2:57089/251‑1 (MQ=255)                                                                                                                                                                                                                                   cATCACATAGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTa  <  1:124309/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   TGTGGTCATCATCACTCATTTATCAACCGTTTGATGCCGGCATATCGCGGGATGATCAATAAGATGAGTGCGGATTTGCTCATCGTGCCGTTTATCGATAAGTAACCGCTTACCAGTAATGCTCCGCTGTCATATGGCCCGGTCGGCGACGTAAATGTTTCGTCATCTGCCGGGTCTCTTTCAGCAACTGTTGTGTATCGCGCACCATCTGTGGATTGCCGCACAGCATCACATGGCTGGTTTCTTTATTCATCGGCAGGCCAATCGTGCTTTCCAGTTCCCCACTTTCAATTAATGCCGGTATCCGTCCGGTGAGCGACCCCGCTGCCGTTTCCCGACTGACCACCGTCTGAATGCGCAGTTTTCCTTCGTAGCGTTTTTCCAGTTCCTGCATCAGTGGCAAATAGCTTAAGTCGGCGGCATAACGTGCGGCGTGCACCAGGACCAGATTTTTGAAGCGATCTAAATCTTTACCTA  >  NC_000913/4113617‑4114093

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