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

GGATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAA  >  NC_000913/915600‑915845                                                                                                                                   |                                                                                                                      ggATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                               <  2:1325512/139‑1 (MQ=255)   aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAGCCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCGTACCACTgg                                                                                                             <  1:1141898/139‑1 (MQ=255)   aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                             <  2:914644/139‑1 (MQ=255)   aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                             <  1:407670/139‑1 (MQ=255)                     gCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAgcg                                                                                         <  2:1060315/139‑1 (MQ=255)                         tCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGcgccgc                                                                                     <  2:1392307/139‑1 (MQ=255)                                   ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                             <  2:43384/107‑1 (MQ=255)                                   ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                             >  1:43384/1‑107 (MQ=255)                                   ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATAcagcag                                                                                           >  1:1113321/1‑123 (MQ=255)                                   ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATAcagcag                                                                                           <  2:1113321/123‑1 (MQ=255)                                   ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATa                                                                           >  1:1485490/1‑139 (MQ=255)                                   ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATa                                                                           >  1:421194/1‑139 (MQ=255)                                    tCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATAc                                                                          <  1:106612/139‑1 (MQ=255)                                     cAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                             >  1:1051241/1‑105 (MQ=255)                                     cAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                             <  2:1051241/105‑1 (MQ=255)                                     cAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGcaa                                                                                <  2:16453/132‑1 (MQ=255)                                     cAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGcaa                                                                                >  1:16453/1‑132 (MQ=255)                                       aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACagcagc                                                                                          >  1:292525/1‑120 (MQ=255)                                       aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACagcagc                                                                                          <  2:292525/120‑1 (MQ=255)                                          aCCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGa                                                                    <  2:1485490/139‑1 (MQ=255)                                            cAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                             >  2:1263515/1‑98 (MQ=255)                                            cAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                             <  1:1263515/98‑1 (MQ=255)                                             aGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACaa                                                                 >  2:202623/1‑139 (MQ=255)                                                 cgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTg                                                             <  1:202623/139‑1 (MQ=255)                                                           ggAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACg                                                   <  2:421194/139‑1 (MQ=255)                                                                  ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCaa                                            >  2:579625/1‑139 (MQ=255)                                                                   cGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTAcc                                                                                                                  >  1:1349927/1‑70 (MQ=255)                                                                   cGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTAcc                                                                                                                  <  2:1349927/70‑1 (MQ=255)                                                                      cTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCaa                                            >  1:69802/1‑135 (MQ=255)                                                                      cTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCaa                                            <  2:69802/135‑1 (MQ=255)                                                                      cTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGAc                                        >  2:262949/1‑139 (MQ=255)                                                                       tagTGAATGCTCGCCATAACCGATCGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCGTACCACTGGCAAAAA‑AAATACAGCAGCGCCGCTGCAACAATACCGCAGACAACCTGGGCAATTACGTAGCCAACGAc                                        >  2:205600/3‑139 (MQ=255)                                                                         gTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATAccgccg                                                                     >  1:1579611/1‑107 (MQ=255)                                                                         gTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATAccgccg                                                                     <  2:1579611/107‑1 (MQ=255)                                                                           gAATGCTCGCCATAACCGTTAGAAGCAATGCCGCTGGCTGCCGCGTCAAACCCCGTCTTACCACTGAC‑AATTAAAATACAGCAGCGCCGCTGCAACAATACCGCCGTCAACCTGGTTAATTACGTAGCCAACGActtct                                    <  1:205600/139‑1 (MQ=255)                                                                            aaTGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AAtt                                                                                                       >  1:845531/1‑71 (MQ=255)                                                                            aaTGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AAtt                                                                                                       <  2:845531/71‑1 (MQ=255)                                                                                cTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                               <  2:1312590/60‑1 (MQ=255)                                                                                cTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                               >  1:1312590/1‑60 (MQ=255)                                                                                  cGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCgg                            <  1:1275720/139‑1 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAcc                                                               <  2:1457676/103‑1 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAcc                                                               <  2:463203/103‑1 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAcc                                                               >  1:463203/1‑103 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAcc                                                               >  1:1457676/1‑103 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAc                                                                <  2:1465558/102‑1 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAc                                                                >  1:1465558/1‑102 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGActtctt                                   <  1:984827/131‑1 (MQ=255)                                                                                   gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGActtctt                                   >  2:984827/1‑131 (MQ=255)                                                                                                   gCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAgcg                                                                                         >  2:462253/1‑61 (MQ=255)                                                                                                   gCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAgcg                                                                                         <  1:462253/61‑1 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGa                                                                    >  1:5799/1‑81 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGa                                                                    <  2:5799/81‑1 (MQ=255)                                                                                                    cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa          <  1:262949/139‑1 (MQ=255)                                                                                                          cGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTg                                                             <  2:42294/82‑1 (MQ=255)                                                                                                          cGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTg                                                             >  1:42294/1‑82 (MQ=255)                                                                                                            ctggctgCCGCGTCAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCaa  >  2:29467/1‑139 (MQ=255)                                                                                                                         cAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAgcg                                                                                         >  1:632135/1‑39 (MQ=255)                                                                                                                         cAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAgcg                                                                                         <  2:632135/39‑1 (MQ=255)                                                                                                                         cAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGAc                                                                   <  1:1691632/61‑1 (MQ=255)                                                                                                                         cAAAACCCGTCTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGAc                                                                   >  2:1691632/1‑61 (MQ=255)                                                                                                                                   |                                                                                                                      GGATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGC‑AATT‑AAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAA  >  NC_000913/915600‑915845

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