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

ATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCA  >  NC_000913/915602‑915837                                                                                                                                 |                                                                                                            aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAATCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                   <  1:772672/139‑1 (MQ=255) aTCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                   <  2:1552750/139‑1 (MQ=255)     ccAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAAt                                                                                               >  1:750932/1‑139 (MQ=255)                ccTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATAcagca                                                                                    >  1:819311/1‑139 (MQ=255)                ccTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATAcagca                                                                                    >  2:414470/1‑139 (MQ=255)                       tCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                             <  1:1504070/139‑1 (MQ=255)                                 ttCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATa                                                                   >  2:1478337/1‑139 (MQ=255)                                 ttCAACTACCAGCGCCGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                   >  1:1216521/1‑107 (MQ=255)                                 ttCAACTACCAGCGCCGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                   <  2:1216521/107‑1 (MQ=255)                                  tCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAc                                                                  <  2:750932/139‑1 (MQ=255)                                     aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                     <  1:1550370/101‑1 (MQ=255)                                     aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                     >  2:1550370/1‑101 (MQ=255)                                     aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                             <  1:163127/125‑1 (MQ=255)                                     aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                             >  2:1306097/1‑125 (MQ=255)                                     aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                             <  1:1306097/125‑1 (MQ=255)                                     aCTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                             >  2:163127/1‑125 (MQ=255)                                            gcgcGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAccg                                                                >  1:1416169/1‑131 (MQ=255)                                            gcgcGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATAccg                                                                <  2:1416169/131‑1 (MQ=255)                                            gcgcGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAc                                                        >  2:359848/1‑139 (MQ=255)                                             cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCAc                                                                                                      >  2:329401/1‑92 (MQ=255)                                             cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCAc                                                                                                      <  1:329401/92‑1 (MQ=255)                                             cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                             <  1:1535404/117‑1 (MQ=255)                                             cgcgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                             >  2:1535404/1‑117 (MQ=255)                                               cgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCATTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTg                                                     <  2:1520266/139‑1 (MQ=255)                                               cgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGc                                                                          <  1:1537463/118‑1 (MQ=255)                                               cgGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGc                                                                          >  2:1537463/1‑118 (MQ=255)                                                   aaGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCa                                                 <  1:359848/139‑1 (MQ=255)                                                                ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                   >  1:1312811/1‑76 (MQ=255)                                                                ccGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTgg                                                                                                   <  2:1312811/76‑1 (MQ=255)                                                                     tGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGAct                               >  1:262908/1‑139 (MQ=255)                                                                       gTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACGACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAcc                                                       >  1:877357/1‑113 (MQ=255)                                                                       gTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAAcc                                                       <  2:877357/113‑1 (MQ=255)                                                                          aaTGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccg                                                                              <  1:1008545/87‑1 (MQ=255)                                                                          aaTGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccg                                                                              >  2:1008545/1‑87 (MQ=255)                                                                                 gCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGa                   >  2:393353/1‑139 (MQ=255)                                                                                           ttAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGAc                                                           >  2:1273134/1‑89 (MQ=255)                                                                                           ttAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGAc                                                           <  1:1273134/89‑1 (MQ=255)                                                                                            tAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                     >  1:1308934/1‑46 (MQ=255)                                                                                            tAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACt                                                                                                     <  2:1308934/46‑1 (MQ=255)                                                                                             aGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGa                   <  2:17681/127‑1 (MQ=255)                                                                                             aGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGa                   >  1:17681/1‑127 (MQ=255)                                                                                             aGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCa       >  2:1596881/1‑139 (MQ=255)                                                                                                  cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                             >  2:70971/1‑64 (MQ=255)                                                                                                  cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGcgccgc                                                                             <  1:70971/64‑1 (MQ=255)                                                                                                  cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCa       >  1:391661/1‑134 (MQ=255)                                                                                                  cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCa       <  2:391661/134‑1 (MQ=255)                                                                                                  cAAAACCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCa  <  1:393353/139‑1 (MQ=255)                                                                                                    aaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGAc                                                           <  1:70969/80‑1 (MQ=255)                                                                                                    aaaCCGCTGGCTGCCGCGTCAAAACCCGTCTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGAc                                                           >  2:70969/1‑80 (MQ=255)                                                                                                                                 |                                                                                                            ATCACCAACAGGAAACCTGCACTCAATACCAGTTCAACTACCAGCGCGGAAAGCATGGAATAACCGCCTGGTGAATGCTCGCCATAACCGTTAGAAGCAAAACCGCTGGCTGCCGCGTCAAAACCCGTTTTACCACTGGCAATTAAATACAGCAGCGCCGCTGCAACAATACCGCCGACAACCTGGGCAATTACGTAGCCAACGACTTCTTTTGCCGGAAAACGTCCGCCAGCCCA  >  NC_000913/915602‑915837

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