Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	915,730	T→C	K106K (AAA→AAG)	aqpZ ←	aquaporin Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	915,730	0	T	C	100.0%	149.8 / NA	43	K106K (AAA→AAG)	aqpZ	aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/25);  total (18/25)

GGATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGACCG  >  NC_000913/915600‑915854                                                                                                                                   |                                                                                                                             ggATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                      <  2:174059/139‑1 (MQ=255)   aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGTGTCAAAACCCGTCTTACCACTgg                                                                                                                    <  1:137614/139‑1 (MQ=255)   aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                    <  1:191731/139‑1 (MQ=255)      aCCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                 <  1:47073/139‑1 (MQ=255)           cAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTaaa                                                                                                            <  1:105208/139‑1 (MQ=255)                    tGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACagcagc                                                                                                   <  2:741674/139‑1 (MQ=255)                             tACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                 >  2:241073/1‑116 (MQ=255)                             tACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                 <  1:241073/116‑1 (MQ=255)                                    tCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAc                                                                                   <  2:1127036/139‑1 (MQ=255)                                       aCTACCAGCGCGGAAAGAATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAccgc                                                                                >  2:530087/1‑139 (MQ=255)                                         tACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCTGCGCCGCTGCAACAATAccgccg                                                                              <  1:530087/139‑1 (MQ=255)                                          aCCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              >  2:275604/1‑122 (MQ=255)                                          aCCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              <  1:275604/122‑1 (MQ=255)                                             aGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACaa                                                                          >  2:1036326/1‑139 (MQ=255)                                               cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                      >  1:1180555/1‑93 (MQ=255)                                               cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                      <  2:1180555/93‑1 (MQ=255)                                               cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                    >  1:261293/1‑95 (MQ=255)                                               cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                    <  2:261293/95‑1 (MQ=255)                                               cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              >  1:724662/1‑117 (MQ=255)                                               cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              <  2:724662/117‑1 (MQ=255)                                                gcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCt                                                                       >  1:1241908/1‑139 (MQ=255)                                                 cgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAgcgc                                                                                                 >  2:626897/1‑112 (MQ=255)                                                 cgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAgcgc                                                                                                 <  1:626897/112‑1 (MQ=255)                                                   gAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTggg                                                                    <  1:1036326/139‑1 (MQ=255)                                                          tGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              <  2:295928/106‑1 (MQ=255)                                                          tGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              >  1:295928/1‑106 (MQ=255)                                                                     ccTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCAc                                                                                                                       <  2:799080/70‑1 (MQ=255)                                                                     ccTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCAc                                                                                                                       >  1:799080/1‑70 (MQ=255)                                                                     ccTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              >  1:263775/1‑95 (MQ=255)                                                                     ccTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              <  2:263775/95‑1 (MQ=255)                                                                                             ttAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGc                          >  1:1058117/1‑139 (MQ=255)                                                                                                  aGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGcc                         >  2:1005951/1‑135 (MQ=255)                                                                                                  aGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGcc                         <  1:1005951/135‑1 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTg                                                                      <  2:1075205/88‑1 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTg                                                                      >  1:1075205/1‑88 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa                   <  1:1303316/139‑1 (MQ=255)                                                                                                      aaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                 <  2:1343020/43‑1 (MQ=255)                                                                                                      aaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                 >  1:1343020/1‑43 (MQ=255)                                                                                                       aaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa                   <  2:1165905/136‑1 (MQ=255)                                                                                                       aaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa                   >  1:1165905/1‑136 (MQ=255)                                                                                                            ctggctgCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              >  1:668243/1‑56 (MQ=255)                                                                                                            ctggctgCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                              <  2:668243/56‑1 (MQ=255)                                                                                                                     gcgTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGAccg  <  1:1293288/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                             GGATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGACCG  >  NC_000913/915600‑915854

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