Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,303,421	C→A	W26L (TGG→TTG)	napF ←	ferredoxin‑type protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,303,421	0	C	A	100.0%	118.3 / NA	35	W26L (TGG→TTG)	napF	ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (19/16);  total (19/16)

AGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAATTGAGAGCGACCAAATAAACC  >  NC_000913/2303287‑2303550                                                                                                                                       |                                                                                                                                  aGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAggg                                                                                                                               >  1:152188/1‑139 (MQ=255) aGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAggg                                                                                                                               <  2:389680/139‑1 (MQ=255) aGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAggg                                                                                                                               >  1:389680/1‑139 (MQ=255) aGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAggg                                                                                                                               <  2:152188/139‑1 (MQ=255) aGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGc                                                                                                                              >  1:323464/1‑140 (MQ=255)       aCTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGcggacgg                                                                                                                        >  1:205657/1‑140 (MQ=255)        cTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGcggacgga                                                                                                                       <  1:259490/140‑1 (MQ=255)                 tttGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTAc                                                                                                              >  2:301267/1‑140 (MQ=255)                         tAACGCTTGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTg                                                                                                      >  1:219507/1‑140 (MQ=255)                                cGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACTGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGc                                                                                               <  2:54184/140‑1 (MQ=255)                                   aTAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGAc                                                                                            <  1:64406/140‑1 (MQ=255)                                       ccgcccgcGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGt                                                                                        >  2:327627/1‑140 (MQ=255)                                            cgcgccgcgTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTa                                                                                   >  2:315895/1‑140 (MQ=255)                                                        cAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACggg                                                                       >  1:98233/1‑140 (MQ=255)                                                         aGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGa                                                                      <  2:205657/140‑1 (MQ=255)                                                             tgttgtTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCa                                                                  >  1:274415/1‑140 (MQ=255)                                                             tgttgtTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCa                                                                  >  2:276246/1‑140 (MQ=255)                                                                tgttTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCa                                                               >  1:447072/1‑140 (MQ=255)                                                                   tttCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATc                                                            <  1:128745/140‑1 (MQ=255)                                                                     tCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCtt                                                          <  2:219507/140‑1 (MQ=255)                                                                       acaAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTca                                                        <  2:362827/140‑1 (MQ=255)                                                                         aaGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTcaca                                                      >  1:449614/1‑140 (MQ=255)                                                                              ttAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGAc                                                 <  2:154213/140‑1 (MQ=255)                                                                                  tACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTc                                             <  1:301267/140‑1 (MQ=255)                                                                                      ggCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCcatc                                         >  1:67642/1‑140 (MQ=255)                                                                                       gCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCcatca                                        >  2:259395/1‑140 (MQ=255)                                                                                       gCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCcatca                                        >  1:175014/1‑140 (MQ=255)                                                                                           cGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTaa                                    >  1:61703/1‑140 (MQ=255)                                                                                           cGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTaa                                    >  2:3944/1‑140 (MQ=255)                                                                                                  ggTGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGctc                             <  1:276246/140‑1 (MQ=255)                                                                                                    tGCAATGGGTCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTg                           <  2:61703/140‑1 (MQ=255)                                                                                                             tCAGAAAATGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAATTgag                  >  2:385720/1‑140 (MQ=255)                                                                                                                 aaaaTGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAATTGAGAGCg              <  1:259395/140‑1 (MQ=255)                                                                                                                   aaTGAGATTCATCACCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAATTGAGAGCGAc            <  2:86265/140‑1 (MQ=255)                                                                                                                             atcaCCGGACAAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAATTGAGAGCGACCAAATAAAcc  <  2:449614/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                  AGCTGCACTCGTTATTTTTGAAATTAACGCTCGGATAGCCGCCCGCGCCGCGTTGCAGAATGTTGTTTTCACAAGCATTAATACAGGCGTCGCAGCGGGTGCAATGGGTCAGAAAATGAGATTCATCACCGGACCAGGGCGGACGGATACCGTTACTGGCTTTGCGCCAGCGACCAGTGAGTATGCCCCGACGGGATGCATCAATCTTCACATTGACCTTCCATCATTAACGCGCTCTGAAAATTGAGAGCGACCAAATAAACC  >  NC_000913/2303287‑2303550

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