Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,299,499	A→G	L27S (TTA→TCA)	napH ←	ferredoxin‑type protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,299,499	0	A	G	100.0%	193.6 / NA	55	L27S (TTA→TCA)	napH	ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/24);  total (31/24)

TAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAG  >  NC_000913/2299402‑2299595                                                                                                  |                                                                                                 tAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGac                                                                                                <  1:53331/100‑1 (MQ=255)  aaCAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGAACCGCTCAACAACATCCCCAGAACGAAGAAATGACAAAGGCGACGTGaca                                                                                               <  2:791327/100‑1 (MQ=255)   aCAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACCCCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGacac                                                                                              >  2:1212695/1‑100 (MQ=255)     aGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCa                                                                                            >  2:578681/1‑100 (MQ=255)        cTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGcc                                                                                         >  1:948017/1‑100 (MQ=255)        cTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGcc                                                                                         <  1:347849/100‑1 (MQ=255)             gctgTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACgg                                                                                    <  1:54917/100‑1 (MQ=255)                     tGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTgcg                                                                            <  1:574908/100‑1 (MQ=255)                       ccGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGcc                                                                          >  1:1107608/1‑100 (MQ=255)                       ccGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGcc                                                                          >  2:771590/1‑100 (MQ=255)                       ccGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGcc                                                                          <  2:448577/100‑1 (MQ=255)                         gTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGccac                                                                        >  1:402406/1‑100 (MQ=255)                            cAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAg                                                                     <  2:1211498/100‑1 (MQ=255)                             aaGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGc                                                                    >  2:976502/1‑100 (MQ=255)                             aaGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGc                                                                    >  2:976308/1‑100 (MQ=255)                                 tCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCttt                                                                <  1:529351/100‑1 (MQ=255)                                  ccACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCtttt                                                               >  2:561612/1‑100 (MQ=255)                                    acacacCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCtttttt                                                             <  2:639197/100‑1 (MQ=255)                                      acacCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTtct                                                           <  1:1036393/100‑1 (MQ=255)                                        acCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTcc                                                         <  1:237858/100‑1 (MQ=255)                                        acCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTcc                                                         <  2:132978/100‑1 (MQ=255)                                        acCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTcc                                                         <  2:371196/100‑1 (MQ=255)                                        acCAAACCATGAACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCACCCTTTTTTCTcc                                                         >  2:823985/1‑100 (MQ=255)                                         ccAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCa                                                        >  2:912337/1‑100 (MQ=255)                                          cAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAg                                                       >  1:1179934/1‑100 (MQ=255)                                          cAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAg                                                       <  2:64664/100‑1 (MQ=255)                                            aaCCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAgcg                                                     <  1:344644/100‑1 (MQ=255)                                                 tGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTcg                                                >  1:316606/1‑100 (MQ=255)                                                   gACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTcgcg                                              >  2:787555/1‑100 (MQ=255)                                                       gCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCg                                          >  2:260448/1‑100 (MQ=255)                                                          cAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCg                                       <  2:306942/100‑1 (MQ=255)                                                          cAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCg                                       >  1:1162405/1‑100 (MQ=255)                                                           aaaaaCATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGt                                      >  2:475261/1‑100 (MQ=255)                                                            aaaaCATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTc                                     >  1:67325/1‑100 (MQ=255)                                                             aaaCATCCCCAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCa                                    <  1:561616/100‑1 (MQ=255)                                                                    cccAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTAc                             >  1:1256212/1‑100 (MQ=255)                                                                    cccAGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTAc                             <  1:322566/100‑1 (MQ=255)                                                                       aGCACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACg                            >  2:899059/1‑98 (MQ=255)                                                                         cACGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGAttt                        >  2:334197/1‑100 (MQ=255)                                                                           cGAAGAACTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGAtt                         >  2:504948/1‑97 (MQ=255)                                                                            gaagaaCTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGcc                     <  2:727030/100‑1 (MQ=255)                                                                            gaagaaCTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTCTGCAAGCGCTTAGCGACCTGCGCGGCGTGTACGCTTTGcc                     >  2:490750/1‑100 (MQ=255)                                                                             aagaaCTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCa                    >  2:590788/1‑100 (MQ=255)                                                                             aagaaCTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCa                    <  1:426084/100‑1 (MQ=255)                                                                               gaaCTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCAtt                  >  1:694147/1‑100 (MQ=255)                                                                                aaCTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTTTTTTTCTTCTCGCGCCCGGCGTCACGTTTACGTTTTGCCATTg                 >  2:882828/1‑100 (MQ=255)                                                                                aaCTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTg                 <  1:767357/100‑1 (MQ=255)                                                                                 aCTGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGt                <  2:523883/100‑1 (MQ=255)                                                                                   tGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGCTTTGCCATTGTTc              >  2:960062/1‑100 (MQ=255)                                                                                   tGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTc              >  2:602537/1‑100 (MQ=255)                                                                                   tGACAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTc              >  2:651073/1‑100 (MQ=255)                                                                                      cAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTcccc           >  2:394798/1‑100 (MQ=255)                                                                                      cAAAGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTcccc           <  1:1112574/100‑1 (MQ=255)                                                                                       aaaGGCGACGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCt          <  1:61718/100‑1 (MQ=255)                                                                                               cGTGACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTccagccag  >  1:555159/1‑100 (MQ=255)                                                                                                  |                                                                                                 TAACAGGCTACTGCTGTAGTTGCCGTGCAAGATCCACACACCAAACCATGGACCGCTCAAAAACATCCCCAGCACGAAGAACTGACAAAGGCGACGTAACACCAGCCAACGGTGACTGCGCCACCAGCCTTTTTTCTCCAGCGCCTCGCGCCCGGCGTCACGTTTACGATTTGCCATTGTTCCCCTCCAGCCAG  >  NC_000913/2299402‑2299595

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