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

TGGATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGACCGC  >  NC_000913/915599‑915855                                                                                                                                    |                                                                                                                              tGGATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCAc                                                                                                                        <  2:1359834/139‑1 (MQ=255)  ggATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                       <  1:441327/139‑1 (MQ=255)    aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGTAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                     <  2:1094261/139‑1 (MQ=255)    aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                     <  2:389801/139‑1 (MQ=255)             aGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAAt                                                                                                            >  1:1158697/1‑139 (MQ=255)                  aCCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                  >  2:661353/1‑128 (MQ=255)                  aCCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                  <  1:661353/128‑1 (MQ=255)                  aCCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATAcagc                                                                                                       >  1:667497/1‑139 (MQ=255)                   ccTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATAcagca                                                                                                      >  2:1225795/1‑139 (MQ=255)                        aCTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgcc                                                                                                 <  2:148263/139‑1 (MQ=255)                          tCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               <  2:1076551/139‑1 (MQ=255)                          tCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               <  1:1225795/139‑1 (MQ=255)                                 cAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGcaaca                                                                                        <  1:893645/139‑1 (MQ=255)                                    ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATa                                                                                     >  1:1037315/1‑139 (MQ=255)                                      cAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                     <  1:1482300/105‑1 (MQ=255)                                      cAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTACCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                     >  2:1482300/1‑105 (MQ=255)                                        aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                  >  1:474032/1‑106 (MQ=255)                                        aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCaa                                                                                                                  <  2:474032/106‑1 (MQ=255)                                        aCTACCAGCGCGGAAAGCATGGAATAA‑‑‑CCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               >  2:1253464/1‑122 (MQ=255)                                        aCTACCAGCGCGGAAAGCATGGAATAA‑‑‑CCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               <  1:1253464/122‑1 (MQ=255)                                               gcgcGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCACTGCAACAATACCGCCGACAAc                                                                          <  2:667497/139‑1 (MQ=255)                                                cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAcc                                                                         <  2:1158697/139‑1 (MQ=255)                                                  cgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTg                                                                       <  2:842479/139‑1 (MQ=255)                                                    gAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAGACCGCTGGCGGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTtgg                                                                     <  2:111308/139‑1 (MQ=255)                                                      aaGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCa                                                                   <  2:1037315/139‑1 (MQ=255)                                                        gCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAAt                                                                 >  2:1427692/1‑139 (MQ=255)                                                           tGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               <  2:312765/106‑1 (MQ=255)                                                           tGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               <  2:1486566/106‑1 (MQ=255)                                                           tGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               >  1:1486566/1‑106 (MQ=255)                                                           tGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               >  1:312765/1‑106 (MQ=255)                                                           tGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTC‑AAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACg                                                             >  1:368960/1‑139 (MQ=255)                                                                      ccTGGTGAATGCTCGCCATAACCGTTAGAAGCAAATCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGa                                                   <  1:838516/139‑1 (MQ=255)                                                                        tGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGAct                                                 >  1:215856/1‑139 (MQ=255)                                                                                    gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATAcagcag                                                                                                     <  2:72059/75‑1 (MQ=255)                                                                                    gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATAcagcag                                                                                                     >  1:72059/1‑75 (MQ=255)                                                                                    gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               <  2:1296278/81‑1 (MQ=255)                                                                                    gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                                               >  1:1296278/1‑81 (MQ=255)                                                                                     ccATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                     >  1:883098/1‑58 (MQ=255)                                                                                     ccATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                                     <  2:883098/58‑1 (MQ=255)                                                                                         aaCCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACg                                >  1:33635/1‑139 (MQ=255)                                                                                              ttAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                       <  2:477888/47‑1 (MQ=255)                                                                                              ttAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                                       >  1:477888/1‑47 (MQ=255)                                                                                              ttAGAAGCAAAACCGCTGGCTGCCGCGTC‑AAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGcc                          <  2:368960/139‑1 (MQ=255)                                                                                               tAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGcc                          <  2:215856/139‑1 (MQ=255)                                                                                               tAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGcc                          <  2:33635/139‑1 (MQ=255)                                                                                                     cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa                    >  1:162261/1‑139 (MQ=255)                                                                                                        aaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGc                                                                                            <  1:887786/64‑1 (MQ=255)                                                                                                        aaCCGCTGGCTGCCGCGTCAAAACCCGTCCTACCACTGGCAATTAAATACAGTAGCGCCGCTGc                                                                                            >  2:887786/1‑64 (MQ=255)                                                                                                               ggctgCCGCGTCAAAACCCGTCTTACCACTAGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAg                                                          >  1:782707/1‑91 (MQ=255)                                                                                                               ggctgCCGCGTCAAAACCCGTCTTACCACTAGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAg                                                          <  2:782707/91‑1 (MQ=255)                                                                                                                       cgTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGAccgc  <  1:1225828/139‑1 (MQ=255)                                                                                                                          cAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccg                                                                                                >  2:1446095/1‑42 (MQ=255)                                                                                                                          cAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccg                                                                                                <  1:1446095/42‑1 (MQ=255)                                                                                                                                    |                                                                                                                              TGGATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCATAAACCAATAGTGACCGC  >  NC_000913/915599‑915855

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