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

CCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATAT  >  NC_000913/1491259‑1491518                                                                                                                                   |                                                                                                                                  ccTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTTCACTGGCACACCACTATTTCTGGAGGGacaaca                                                                                                                           <  2:201422/139‑1 (MQ=255) ccTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGacaaca                                                                                                                           <  1:567881/139‑1 (MQ=255)                   aaTTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATGTTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatgatg                                                                                                          >  1:344509/1‑138 (MQ=255)                   aaTTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATGTTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatgatg                                                                                                          <  2:344509/138‑1 (MQ=255)                   aaTTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGc                                                                                                         >  2:416586/1‑139 (MQ=255)                               gTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGtt                                                                                             >  1:523872/1‑139 (MQ=255)                                      gCTTATAGTCGTGTCATGAAACGCTGGCGAATCTGTGATATTTTGTCATCGGGTTGCACGGTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCGACGCATGATGCCGCGTAAACGTTCCTGAAg                                                                                      <  2:782518/139‑1 (MQ=255)                                         tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGTGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGTGGGACAactca                                                                                                                         <  2:486205/101‑1 (MQ=255)                                         tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacaca                                                                                                                         >  1:486205/1‑101 (MQ=255)                                                gTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTaaa                                                                            >  1:810229/1‑139 (MQ=255)                                                 tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatgatg                                                                                                          >  1:550504/1‑108 (MQ=255)                                                 tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCatgatg                                                                                                          <  2:550504/108‑1 (MQ=255)                                                         aaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                   <  2:523872/139‑1 (MQ=255)                                                             cTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                   <  2:493100/135‑1 (MQ=255)                                                             cTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                   >  1:493100/1‑135 (MQ=255)                                                                      cTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCg                                                      <  2:810229/139‑1 (MQ=255)                                                                         ggATATTTTGTCATCGGGTTGCACGCTGCTGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGAcc                                                   <  2:161128/139‑1 (MQ=255)                                                                                   tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCt                                                                                                                    >  2:421023/1‑64 (MQ=255)                                                                                   tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCt                                                                                                                    <  1:421023/64‑1 (MQ=255)                                                                                     aTCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGacaaca                                                                                                                           <  2:652935/55‑1 (MQ=255)                                                                                     aTCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGacaaca                                                                                                                           >  1:652935/1‑55 (MQ=255)                                                                                        gggTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCgaa                                    >  1:402334/1‑139 (MQ=255)                                                                                                   tGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATa                               >  1:27932/1‑133 (MQ=255)                                                                                                   tGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATa                               <  2:27932/133‑1 (MQ=255)                                                                                                    gCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTc                         >  1:682145/1‑138 (MQ=255)                                                                                                    gCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTc                         <  2:682145/138‑1 (MQ=255)                                                                                                               ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgc             >  1:678730/1‑139 (MQ=255)                                                                                                                      cacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTc                                                                                            >  2:329038/1‑53 (MQ=255)                                                                                                                      cacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                            <  2:271798/101‑1 (MQ=255)                                                                                                                      cacTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                            >  1:271798/1‑101 (MQ=255)                                                                                                                      cacTATTGCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTc                                                                                            <  1:329038/53‑1 (MQ=255)                                                                                                                         tATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATtata   <  2:97512/139‑1 (MQ=255)                                                                                                                          aTTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTatat  >  2:676297/1‑139 (MQ=255)                                                                                                                                   |                                                                                                                                  CCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATAT  >  NC_000913/1491259‑1491518

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