Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,392	A→G	D164G (GAC→GGC)	cybB →	cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,392	0	A	G	100.0%	159.0 / NA	44	D164G (GAC→GGC)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/23);  total (21/23)

GACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGATAA  >  NC_000913/1491252‑1491533                                                                                                                                             |                                                                                                                                              gACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACtt                                                                                                                                     <  2:537282/151‑1 (MQ=255) gACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACt                                                                                                                                      <  2:636724/150‑1 (MQ=255)     aTGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTAc                                                                                                                                 <  2:180496/151‑1 (MQ=255)          ttACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCatg                                                                                                                            >  1:267844/1‑151 (MQ=255)             cGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCatgatg                                                                                                                         <  1:435541/151‑1 (MQ=255)                  cAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCcgc                                                                                                                     <  2:203195/150‑1 (MQ=255)                     agGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTa                                                                                                                  >  2:749727/1‑150 (MQ=255)                              tCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTcc                                                                                                          <  1:611578/149‑1 (MQ=255)                              tCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTg                                                                                                        >  2:774229/1‑151 (MQ=255)                                        aGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGTTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACAAACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATAttt                                                                                              >  2:771537/1‑151 (MQ=255)                                        aGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATAttt                                                                                              <  1:186502/151‑1 (MQ=255)                                        aGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATAttt                                                                                              >  2:420748/1‑151 (MQ=255)                                         gATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATAttt                                                                                              >  1:772221/1‑150 (MQ=255)                                           tAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTaa                                                                                            >  2:154245/1‑150 (MQ=255)                                                tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCCGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaa                                                                                        >  1:323445/1‑149 (MQ=255)                                                  aaGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACg                                                                                    >  2:214741/1‑151 (MQ=255)                                                   aGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGc                                                                                   >  1:288813/1‑151 (MQ=255)                                                         ggCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTa                                                                              <  2:94/150‑1 (MQ=255)                                                                     tGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGAccc                                                                 >  1:156018/1‑151 (MQ=255)                                                                     tGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGAccc                                                                 >  1:421583/1‑151 (MQ=255)                                                                           aTCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGt                                                            >  2:742924/1‑150 (MQ=255)                                                                               gggATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACa                                                         <  1:84396/149‑1 (MQ=255)                                                                                ggATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGc                                                      >  2:539100/1‑151 (MQ=255)                                                                                 gATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCg                                                     <  1:774153/151‑1 (MQ=255)                                                                                         gTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATa                                              <  2:493174/150‑1 (MQ=255)                                                                                            aTCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACggg                                          >  1:473034/1‑151 (MQ=255)                                                                                                ggTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTa                                      <  1:643690/151‑1 (MQ=255)                                                                                                 gTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTAc                                     >  1:611792/1‑151 (MQ=255)                                                                                                  ttGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTAc                                     >  2:336139/1‑150 (MQ=255)                                                                                                    gCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACAt                                   <  1:8650/150‑1 (MQ=255)                                                                                                       cGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGa                               <  1:539048/151‑1 (MQ=255)                                                                                                        gCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGaa                              <  1:572755/151‑1 (MQ=255)                                                                                                        gCTGCGGCAGCACTGGCACACCACTATTTCTGGAAGGGCAACACACGTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTTTACTAAAACCGACCCGTGGTCCAGGCGAAGAACACGGGTCTACATCGGaa                              <  1:741119/151‑1 (MQ=255)                                                                                                                tGCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCtat                      >  1:80226/1‑151 (MQ=255)                                                                                                                 gCACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCtatt                     <  1:344504/151‑1 (MQ=255)                                                                                                                  cACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCtatta                    <  2:662418/151‑1 (MQ=255)                                                                                                                  cACTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCtatta                    <  1:170920/151‑1 (MQ=255)                                                                                                                   aCTGGCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCtattat                   <  2:48916/151‑1 (MQ=255)                                                                                                                      ggCACACCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTAtatt                <  1:541213/151‑1 (MQ=255)                                                                                                                         acacCACTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAATCCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTAtatttat             >  2:281953/1‑151 (MQ=255)                                                                                                                          caccacTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTAtatttatt            >  2:699331/1‑151 (MQ=255)                                                                                                                           accacTATTTCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATatttattt           <  1:2628/151‑1 (MQ=255)                                                                                                                                    tCTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGata   <  1:194955/150‑1 (MQ=255)                                                                                                                                     cTGGAAGGGCAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGataa  >  2:167776/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                              GACGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGATAA  >  NC_000913/1491252‑1491533

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