Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,389	A→G	K163R (AAG→AGG)	cybB →	cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,389	0	A	G	100.0%	234.8 / NA	66	K163R (AAG→AGG)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/38);  total (28/38)

GCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATA  >  NC_000913/1491258‑1491517                                                                                                                                    |                                                                                                                                 gCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGacaac                                                                                                                           >  2:160860/1‑139 (MQ=255) gCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGacaac                                                                                                                           >  1:191078/1‑139 (MQ=255)   cTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacac                                                                                                                         <  2:110175/139‑1 (MQ=255)                 gCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacaca                                                                                                                        <  2:353370/126‑1 (MQ=255)                 gCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacaca                                                                                                                        >  1:353370/1‑126 (MQ=255)                    aaTTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                               <  2:553123/116‑1 (MQ=255)                    aaTTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                               >  1:553123/1‑116 (MQ=255)                            aCGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTaaa                                                                                                <  1:71955/139‑1 (MQ=255)                                gTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCa                                                                                                              >  1:522976/1‑121 (MQ=255)                                gTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCa                                                                                                              <  2:522976/121‑1 (MQ=255)                                       gCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAg                                                                                     <  1:1311444/139‑1 (MQ=255)                                         ttAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCt                                                                                                                   <  2:1220192/107‑1 (MQ=255)                                         ttAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCt                                                                                                                   >  1:1220192/1‑107 (MQ=255)                                          tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                               <  2:1215642/94‑1 (MQ=255)                                          tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                               >  1:1215642/1‑94 (MQ=255)                                          tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGat                                                                                  >  2:701232/1‑139 (MQ=255)                                           aaaGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGc                                                                                                        >  1:153402/1‑116 (MQ=255)                                           aaaGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGc                                                                                                        <  2:153402/116‑1 (MQ=255)                                             aGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatgatg                                                                                                         <  2:661475/113‑1 (MQ=255)                                             aGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatgatg                                                                                                         >  1:661475/1‑113 (MQ=255)                                                    gCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaa                                                                        <  2:1371809/139‑1 (MQ=255)                                                         aaaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACT‑GCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCcg                                                                                                      >  1:835424/1‑103 (MQ=255)                                                         aaaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACT‑GCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCcg                                                                                                      <  2:835424/103‑1 (MQ=255)                                                          aaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                  <  1:1401908/139‑1 (MQ=255)                                                          aaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                  <  2:1095937/139‑1 (MQ=255)                                                          aaCGCTGGCGAATATGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCAAACCACTATTTCTGGAGGAACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTATAGAAAACGcc                                                                  <  2:1000715/139‑1 (MQ=255)                                                                 gCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTa                                                           >  2:488018/1‑139 (MQ=255)                                                                    aaTCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTaaaa                                                        <  2:511040/139‑1 (MQ=255)                                                                      tCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAAcc                                                      <  2:191078/139‑1 (MQ=255)                                                                          ggATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGAcc                                                  <  1:488018/139‑1 (MQ=255)                                                                             tatTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                  <  2:28907/120‑1 (MQ=255)                                                                             tatTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                  >  1:28907/1‑120 (MQ=255)                                                                                tttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatgatg                                                                                                         <  2:334664/78‑1 (MQ=255)                                                                                tttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatgatg                                                                                                         >  1:334664/1‑78 (MQ=255)                                                                                tttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGt                                            >  1:143083/1‑139 (MQ=255)                                                                                 ttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                           <  2:143083/139‑1 (MQ=255)                                                                                    tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAg                                        >  2:1333942/1‑139 (MQ=255)                                                                                    tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAg                                        <  2:702051/139‑1 (MQ=255)                                                                                       tCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                               >  1:1380139/1‑49 (MQ=255)                                                                                       tCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                               <  2:1380139/49‑1 (MQ=255)                                                                                       tCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCg                                     <  1:1333942/139‑1 (MQ=255)                                                                                              gCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATa                              <  1:701232/139‑1 (MQ=255)                                                                                                 cGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTa                                                                             >  2:374377/1‑89 (MQ=255)                                                                                                 cGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTa                                                                             <  1:374377/89‑1 (MQ=255)                                                                                                         cTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGt                                                                                             <  1:1096473/65‑1 (MQ=255)                                                                                                         cTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGt                                                                                             >  2:1096473/1‑65 (MQ=255)                                                                                                            cACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTaaa                                                                                                >  2:1197515/1‑59 (MQ=255)                                                                                                            cACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTaaa                                                                                                <  1:1197515/59‑1 (MQ=255)                                                                                                                ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATa                              >  1:1429724/1‑121 (MQ=255)                                                                                                                ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATa                              <  2:1429724/121‑1 (MQ=255)                                                                                                                ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgc            >  2:488643/1‑139 (MQ=255)                                                                                                                ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgc            >  1:96674/1‑139 (MQ=255)                                                                                                                  cacaCCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgcgc          <  1:488643/139‑1 (MQ=255)                                                                                                                     accacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCta       <  1:160860/139‑1 (MQ=255)                                                                                                                      ccacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCcgc                                                                                                     >  1:516943/1‑44 (MQ=255)                                                                                                                      ccacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCcgc                                                                                                     <  2:516943/44‑1 (MQ=255)                                                                                                                      ccacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACg                                                                                              >  1:227507/1‑51 (MQ=255)                                                                                                                      ccacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACg                                                                                              <  2:227507/51‑1 (MQ=255)                                                                                                                      ccacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGCAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTAc                                                             >  2:352153/1‑84 (MQ=255)                                                                                                                      ccacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGCAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTAc                                                             <  1:352153/84‑1 (MQ=255)                                                                                                                       cacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAAc                                                       <  1:231854/89‑1 (MQ=255)                                                                                                                       cacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAAc                                                       >  2:231854/1‑89 (MQ=255)                                                                                                                       cacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCtatt     <  2:96674/139‑1 (MQ=255)                                                                                                                          tATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATtata  <  2:346766/139‑1 (MQ=255)                                                                                                                           aTTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTTTACATCGGAAgc            >  1:1030345/1‑128 (MQ=255)                                                                                                                           aTTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgc            <  2:1030345/128‑1 (MQ=255)                                                                                                                                    |                                                                                                                                 GCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATA  >  NC_000913/1491258‑1491517

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