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

ATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAAT  >  NC_000913/915602‑915846                                                                                                                                 |                                                                                                                     aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                            <  1:612370/139‑1 (MQ=255) aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                            <  1:1512653/139‑1 (MQ=255)                    cACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAACAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAgcgc                                                                                         <  1:30174/139‑1 (MQ=255)                          aTACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGc                                                                                   <  1:968873/139‑1 (MQ=255)                              cAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGcaaca                                                                               >  1:156709/1‑139 (MQ=255)                                 ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                              >  2:1591668/1‑105 (MQ=255)                                 ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                              <  1:1591668/105‑1 (MQ=255)                                 ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATa                                                                            >  2:1476925/1‑139 (MQ=255)                                 ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATa                                                                            >  2:61378/1‑139 (MQ=255)                                   cAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAcc                                                                          >  1:1146915/1‑139 (MQ=255)                                       tACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAccgccg                                                                      >  1:383311/1‑139 (MQ=255)                                             cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgcc                                                                                        <  2:60925/115‑1 (MQ=255)                                             cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgcc                                                                                        >  1:60925/1‑115 (MQ=255)                                                   aaGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCa                                                          <  2:1146915/139‑1 (MQ=255)                                                    aGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCaa                                                         >  2:313345/1‑139 (MQ=255)                                                       aTGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGTAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTa                                                      >  1:1393310/1‑139 (MQ=255)                                                       aTGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGACTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                            >  1:1286247/1‑85 (MQ=255)                                                       aTGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGACTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                            <  2:1286247/85‑1 (MQ=255)                                                                ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAgcgc                                                                                         <  2:820570/95‑1 (MQ=255)                                                                ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAgcgc                                                                                         >  1:820570/1‑95 (MQ=255)                                                                     tGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                            >  1:1409087/1‑71 (MQ=255)                                                                     tGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                            <  2:1409087/71‑1 (MQ=255)                                                                     tGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGAct                                        <  1:61378/139‑1 (MQ=255)                                                                                 gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTg                                                              <  1:281662/105‑1 (MQ=255)                                                                                 gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTg                                                              >  2:281662/1‑105 (MQ=255)                                                                                 gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGa                            >  1:316479/1‑139 (MQ=255)                                                                                      aaCCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGcaa                                                                                 >  2:467112/1‑81 (MQ=255)                                                                                      aaCCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGcaa                                                                                 <  1:467112/81‑1 (MQ=255)                                                                                      aaCCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACg                       >  1:757287/1‑139 (MQ=255)                                                                                            tAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGcc                 <  2:757287/139‑1 (MQ=255)                                                                                             aGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCaa                                             <  1:1606239/110‑1 (MQ=255)                                                                                             aGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCaa                                             >  2:1606239/1‑110 (MQ=255)                                                                                             aGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCa                <  2:383311/139‑1 (MQ=255)                                                                                                   aaaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                              <  1:1300580/39‑1 (MQ=255)                                                                                                   aaaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                              >  2:1300580/1‑39 (MQ=255)                                                                                                    aaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                            >  2:996028/1‑40 (MQ=255)                                                                                                    aaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                            <  1:996028/40‑1 (MQ=255)                                                                                                         gctggctgCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCa    >  2:7007/1‑139 (MQ=255)                                                                                                          ctggctgCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGcaaca                                                                               <  1:1160618/63‑1 (MQ=255)                                                                                                          ctggctgCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCaa   >  1:1210531/1‑139 (MQ=255)                                                                                                          ctggctgCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGATGcaaca                                                                               >  2:1160618/1‑63 (MQ=255)                                                                                                           tggctgCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAAt  >  2:762270/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                     ATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAAT  >  NC_000913/915602‑915846

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