Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,058,514	A→G	Y144C (TAC→TGC)	torC →	cytochrome c menaquinol dehydrogenase TorC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,058,514	0	A	G	100.0%	153.3 / NA	43	Y144C (TAC→TGC)	torC	cytochrome c menaquinol dehydrogenase TorC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/20);  total (23/20)

TGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGT  >  NC_000913/1058377‑1058650                                                                                                                                          |                                                                                                                                         agATATCTCCCAGACCTTTATGGCTCACTCCTTTAATACACCTGAAAAATTCGAAGCCAAACGCGGGGAACTTGCCGAGCGTGAATGGGCGCGAAAGAAAGAAACCAACTGGCCACCCTGCCGCTCCTGCCATAACTGcg                                                                                                                                        <  2:58799/139‑1 (MQ=255)       cTACCAGCCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAGCAACTCGGCAACCTGCCGCTCCTGCCATAACtgcgatgcga                                                                                                                                  <  1:651909/140‑1 (MQ=255)           cAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGa                                                                                                                              >  1:37336/1‑140 (MQ=255)           agacagtttATTGCTCCCCCCATTAATACCCTTGAAAAATTCGAAGCCAACCGCGCGGAACTTGCCGAGGGTGAAGGGGCGCGAATAAAAAAAAACAACCCGGCACCCTGCCGCTCCTGCCATAACTGCGATGCGATGGa                                                                                                                              <  2:529667/134‑1 (MQ=255)              aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATct                                                                                                                           >  1:58799/1‑139 (MQ=255)              aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCa                                                                                                                           >  1:194911/1‑140 (MQ=255)              aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCa                                                                                                                           >  1:74590/1‑140 (MQ=255)                  ttATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCg                                                                                                                       >  2:452753/1‑140 (MQ=255)                  ttATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCg                                                                                                                       >  1:737265/1‑140 (MQ=255)                     tggCTCCCCCCATTGATACCCCTAAAAAATTCGAAGCCAACCGCGCGGAACTTGCCGAGCGTGAAGGGGCGCAAATGAAAGAAAACACCCCGGCACCCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAag                                                                                                                    <  2:37336/138‑1 (MQ=255)                       gCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAagca                                                                                                                  >  1:450498/1‑140 (MQ=255)                        cTCACTCCATTGATACACCTGAAAAATTCGAAGCAAAACGCGCGAAACTTGCCGAGCGTGAAGGGGCGCGAATGAAAAAAAACACCTCGGCACCCTGCCGCCCCTGCCATAACTGCGATGCGATGGATCATGCGAagcag                                                                                                                 <  2:728164/140‑1 (MQ=255)                            cTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAagcctgtct                                                                                                             >  1:529667/1‑134 (MQ=255)                             tCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATcc                                                                                                            >  2:13252/1‑140 (MQ=255)                                       cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAgcagca                                                                                                  >  2:770491/1‑140 (MQ=255)                                       cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAgcagca                                                                                                  >  2:321141/1‑140 (MQ=255)                                       cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAgcagca                                                                                                  <  2:194911/140‑1 (MQ=255)                                       cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAgcagca                                                                                                  >  1:758513/1‑140 (MQ=255)                                         ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACg                                                                                                >  1:56196/1‑140 (MQ=255)                                          ctaaaaaaTTCGAACCAAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATAAAAAAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGt                                                                                               <  2:445152/137‑1 (MQ=255)                                          cTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGt                                                                                               >  1:302302/1‑140 (MQ=255)                                           tGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTc                                                                                              >  1:508517/1‑140 (MQ=255)                                            gAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCa                                                                                             >  2:605510/1‑140 (MQ=255)                                            gAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCa                                                                                             >  2:347168/1‑140 (MQ=255)                                             aaaaaTTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATAAAAAAAAACAACCCGGCACCCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCACCATCCTGAAGCAGCACGTCAg                                                                                            <  2:574638/140‑1 (MQ=255)                                                    cGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAgg                                                                                     >  1:695840/1‑140 (MQ=255)                                                           aaaCGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGc                                                                              <  1:13252/140‑1 (MQ=255)                                                                     atcTTGCCGGGCGTAAATGGGCGCGAATGAAAAAAAACAACTCGGCACCCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATaat                                                                    <  2:764102/138‑1 (MQ=255)                                                                     aaCTGGCCGGGGGTAATTGGGCGCAAATGAAAGAAAACACCTCGGCAACCTGCCGCTCCTGCCATACCTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATaat                                                                    <  2:501000/140‑1 (MQ=255)                                                                               gCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCa                                                          >  1:792846/1‑140 (MQ=255)                                                                                    aaTGGGCGCAAATAAAAAAAAACAACTCGGCACCCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                     <  1:188235/140‑1 (MQ=255)                                                                                     aggggCGCGAATAAAAAAAAACACCCCGGCACCCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACt                                                    <  2:792846/138‑1 (MQ=255)                                                                                           gcgAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATa                                              >  2:55983/1‑140 (MQ=255)                                                                                              aaTGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAg                                           >  2:415388/1‑140 (MQ=255)                                                                                              aaTGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAg                                           >  2:606356/1‑140 (MQ=255)                                                                                               aTGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAgg                                          <  1:516920/140‑1 (MQ=255)                                                                                                      aaaaCAACTCGGCAACCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGcc                                   <  2:74590/140‑1 (MQ=255)                                                                                                                ggCACCCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                         <  2:496319/140‑1 (MQ=255)                                                                                                                       tGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGatat                  <  1:606356/140‑1 (MQ=255)                                                                                                                       tGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGatat                  >  1:205838/1‑140 (MQ=255)                                                                                                                         ccGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGa                <  2:56196/140‑1 (MQ=255)                                                                                                                                tGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCTCCAGTTACCGGATATGAGTAGCgg         <  2:522882/140‑1 (MQ=255)                                                                                                                                  ccATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCt       <  2:508517/140‑1 (MQ=255)                                                                                                                                     tAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTcc    <  2:597224/140‑1 (MQ=255)                                                                                                                                       aCTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGt  <  1:415388/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         TGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGT  >  NC_000913/1058377‑1058650

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