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%	213.2 / NA	59	K106K (AAA→AAG)	aqpZ	aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/29);  total (30/29)

GGATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGACCGCCGGGTTAAA  >  NC_000913/915600‑915864                                                                                                                                   |                                                                                                                                       ggATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                                <  2:620949/139‑1 (MQ=255)        cAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                           <  1:1256657/137‑1 (MQ=255)        cAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                           >  2:1256657/1‑137 (MQ=255)             ggAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATa                                                                                                                    >  2:263863/1‑139 (MQ=255)             ggAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATa                                                                                                                    >  1:655606/1‑139 (MQ=255)              gAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                              >  2:1565987/1‑128 (MQ=255)              gAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                              <  1:1565987/128‑1 (MQ=255)                   cTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATAcagcag                                                                                                              <  2:1380381/139‑1 (MQ=255)                      cACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAgcgc                                                                                                           <  2:922657/139‑1 (MQ=255)                       aCTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgcc                                                                                                          <  2:665407/139‑1 (MQ=255)                         tCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                                        <  2:655606/139‑1 (MQ=255)                                   ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAgcgc                                                                                                           >  1:466964/1‑126 (MQ=255)                                   ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAgcgc                                                                                                           <  2:466964/126‑1 (MQ=255)                                     cAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                                <  1:1026733/103‑1 (MQ=255)                                     cAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                                >  2:1026733/1‑103 (MQ=255)                                       aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATAc                                                                                                                   >  2:446025/1‑114 (MQ=255)                                       aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATAc                                                                                                                   <  1:446025/114‑1 (MQ=255)                                       aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAccgc                                                                                          >  1:989495/1‑139 (MQ=255)                                       aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAccgc                                                                                          >  2:886036/1‑139 (MQ=255)                                               cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAcc                                                                                  <  1:263863/139‑1 (MQ=255)                                                         aTGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAc                                                                                   <  1:412008/128‑1 (MQ=255)                                                         aTGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAc                                                                                   >  2:412008/1‑128 (MQ=255)                                                         aTGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCa                                                                            >  2:1260897/1‑135 (MQ=255)                                                         aTGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCa                                                                            <  1:1260897/135‑1 (MQ=255)                                                         aTGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTa                                                                        >  2:1368509/1‑139 (MQ=255)                                                           ggAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                           >  1:670132/1‑86 (MQ=255)                                                           ggAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                           <  2:670132/86‑1 (MQ=255)                                                                aaCCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAc                                                                                             >  1:1493930/1‑111 (MQ=255)                                                                aaCCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAc                                                                                             <  2:1493930/111‑1 (MQ=255)                                                                aaCCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGcc                                                                 >  2:727940/1‑139 (MQ=255)                                                                 aCCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGTCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCa                                                                <  1:1609842/139‑1 (MQ=255)                                                                  ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                                >  1:1424215/1‑74 (MQ=255)                                                                  ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                                <  2:1424215/74‑1 (MQ=255)                                                                  ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGcaacaa                                                                                                <  2:542969/106‑1 (MQ=255)                                                                  ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGcaacaa                                                                                                >  1:542969/1‑106 (MQ=255)                                                                  ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCaa                                                               >  2:296704/1‑139 (MQ=255)                                                                  ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCaa                                                               >  1:1530990/1‑139 (MQ=255)                                                                     ccTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGa                                                            >  2:876620/1‑139 (MQ=255)                                                                     ccTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGa                                                            >  2:902866/1‑139 (MQ=255)                                                                       tGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGAct                                                          <  2:1530990/139‑1 (MQ=255)                                                                            aaTGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCttt                                                     >  2:1111311/1‑139 (MQ=255)                                                                             aTGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCtttt                                                    >  2:1438932/1‑139 (MQ=255)                                                                                   gcgaTAACCGTGAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTAACACTgg                                                                                                                              >  1:932789/4‑59 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                              <  2:932789/59‑1 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                           <  2:1528995/62‑1 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                           >  1:1528995/1‑62 (MQ=255)                                                                                               aGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCa                                  >  1:227013/1‑139 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa                             <  2:467969/139‑1 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa                             >  1:467969/1‑139 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa                             <  1:296704/139‑1 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa                             <  1:1111311/139‑1 (MQ=255)                                                                                                     aaaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCAC‑GGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATa                           >  1:258932/1‑139 (MQ=255)                                                                                                      aaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATa                           >  1:788503/1‑139 (MQ=255)                                                                                                         ccGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAAc                        <  1:876620/139‑1 (MQ=255)                                                                                                                     gcgTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGAccg            <  1:1438932/139‑1 (MQ=255)                                                                                                                      cgTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGAccgc           <  1:902866/139‑1 (MQ=255)                                                                                                                         cAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGACCGCCgg        >  1:192962/1‑139 (MQ=255)                                                                                                                               ccGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGACCGCCGGGTTaaa  <  2:227013/139‑1 (MQ=255)                                                                                                                               ccGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGACCGCCGGGTTaaa  <  2:192962/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       GGATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGACCGCCGGGTTAAA  >  NC_000913/915600‑915864

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