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%	258.0 / NA	70	K163R (AAG→AGG)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/39);  total (31/39)

CGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGC  >  NC_000913/1491254‑1491509                                                                                                                                        |                                                                                                                         cGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAgggc                                                                                                                       >  2:954720/1‑138 (MQ=255) cGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                       <  2:966183/139‑1 (MQ=255)    tGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGacaa                                                                                                                    <  2:215159/139‑1 (MQ=255)     gCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGacaac                                                                                                                   >  2:1196602/1‑139 (MQ=255)     gCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGacaac                                                                                                                   >  2:1333736/1‑139 (MQ=255)          aCGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACt                                                                                                              <  1:1196602/139‑1 (MQ=255)           cGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACtt                                                                                                             <  1:1263442/139‑1 (MQ=255)              ttCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTa                                                                                                          <  2:844222/139‑1 (MQ=255)                cAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACg                                                                                                        >  2:451105/1‑139 (MQ=255)                  gagGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCa                                                                                                      <  2:376593/139‑1 (MQ=255)                   agGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCat                                                                                                     >  2:1652513/1‑139 (MQ=255)                   agGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCat                                                                                                     >  2:262478/1‑139 (MQ=255)                   agGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCat                                                                                                     >  2:800454/1‑139 (MQ=255)                     gCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatga                                                                                                   >  2:1722113/1‑139 (MQ=255)                       cAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatgatg                                                                                                 <  1:262478/139‑1 (MQ=255)                        aaTTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCa                                                                                                      <  1:1729386/133‑1 (MQ=255)                        aaTTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCa                                                                                                      >  2:1729386/1‑133 (MQ=255)                          tttCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacac                                                                                                                 >  2:734745/1‑120 (MQ=255)                          tttCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacac                                                                                                                 <  1:734745/120‑1 (MQ=255)                                 cGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAAc                                                                                       <  1:1722113/139‑1 (MQ=255)                                  gggTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACg                                                                                      <  1:1652513/139‑1 (MQ=255)                                  gggTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACg                                                                                      <  1:1220766/139‑1 (MQ=255)                                            cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacac                                                                                                                 >  1:342356/1‑102 (MQ=255)                                            cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacac                                                                                                                 <  2:342356/102‑1 (MQ=255)                                             ttAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                       <  1:1696178/95‑1 (MQ=255)                                             ttAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                       >  2:96291/1‑95 (MQ=255)                                             ttAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                       <  1:96291/95‑1 (MQ=255)                                             ttAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa                                                                                                                       >  2:1696178/1‑95 (MQ=255)                                              tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTc                                                                                   >  1:397215/1‑130 (MQ=255)                                              tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTc                                                                                   <  2:397215/130‑1 (MQ=255)                                              tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGat                                                                          >  2:725319/1‑139 (MQ=255)                                                        gCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaa                                                                <  1:451105/139‑1 (MQ=255)                                                             aaaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAg                                                                             >  1:849874/1‑121 (MQ=255)                                                             aaaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAg                                                                             <  2:849874/121‑1 (MQ=255)                                                             aaaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGc                                                           >  1:708499/1‑139 (MQ=255)                                                                        aaTCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTaaaa                                                >  1:1723822/1‑139 (MQ=255)                                                                              ggttatTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTATAGAAAACGCCTGTACTAAAACCGAcc                                          <  2:1723822/136‑1 (MQ=255)                                                                                atatTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCg                                        <  1:800454/139‑1 (MQ=255)                                                                                atatTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCg                                        <  1:988669/139‑1 (MQ=255)                                                                                     ttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                   <  1:214517/139‑1 (MQ=255)                                                                                      tGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                   >  2:1492500/1‑138 (MQ=255)                                                                                      tGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                   <  1:1492500/138‑1 (MQ=255)                                                                                      tGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTAc                                  <  1:1269764/139‑1 (MQ=255)                                                                                        tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGa                                                                               <  2:1400012/92‑1 (MQ=255)                                                                                        tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGa                                                                               >  1:1400012/1‑92 (MQ=255)                                                                                        tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                          <  2:1444148/113‑1 (MQ=255)                                                                                        tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                          >  1:1444148/1‑113 (MQ=255)                                                                                           tCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCg                             <  1:1333736/139‑1 (MQ=255)                                                                                                        tGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTc                <  1:725319/139‑1 (MQ=255)                                                                                                             cTGCACTGGCACACCACTATTTCTGGAGGGacaaca                                                                                                                  <  1:445585/36‑1 (MQ=255)                                                                                                             cTGCACTGGCACACCACTATTTCTGGAGGGacaaca                                                                                                                  >  2:445585/1‑36 (MQ=255)                                                                                                             cTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                          >  1:1233465/1‑92 (MQ=255)                                                                                                             cTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                          <  2:1233465/92‑1 (MQ=255)                                                                                                                cACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCgg        >  2:1643302/1‑139 (MQ=255)                                                                                                                 aCTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGa       <  2:708499/139‑1 (MQ=255)                                                                                                                  cTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACgg                   <  2:50679/126‑1 (MQ=255)                                                                                                                  cTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACgg                   >  1:50679/1‑126 (MQ=255)                                                                                                                    ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaaa                                                               >  1:797429/1‑80 (MQ=255)                                                                                                                    ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaaa                                                               <  2:797429/80‑1 (MQ=255)                                                                                                                    ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTc                <  1:1186881/127‑1 (MQ=255)                                                                                                                    ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTc                >  2:1186881/1‑127 (MQ=255)                                                                                                                     gCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                          >  1:869862/1‑84 (MQ=255)                                                                                                                     gCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                          <  2:869862/84‑1 (MQ=255)                                                                                                                           cacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAg                                                                             <  2:1172454/59‑1 (MQ=255)                                                                                                                           cacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAg                                                                             >  1:1172454/1‑59 (MQ=255)                                                                                                                                   cTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTGAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgcgc  >  2:1241252/1‑126 (MQ=255)                                                                                                                                   cTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgcgc  <  1:1241252/126‑1 (MQ=255)                                                                                                                                     ggAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCt                                                                                 >  2:1484523/1‑45 (MQ=255)                                                                                                                                     ggAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCt                                                                                 <  1:1484523/45‑1 (MQ=255)                                                                                                                                     ggAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                          >  2:463588/1‑68 (MQ=255)                                                                                                                                     ggAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                          <  1:463588/68‑1 (MQ=255)                                                                                                                                        |                                                                                                                         CGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGC  >  NC_000913/1491254‑1491509

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