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

AATTTTCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGCCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTC  >  NC_000913/450628‑450897                                                                                                                                         |                                                                                                                                      aaTTTTCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCa                                                                                                                                     <  1:1274222/139‑1 (MQ=255)      tCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGTTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                <  2:1575662/139‑1 (MQ=255)        cGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACtt                                                                                                                              <  1:827942/139‑1 (MQ=255)                tttATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTa                                                                                                                      <  2:1491902/139‑1 (MQ=255)                       ttcCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGc                                                                                                               >  1:1154846/1‑139 (MQ=255)                         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                             >  1:337577/1‑139 (MQ=255)                         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                             >  1:1433774/1‑139 (MQ=255)                         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                             >  2:1356069/1‑139 (MQ=255)                         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                             >  1:831149/1‑139 (MQ=255)                         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                             >  1:973247/1‑139 (MQ=255)                         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                             >  1:838086/1‑139 (MQ=255)                         ccTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaa                                                                                                             >  2:790857/1‑139 (MQ=255)                          cTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCaaa                                                                                                            >  2:1318128/1‑139 (MQ=255)                               ccccTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGt                                                                                                       >  2:456914/1‑139 (MQ=255)                               ccccTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGt                                                                                                       <  2:186275/139‑1 (MQ=255)                                    ttAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCTTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGt                                                                                                  >  1:286266/1‑139 (MQ=255)                                     tAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGtt                                                                                                 >  2:479750/1‑139 (MQ=255)                                      aaTGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGttt                                                                                                <  2:808740/139‑1 (MQ=255)                                         gggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCAc                                                                                             >  2:402200/1‑139 (MQ=255)                                          ggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCa                                                                                                                                     >  1:317370/1‑98 (MQ=255)                                          ggCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACg                                                                                            <  2:431255/139‑1 (MQ=255)                                                       tGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATAtt                                                                               <  1:1356069/139‑1 (MQ=255)                                                         gCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTAc                                                                                                                     <  1:374615/99‑1 (MQ=255)                                                         gCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTAc                                                                                                                     >  2:374615/1‑99 (MQ=255)                                                         gCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTcc                                                                             >  1:210755/1‑139 (MQ=255)                                                          cTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCa                                                                            <  2:1154846/139‑1 (MQ=255)                                                                    tGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTa                                                                  <  2:337577/139‑1 (MQ=255)                                                                    tGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTa                                                                  >  2:713342/1‑139 (MQ=255)                                                                         tCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGc                                                             <  1:398803/139‑1 (MQ=255)                                                                                cTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGgcgc                                                      <  1:713342/139‑1 (MQ=255)                                                                                       tgCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGt                                               >  1:602561/1‑139 (MQ=255)                                                                                         cgcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                <  2:335889/56‑1 (MQ=255)                                                                                         cgcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                >  1:335889/1‑56 (MQ=255)                                                                                         cgcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCg                                                                                                                 >  1:215058/1‑71 (MQ=255)                                                                                         cgcTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCg                                                                                                                 <  2:215058/71‑1 (MQ=255)                                                                                             gtgCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCa                                                                                              >  1:775502/1‑86 (MQ=255)                                                                                             gtgCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCa                                                                                              <  2:775502/86‑1 (MQ=255)                                                                                              tgCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGa                                        <  1:456914/139‑1 (MQ=255)                                                                                               gCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGaa                                       <  2:1433774/139‑1 (MQ=255)                                                                                                cTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAAc                                      <  2:831149/139‑1 (MQ=255)                                                                                                 tCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACg                                     <  2:1579849/139‑1 (MQ=255)                                                                                                    ccTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAg                                  <  2:838086/139‑1 (MQ=255)                                                                                                    ccTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAg                                  <  1:479750/139‑1 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATa                                                                                                                                   >  1:32449/1‑38 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATa                                                                                                                                   <  1:542654/38‑1 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATa                                                                                                                                   <  2:32449/38‑1 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                <  2:481992/41‑1 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAAc                                                                                                                                >  1:481992/1‑41 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGt                                                                                                       <  1:1208648/66‑1 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGt                                                                                                       >  2:1208648/1‑66 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAc                                                                           <  2:351317/94‑1 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCAc                                                                           >  1:351317/1‑94 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTgg                                                                       <  2:1673482/98‑1 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTgg                                                                       >  1:1673482/1‑98 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTa                                                                  <  1:806056/103‑1 (MQ=255)                                                                                                        ctggctggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTa                                                                  >  2:806056/1‑103 (MQ=255)                                                                                                        ctggatggGTCATGTCCATGCTCTTACCGTGAGTCATa                                                                                                                                   >  2:542654/1‑38 (MQ=255)                                                                                                             tggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTg                                                                                                                             >  2:1611787/1‑39 (MQ=255)                                                                                                             tggGTCATGTCCATGCTCTTACCGTGAGTCATAAACTTg                                                                                                                             <  1:1611787/39‑1 (MQ=255)                                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCa                                                                            >  1:459892/1‑86 (MQ=255)                                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCa                                                                            <  2:459892/86‑1 (MQ=255)                                                                                                               ggTCATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGAc                       >  1:1552319/1‑139 (MQ=255)                                                                                                                  cATGTCCATGCTCTTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATg                    <  2:210755/139‑1 (MQ=255)                                                                                                                          tgatctTACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGGTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCgg            >  2:1004111/4‑139 (MQ=255)                                                                                                                               ttACCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACt       >  2:1344926/1‑139 (MQ=255)                                                                                                                                 aCCGTGAGTCATAAACTTGTTAATTACATCGGCaaa                                                                                                            <  1:853142/36‑1 (MQ=255)                                                                                                                                 aCCGTGAGTCATAAACTTGTTAATTACATCGGCaaa                                                                                                            >  2:853142/1‑36 (MQ=255)                                                                                                                                 aCCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTgg                                                                       <  2:366615/73‑1 (MQ=255)                                                                                                                                 aCCGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTgg                                                                       >  1:366615/1‑73 (MQ=255)                                                                                                                                  ccGTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCt    <  1:402200/139‑1 (MQ=255)                                                                                                                                    gTGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTc  <  1:1318128/139‑1 (MQ=255)                                                                                                                                     tGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGcgcg                                                     <  2:287637/87‑1 (MQ=255)                                                                                                                                     tGAGTCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGcgcg                                                     >  1:287637/1‑87 (MQ=255)                                                                                                                                         |                                                                                                                                      AATTTTCCGAACATCTTTATTCTTCCTCAACCCCTTTAATGGGCGGATTCCGCGTGGCTCATGTCCATGCCTTCCATACCTTCGTGTGCGCTGTGCTCACCTTCTGGCTGGGTCATGTCCATGCTCTTACCGTGAGCCATAAACTTGTTAATTACATCGGCAAACAAGTCTGGTTTCACGTTGGAGAAATATTCCACCTGGTTGTATTCGCTAGGCGCGGCCAGTTTTTCGAACGCAGCCATGTCAGACATGGTGTTCGGCGACTGCTTC  >  NC_000913/450628‑450897

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: ‑