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%	179.1 / NA	56	A283T (GCT→ACT)	cyoA	cytochrome o ubiquinol oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/26);  total (30/26)

TTTTCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGCCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTCG  >  NC_000913/450630‑450898                                                                                                                                       |                                                                                                                                       ttttCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATa                                                                                                                                    <  2:5461/139‑1 (MQ=255)    tccagaCATATTTTTTCTTCCTCTACCCCTTTGATGGGTTGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTTGTTTGCGCTGTGCTTACTTTCTGGCTGGGTCAAGTCCATGCTCTAACCGTGAGTCATAGAc                                                                                                                                 <  2:378886/134‑1 (MQ=255)        aaCATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACGTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGt                                                                                                                             >  2:1391035/1‑139 (MQ=255)        aaCATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGt                                                                                                                             >  1:695656/1‑139 (MQ=255)                       ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                 <  2:928336/120‑1 (MQ=255)                       ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                 >  1:928336/1‑120 (MQ=255)                       ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                 <  2:1626592/120‑1 (MQ=255)                       ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                 >  1:1626592/1‑120 (MQ=255)                       ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                              >  2:1133532/1‑139 (MQ=255)                       ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                              >  1:1262156/1‑139 (MQ=255)                       ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                              >  1:597439/1‑139 (MQ=255)                       ccTCAACCCCTTTAATGGGCGGATTCCGAGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                              >  1:1026018/1‑139 (MQ=255)                        cTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGACTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTAGATCGGCaaa                                                                                                             >  2:978564/1‑139 (MQ=255)                           aaCCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACaa                                                                                                          >  1:549610/1‑139 (MQ=255)                              cccTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTc                                                                                                       <  2:1330478/139‑1 (MQ=255)                                     aTGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATa                                                                                                                                    <  1:745443/103‑1 (MQ=255)                                     aTGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATa                                                                                                                                    >  2:745443/1‑103 (MQ=255)                                       gggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCAc                                                                                              >  2:90421/1‑139 (MQ=255)                                        ggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACg                                                                                             <  2:175767/139‑1 (MQ=255)                                        ggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACg                                                                                             <  1:1114842/139‑1 (MQ=255)                                          cGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGtt                                                                                           >  2:588687/1‑139 (MQ=255)                                           ggATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTg                                                                                          <  2:743448/139‑1 (MQ=255)                                                   cgTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAata                                                                                  >  2:1332808/1‑139 (MQ=255)                                                   cgTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAata                                                                                  <  1:1311972/139‑1 (MQ=255)                                                         tCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAc                                                                            <  1:588687/139‑1 (MQ=255)                                                            tGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACAt                                                                                                                    >  2:894057/1‑96 (MQ=255)                                                            tGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACAt                                                                                                                    <  1:894057/96‑1 (MQ=255)                                                            tGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTg                                                                         <  1:1391035/139‑1 (MQ=255)                                                             gTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTgg                                                                        <  2:597439/139‑1 (MQ=255)                                                               ccATGCCTTCCATACCGTCGTGTGCGCTGTGCTCAC‑‑‑CTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTa                                                                   <  1:204696/139‑1 (MQ=255)                                                                  tGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTa                                                                   <  1:978564/139‑1 (MQ=255)                                                                       tCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGc                                                              >  1:84206/1‑139 (MQ=255)                                                                                tCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCgg                                                     <  1:1332808/139‑1 (MQ=255)                                                                                 cGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGc                                                    >  1:1469242/1‑139 (MQ=255)                                                                                        gcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCgg                                                     >  2:470040/1‑131 (MQ=255)                                                                                        gcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCgg                                                     <  1:470040/131‑1 (MQ=255)                                                                                        gcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGtttt                                             >  2:174019/1‑139 (MQ=255)                                                                                          tgtgCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTc                                           >  1:1389740/1‑139 (MQ=255)                                                                                                    ttCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGcc                                 <  1:1133532/139‑1 (MQ=255)                                                                                                    ttCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGcc                                 <  1:471925/139‑1 (MQ=255)                                                                                                    ttCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGcc                                 <  2:412130/139‑1 (MQ=255)                                                                                                      ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCAt                               <  1:90421/139‑1 (MQ=255)                                                                                                      ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCAt                               >  1:1415779/1‑139 (MQ=255)                                                                                                      ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCAt                               >  2:1139396/1‑139 (MQ=255)                                                                                                        ggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACAt                                                                                                                    >  1:1582823/1‑52 (MQ=255)                                                                                                        ggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACAt                                                                                                                    <  2:1582823/52‑1 (MQ=255)                                                                                                           tggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAg                          <  2:695656/139‑1 (MQ=255)                                                                                                             ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGAc                        >  2:28116/1‑139 (MQ=255)                                                                                                                   gTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTa                                                                                                                       >  1:1402406/1‑38 (MQ=255)                                                                                                                   gTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTa                                                                                                                       <  2:1402406/38‑1 (MQ=255)                                                                                                                   gTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGgtg                  >  1:834828/1‑139 (MQ=255)                                                                                                                        tGCTCTTACCGTGAGTCAAAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGATCgg             >  1:1091758/1‑139 (MQ=255)                                                                                                                               aCCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCAt                               <  2:326615/114‑1 (MQ=255)                                                                                                                               aCCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCAt                               >  1:326615/1‑114 (MQ=255)                                                                                                                                  gTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTc   >  2:38305/1‑139 (MQ=255)                                                                                                                                   tGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTCg  >  2:916727/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       TTTTCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGCCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTCG  >  NC_000913/450630‑450898

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