Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	449,344	A→G	F433S (TTC→TCC)	cyoB ←	cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	449,344	0	A	G	100.0%	184.7 / NA	51	F433S (TTC→TCC)	cyoB	cytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/23);  total (28/23)

GACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAGGTGACGATAAAACCG  >  NC_000913/449198‑449487                                                                                                                                                   |                                                                                                                                                gacgGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGagc                                                                                                                                              >  1:111480/1‑150 (MQ=255)        cATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAg                                                                                                                                      <  2:805023/151‑1 (MQ=255)           gCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGAc                                                                                                                                    <  1:168364/150‑1 (MQ=255)             ccATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCa                                                                                                                                  >  1:244360/1‑150 (MQ=255)                  aaGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCAcgccgc                                                                                                                            >  1:260858/1‑151 (MQ=255)                         gcgcATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCgatgat                                                                                                                     >  2:974714/1‑151 (MQ=255)                         gcgcATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCgatga                                                                                                                      <  2:183503/150‑1 (MQ=255)                           gcATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCa                                                                                                                   <  1:17054/151‑1 (MQ=255)                             aTACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACg                                                                                                                 >  2:145776/1‑151 (MQ=255)                               aCAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGtt                                                                                                               >  1:909749/1‑151 (MQ=255)                                cAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTa                                                                                                              <  1:505842/151‑1 (MQ=255)                                      cATAAAGGCACCGAAGAAGCCGATGATCCAGAACCAGACCGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGATCGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCTCGTTATGGAAg                                                                                                        <  2:441357/151‑1 (MQ=255)                                             gCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCaa                                                                                                 <  2:405575/151‑1 (MQ=255)                                               aaCGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATc                                                                                               >  1:678434/1‑151 (MQ=255)                                                 cGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAg                                                                                             <  2:1062425/151‑1 (MQ=255)                                                   aagaagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAgg                                                                                            >  2:98942/1‑150 (MQ=255)                                                    agaagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGa                                                                                           >  1:721900/1‑150 (MQ=255)                                                     gaagCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAAc                                                                                         >  2:642181/1‑151 (MQ=255)                                                          cGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCt                                                                                    >  2:733302/1‑151 (MQ=255)                                                           gatgatCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTg                                                                                   <  2:410939/151‑1 (MQ=255)                                                                   agaaccagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACAGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGc                                                                             >  2:181912/1‑149 (MQ=255)                                                                     aaccagaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGa                                                                          >  1:734634/1‑150 (MQ=255)                                                                         agaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGaa                                                                     <  1:718369/151‑1 (MQ=255)                                                                          gaacGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGa                                                                      >  2:246634/1‑149 (MQ=255)                                                                               cgcgTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCcgc                                                               <  2:721862/151‑1 (MQ=255)                                                                                      aCCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTGTAGGCCACCAGTAGGGCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGAACACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCgg                                                         >  2:419603/1‑150 (MQ=255)                                                                                       ccccAGGTTTCGTTCAGTTTGAAACCGAACTCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCgg                                                         <  2:561990/149‑1 (MQ=255)                                                                                       ccccAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTa                                                       <  1:705333/151‑1 (MQ=255)                                                                                            ggTTTCGTACAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGcc                                                  <  2:313533/151‑1 (MQ=255)                                                                                                tCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagc                                               >  1:668812/1‑150 (MQ=255)                                                                                                 cGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCACTCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagc                                               >  2:866794/1‑149 (MQ=255)                                                                                                 cGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCcagc                                               >  2:865344/1‑149 (MQ=255)                                                                                                           gAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTc                                   >  1:73277/1‑151 (MQ=255)                                                                                                              aCCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATc                                >  2:566660/1‑151 (MQ=255)                                                                                                               ccGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATcc                               >  2:1011827/1‑151 (MQ=255)                                                                                                                   aCGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCg                             <  1:881297/149‑1 (MQ=255)                                                                                                                   aCGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGcc                           <  1:880176/151‑1 (MQ=255)                                                                                                                   aCGCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGc                            >  1:113667/1‑150 (MQ=255)                                                                                                                     gCTTTAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGccc                          >  1:110961/1‑150 (MQ=255)                                                                                                                       tttAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCCGCACGCCAGTCATCCCGCCCa                         <  2:233205/149‑1 (MQ=255)                                                                                                                        ttAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCg                      >  2:823855/1‑151 (MQ=255)                                                                                                                         tAGGCCACCAGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCg                      <  1:279427/150‑1 (MQ=255)                                                                                                                             ccaccaGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAg                 <  1:79652/151‑1 (MQ=255)                                                                                                                                 caGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAGGTGa             <  1:974765/151‑1 (MQ=255)                                                                                                                                 caGTAGGTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAGGTGa             <  1:943147/151‑1 (MQ=255)                                                                                                                                       gTCATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAGGTGACGATa        >  2:425075/1‑150 (MQ=255)                                                                                                                                         cATCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAGGTGACGATaaaa     <  1:4077/151‑1 (MQ=255)                                                                                                                                          aTCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAGGTGACGATAAAAc    >  1:668906/1‑151 (MQ=255)                                                                                                                                          aTCCCTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAGGTGACGATAAAAc    >  2:794377/1‑151 (MQ=255)                                                                                                                                            cccTGCGGAGCAGCCGAAGACCACCCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGGGCCCGGTTCGGCCAGCAGCACGACAGTCACCCCGCCCACCGAGAAGGTGACGATAAAAcc   <  2:778011/150‑1 (MQ=255)                                                                                                                                             ccTGCGGAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAGGTGACGATAAAACCg  >  2:912021/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                GACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAACGCGCGTTTACCCCAGGTTTCGTTCAGTTTGAAACCGAACGCTTTAGGCCACCAGTAGGTCATCCCTGCGAAGCAGCCGAAGACCACGCCGCCGATGATCACGTTATGGAAGTGCGCAATCAGGAACAGGCTGTTATGCAGAACGAAGTCCGCGCCCGGTACGGCCAGCAGCACGCCAGTCATCCCGCCCACCGAGAAGGTGACGATAAAACCG  >  NC_000913/449198‑449487

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