Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,491,395	A→G	N165S (AAC→AGC)	cybB →	cytochrome b561

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,395	0	A	G	100.0%	162.8 / NA	45	N165S (AAC→AGC)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/24);  total (21/24)

GCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGA  >  NC_000913/1491265‑1491530                                                                                                                                   |                                                                                                                                        gCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCt                                                                                                                                <  1:983701/140‑1 (MQ=255) gCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCt                                                                                                                                <  1:400317/140‑1 (MQ=255)     cAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCt                                                                                                                                >  2:735744/1‑136 (MQ=255)     cAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCt                                                                                                                                <  1:735744/136‑1 (MQ=255)     cAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGc                                                                                                                            >  1:818610/1‑140 (MQ=255)     cAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGc                                                                                                                            >  1:224638/1‑140 (MQ=255)     cAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGc                                                                                                                            >  2:653268/1‑140 (MQ=255)          gCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTACACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCatgat                                                                                                                       >  1:792208/1‑140 (MQ=255)           ccAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCatgatg                                                                                                                      <  2:224638/140‑1 (MQ=255)           ccAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCatgatg                                                                                                                      <  1:465838/140‑1 (MQ=255)           ccAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTAATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCatgatg                                                                                                                      <  1:63743/140‑1 (MQ=255)               tttCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCcgc                                                                                                                  >  1:499158/1‑140 (MQ=255)                   gAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGGCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTaa                                                                                                              <  1:717491/140‑1 (MQ=255)                    aaCGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTaaa                                                                                                             >  2:582229/1‑140 (MQ=255)                    aaCGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTaaa                                                                                                             >  2:34712/1‑140 (MQ=255)                      cGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACg                                                                                                           <  1:34712/140‑1 (MQ=255)                                   tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACg                                                                                                           <  1:802906/127‑1 (MQ=255)                                   tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACg                                                                                                           >  2:802906/1‑127 (MQ=255)                                   tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGata                                                                                              >  2:422107/1‑140 (MQ=255)                                       gTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTa                                                                                          <  2:942492/140‑1 (MQ=255)                                       gTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATCTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTa                                                                                          <  1:108711/140‑1 (MQ=255)                                            ggCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaa                                                                                     <  1:732562/140‑1 (MQ=255)                                                  aaaCGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGcc                                                                               <  1:274134/140‑1 (MQ=255)                                                        tGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACt                                                                         >  2:981682/1‑140 (MQ=255)                                                             aaTCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAAc                                                                    >  2:174604/1‑140 (MQ=255)                                                               tCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCg                                                                  <  1:174604/140‑1 (MQ=255)                                                                    gATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCg                                                             >  1:214152/1‑140 (MQ=255)                                                                         tttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa                                                        <  2:214152/140‑1 (MQ=255)                                                                             tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAgg                                                    >  1:728660/1‑140 (MQ=255)                                                                             tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAgg                                                    >  2:771196/1‑140 (MQ=255)                                                                               aTCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCg                                                  <  2:332086/140‑1 (MQ=255)                                                                                      tGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGTCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATa                                           <  2:818610/140‑1 (MQ=255)                                                                                       gCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATAc                                          >  2:817889/1‑140 (MQ=255)                                                                                         aCGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACgg                                        <  1:817889/140‑1 (MQ=255)                                                                                             tGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCt                                    >  2:390556/1‑140 (MQ=255)                                                                                             tGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCt                                    >  1:220557/1‑140 (MQ=255)                                                                                                          gCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgcgc                       <  1:653268/140‑1 (MQ=255)                                                                                                          gCACACCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgcgc                       <  1:582229/140‑1 (MQ=255)                                                                                                           cacaCCACTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGcc                      <  2:499158/140‑1 (MQ=255)                                                                                                                  cTATTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATtata               <  2:728660/140‑1 (MQ=255)                                                                                                                    aTTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGTCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGcc                      <  2:32186/131‑1 (MQ=255)                                                                                                                    aTTTCTGGAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGTCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGcc                      >  1:32186/1‑131 (MQ=255)                                                                                                                          ggAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCta                    <  2:346745/127‑1 (MQ=255)                                                                                                                          ggAAGGACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCta                    >  1:346745/1‑127 (MQ=255)                                                                                                                               gACAGCACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGa  >  2:676674/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                        GCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATTATATTTATTTGTATGA  >  NC_000913/1491265‑1491530

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