Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,292,906	A→G	A465A (GCA→GCG)	nrfE →	putative cytochrome c‑type biogenesis protein NrfE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,292,906	0	A	G	100.0%	187.4 / NA	51	A465A (GCA→GCG)	nrfE	putative cytochrome c‑type biogenesis protein NrfE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/19);  total (32/19)

CAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAA  >  NC_000913/4292771‑4293042                                                                                                                                        |                                                                                                                                         cAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGt                                                                                                                                      >  1:781580/1‑140 (MQ=255)    cGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTg                                                                                                                                   <  2:850885/140‑1 (MQ=255)      tCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGAc                                                                                                                                 >  1:43984/1‑140 (MQ=255)         gCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTg                                                                                                                                   >  2:486281/1‑135 (MQ=255)         gCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTg                                                                                                                                   <  1:486281/135‑1 (MQ=255)               aGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGtt                                                                                                                                     >  2:285192/1‑127 (MQ=255)               aGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGtt                                                                                                                                     <  1:285192/127‑1 (MQ=255)                 gagaTCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACg                                                                                                                      <  1:519209/140‑1 (MQ=255)                 gagaTCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACg                                                                                                                      <  1:454153/140‑1 (MQ=255)                         ccTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTg                                                                                                              >  1:719551/1‑140 (MQ=255)                         ccTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTg                                                                                                              >  2:670354/1‑140 (MQ=255)                         ccTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTg                                                                                                              >  1:353908/1‑140 (MQ=255)                                  aCAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGc                                                                                                     >  1:736930/1‑140 (MQ=255)                                       cGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAg                                                                                                <  1:395191/140‑1 (MQ=255)                                                           gCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACg                                                                            <  1:673252/140‑1 (MQ=255)                                                            cAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGc                                                                           >  1:415777/1‑140 (MQ=255)                                                                         ccGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGcc                                                                          <  2:191413/128‑1 (MQ=255)                                                                         ccGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGcc                                                                          >  1:191413/1‑128 (MQ=255)                                                                         ccGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAatga                                                              >  1:278196/1‑140 (MQ=255)                                                                         ccGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAatga                                                              >  2:134327/1‑140 (MQ=255)                                                                          cGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTc                                                                       >  2:75550/1‑130 (MQ=255)                                                                          cGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTc                                                                       <  1:75550/130‑1 (MQ=255)                                                                          cGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAatgat                                                             >  2:744509/1‑140 (MQ=255)                                                                                   tGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTc                                                    >  1:52161/1‑140 (MQ=255)                                                                                      cTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGGACTGACGCCAGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAAt                                                 >  2:825544/1‑140 (MQ=255)                                                                                          aTCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCAGTCAACAAATGATGGAACCGTCAATTCGc                                             <  2:415777/140‑1 (MQ=255)                                                                                           tCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCt                                            >  2:167925/1‑140 (MQ=255)                                                                                               cAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGaa                                        >  2:844765/1‑140 (MQ=255)                                                                                               cAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGaa                                        >  2:2689/1‑140 (MQ=255)                                                                                                 aGCCAAAGGCAATTGCACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACg                                      <  1:670354/140‑1 (MQ=255)                                                                                                    cAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCa                                   >  1:849583/1‑140 (MQ=255)                                                                                                     aaaGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCAt                                  >  2:448905/1‑140 (MQ=255)                                                                                                     aaaGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCAt                                  >  1:96797/1‑140 (MQ=255)                                                                                                     aaaGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCAt                                  >  1:889982/1‑140 (MQ=255)                                                                                                       aGGCAATTACACGAGGGAAAAAGCGATAGAGGCGCTGTATGACCATCAGCAACGCATTGGTGGACTGTCGCCGGAGCGGCGTTTTTATGATGTACGCCGTCAACAAATGATGGAACCGTCAATTCTCTGGTACGGCATcc                                >  2:496840/1‑140 (MQ=255)                                                                                                       aGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATcc                                >  2:651321/1‑140 (MQ=255)                                                                                                       aGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATcc                                >  2:873872/1‑140 (MQ=255)                                                                                                       aGGCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATcc                                >  1:542411/1‑140 (MQ=255)                                                                                                        ggCAATTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCa                               <  1:167925/140‑1 (MQ=255)                                                                                                           aaTTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGa                            <  2:52161/140‑1 (MQ=255)                                                                                                            aTTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGAc                           <  2:353908/140‑1 (MQ=255)                                                                                                            aTTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGAc                           <  2:719551/140‑1 (MQ=255)                                                                                                            aTTACACCAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGAc                           >  1:119568/1‑140 (MQ=255)                                                                                                                  ccAGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTAt                     >  2:66190/1‑140 (MQ=255)                                                                                                                    aGCGAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGc                   <  1:385872/140‑1 (MQ=255)                                                                                                                       gAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGt                >  2:264448/1‑140 (MQ=255)                                                                                                                       gAAAAAGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGt                >  2:133091/1‑140 (MQ=255)                                                                                                                        aaaaaGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTc               >  2:16486/1‑140 (MQ=255)                                                                                                                         aaaaGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCa              <  2:710911/140‑1 (MQ=255)                                                                                                                           aaGCGATAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATg            <  1:264448/140‑1 (MQ=255)                                                                                                                                aTAGTGGCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGga       >  2:235496/1‑140 (MQ=255)                                                                                                                                     ggCGCTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGaaaa  <  2:849583/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         CAGCGTCAGCCGTCAGGAGATCAGCCTGAATTTACAGCCGGGTCAGCAGGTGACGCTGGCAGGATACACCTTCCGTTTTGAGTGCCTCGATCTACAAGCCAAAGGCAATTACACCAGCGAAAAAGCGATAGTGGCACTGTTTGACCATCAGCAACGCATTGGTGAACTGACGCCGGAGCGGCGTTTTTATGAAGCACGCCGTCAACAAATGATGGAACCGTCAATTCGCTGGAACGGCATCCATGACTGGTATGCGGTCATGGGGGAGAAAA  >  NC_000913/4292771‑4293042

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