Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	450,764	C→T	A283T (GCT→ACT)	cyoA ←	cytochrome o ubiquinol oxidase subunit II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	450,764	0	C	T	100.0%	232.8 / NA	71	A283T (GCT→ACT)	cyoA	cytochrome o ubiquinol oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/37);  total (34/37)

TTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGCCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGT  >  NC_000913/450644‑450883                                                                                                                         |                                                                                                                        ttATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTAc                                                                                                       <  2:891058/139‑1 (MQ=255) ttATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTAc                                                                                                       <  1:1785217/139‑1 (MQ=255)       ttcCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGc                                                                                                 >  1:360000/1‑139 (MQ=255)         ccTCAACCCCTTTAATGGGTGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                               >  1:1930505/1‑139 (MQ=255)         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                               >  2:937748/1‑139 (MQ=255)         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                               >  2:917313/1‑139 (MQ=255)         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                               >  2:1679697/1‑139 (MQ=255)            cAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAAc                                                                                             <  2:1467404/138‑1 (MQ=255)            cAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAAc                                                                                             >  1:1467404/1‑138 (MQ=255)               ccccTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACt                                                                                                                 >  1:980519/1‑115 (MQ=255)               ccccTTTAATGGGCGGATTCCGCGTGGCGCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCACTTACCGTGAGTCATAAACt                                                                                                                 <  2:980519/115‑1 (MQ=255)                  cTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTg                                                                                      >  2:988054/1‑139 (MQ=255)                  cTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTg                                                                                      >  2:1967685/1‑139 (MQ=255)                  cTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTg                                                                                      >  1:178269/1‑139 (MQ=255)                   tttAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAAc                                                                                             >  2:469048/1‑131 (MQ=255)                   tttAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAAc                                                                                             <  1:469048/131‑1 (MQ=255)                   tttAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTgg                                                                                     >  2:686635/1‑139 (MQ=255)                     tAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGtt                                                                                   <  2:1580828/139‑1 (MQ=255)                      aaTGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGttt                                                                                  <  1:988054/139‑1 (MQ=255)                      aaTGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGttt                                                                                  <  2:1018513/139‑1 (MQ=255)                          ggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACg                                                                              <  1:686635/139‑1 (MQ=255)                            cGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAAc                                                                                             <  2:1409258/122‑1 (MQ=255)                            cGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAAc                                                                                             >  1:1409258/1‑122 (MQ=255)                             ggATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTg                                                                           <  1:1679697/139‑1 (MQ=255)                                 tCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGgaga                                                                       <  2:178269/139‑1 (MQ=255)                                         gCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                  <  1:663948/88‑1 (MQ=255)                                         gCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                  >  2:663948/1‑88 (MQ=255)                                           tCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAc                                                             <  1:1819334/139‑1 (MQ=255)                                               gTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTgg                                                         <  2:360000/139‑1 (MQ=255)                                               gTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTgg                                                         <  1:937748/139‑1 (MQ=255)                                                  cATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                  <  1:761546/79‑1 (MQ=255)                                                  cATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                  >  2:761546/1‑79 (MQ=255)                                                  cATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTg                                                      <  2:1930505/139‑1 (MQ=255)                                                   aTGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGt                                                     >  2:1662114/1‑139 (MQ=255)                                                                  tCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCgg                                      >  1:1449882/1‑139 (MQ=255)                                                                   cGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACtt                                                                                                                <  2:199406/64‑1 (MQ=255)                                                                   cGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACtt                                                                                                                >  1:199406/1‑64 (MQ=255)                                                                   cGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTa                                                                                                        >  2:1238142/1‑72 (MQ=255)                                                                   cGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTa                                                                                                        <  1:1238142/72‑1 (MQ=255)                                                                     tgtgCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCa                                   >  2:878040/1‑139 (MQ=255)                                                                       tgCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACg                                                                              >  1:248471/1‑94 (MQ=255)                                                                       tgCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACg                                                                              <  2:248471/94‑1 (MQ=255)                                                                              tgCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGa                          >  1:968275/1‑139 (MQ=255)                                                                                   aCCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCa                     <  2:968275/139‑1 (MQ=255)                                                                                    ccTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAg                    <  1:917313/139‑1 (MQ=255)                                                                                     cTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGc                   <  2:1449882/139‑1 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCg                                                                                                   <  2:1758047/56‑1 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCg                                                                                                   >  1:1758047/1‑56 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTg                                                                           >  1:1745445/1‑80 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTg                                                                           <  2:1745445/80‑1 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTa                                                    <  2:1163561/103‑1 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTa                                                    >  1:1163561/1‑103 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTGGCTAGGCGCGGCCAGTTTTTCGAACGCAGc                   <  2:523840/136‑1 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTGGCTAGGCGCGGCCAGTTTTTCGAACGCAGc                   >  1:523840/1‑136 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAAc                        <  2:1925237/131‑1 (MQ=255)                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGCTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAAc                        >  1:1925237/1‑131 (MQ=255)                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAg                                                                                          <  1:767885/58‑1 (MQ=255)                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAg                                                                                          >  2:767885/1‑58 (MQ=255)                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTc                                                                                        >  2:1240501/1‑60 (MQ=255)                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTc                                                                                        <  1:1240501/60‑1 (MQ=255)                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAc                                                             >  1:1520774/1‑87 (MQ=255)                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAc                                                             <  2:1520774/87‑1 (MQ=255)                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCt                                              <  2:1578664/102‑1 (MQ=255)                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCt                                              >  1:1578664/1‑102 (MQ=255)                                                                                                      tCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGgtgt  <  2:605702/139‑1 (MQ=255)                                                                                                               ttACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGt                                                                                         <  1:1983982/43‑1 (MQ=255)                                                                                                               ttACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGt                                                                                         >  2:1983982/1‑43 (MQ=255)                                                                                                                 aCCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAata                                                                   >  1:1652791/1‑63 (MQ=255)                                                                                                                 aCCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAata                                                                   <  2:1652791/63‑1 (MQ=255)                                                                                                                     tGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACg                                                                              <  1:1648559/48‑1 (MQ=255)                                                                                                                     tGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACg                                                                              >  2:1648559/1‑48 (MQ=255)                                                                                                                         |                                                                                                                        TTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGCCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGT  >  NC_000913/450644‑450883

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