Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,866,379	0	T	G	34.0%	82.2 / 12.3	50	N188T (AAC→ACC)	aqpZ	Aquaporin Z
Reads supporting (aligned to +/- strand):  ref base T (4/29);  new base G (16/1);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 1.49e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGC  >  ACB143/4866247‑4866509                                                                                                                                     |                                                                                                                                    ggTTCTTTCGCCAGGCCCCGATACCGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGa                                                                                                                                  <  2:146380/136‑1 (MQ=255) ggTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGa                                                                                                                                  <  3:230225/136‑1 (MQ=255)   ttCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCTGCGACGCGCGGGGGTGGCa                                                                                                                                >  5:309123/1‑136 (MQ=255)        tCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGCGACGCCCGGGGGTGACAGACGt                                                                                                                           >  1:298022/1‑136 (MQ=255)                              gccgccGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGGCCCCGGCTTCGCGCGGGGGTGTCAGAGGTTTTGGTGACCCGGGTCgagct                                                                                                      >  7:158776/1‑133 (MQ=255)                              gccgccGATACCAGCCCCGATCAGCGGCGCCAACCAGAACACCCCCCGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATc                                                                                                     <  5:231381/136‑1 (MQ=255)                                cgccgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGCGGCGCGCGGGGGTGTCCGAGGGTTTGGGGACCGGGCTCGCGATCAg                                                                                                   >  7:225066/1‑136 (MQ=255)                                   cgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTg                                                                                                <  7:84689/136‑1 (MQ=255)                                   cgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTg                                                                                                <  1:269150/136‑1 (MQ=255)                                   cgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTg                                                                                                <  8:151272/136‑1 (MQ=255)                                   cgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTg                                                                                                <  1:84257/136‑1 (MQ=255)                                   cgATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTg                                                                                                <  6:154103/136‑1 (MQ=255)                                    gATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTgg                                                                                               <  2:232923/136‑1 (MQ=255)                                    gATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTgg                                                                                               <  6:318767/136‑1 (MQ=255)                                     aTAGCAGCACCGAGCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGa                                                                                              <  1:154839/136‑1 (MQ=255)                                         cagcaCCGATCAGCTGCGCCACCCAGATCAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcag                                                                                          <  8:123041/136‑1 (MQ=255)                                         cagcaCCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGTTATTGGTCACCGGGATCGAGATCAGGTGGAtcag                                                                                          <  2:258696/136‑1 (MQ=255)                                         cagcaCCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTCTTGGTCACCGGGATCGAGATCAGGTGGAtcag                                                                                          <  5:282901/136‑1 (MQ=255)                                         cagcaCCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcag                                                                                          <  4:151482/136‑1 (MQ=255)                                           gcaCCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGAtcaggg                                                                                        <  4:175275/136‑1 (MQ=255)                                               cGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGGGCGACGCGCGGGGGTGGCAGAGGGATTGGTGACCCGGGTCGAGCTCAGGTGGGtgagggt‑cag                                                                                   >  5:109204/1‑136 (MQ=255)                                               cGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGCGACGCCCGGGGGGGTCAGAGGTATTTGTCACCGGGCTCGAGATCAGGGGGCtcaggggccg                                                                                    >  7:226933/1‑132 (MQ=255)                                                    aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGGTGTCAGAGGTTTTTGTGACCGGGCTCGGGATGAGGTGGGTCAGGGT‑CAGGGCCa                                                                              >  7:45597/1‑136 (MQ=255)                                                    aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCGACGCGCGGGGGTGGCAGACGTATTGGTGACCGGGCTCGAGATCAGGGGGATCAGGGT‑GAGGGccc                                                                              >  2:29608/1‑135 (MQ=255)                                                    aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGGGACCCCCGGGGGTGGCACACGTATTTGTGGACGGGCTAGAGATAAGGGGGGTAAGGGT‑GAGCGCCa                                                                              >  7:176172/1‑136 (MQ=255)                                                    aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGCTACGCGCGGGGGTGACCGAGGTTTTGGTGACCGGGCTGGAGATCAGGGGGATCAGGGT‑CAGGGGCa                                                                              >  6:338696/1‑136 (MQ=255)                                                    aGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGCGACCCGCGGGGTTGGCCGAGGGATTGGTCACCGGGAACGAGCTGACCTGGTTCACGGT‑CACGGccc                                                                              >  8:190968/1‑135 (MQ=255)                                                         cgcCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGGCCCGGGGGGCGCGCGGGGGTGACAGCGGGATTTGGGCCCGGGGTCGAGATCAGGGGGGTGAGGGG‑CACGGCCACGccc                                                                         >  8:198825/1‑135 (MQ=255)                                                              cccAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCcgatggc                                                                    <  7:154910/136‑1 (MQ=255)                                                              cccAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCcgatggc                                                                    <  3:826/136‑1 (MQ=255)                                                              cccAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCGGCTACGCGCGGGGGTGTCAGAGGTATTTGTGACCGGGGTGGAGATCACGTGGGTCAGGGT‑GAGGGCCCGGCcgcaggg                                                                    >  6:172835/1‑131 (MQ=255)                                                               ccAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCcgatggcg                                                                   <  7:108159/136‑1 (MQ=255)                                                                           caGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGGGGGACGCGCGGGGGTGTCAGAGGTTTTGGGCCCCGGGGTCGAGCTCAGGTGGCTGAGGGT‑CAGGGCCCGGCCGCTGGGGATGGGTGCaaa                                                        >  1:198554/1‑135 (MQ=255)                                                                           caGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCCTGCGAAGCCCGGGGGTGACAGAGGTATTTGTCACCCGGGACGAGATCAAGTGGCTAACGGT‑CAGGGCCAGGCCGCTGGCGATGGGTGCaaa                                                        >  2:41344/1‑135 (MQ=255)                                                                                   gcAGGGCCCAGCCACCGACGAACAGGGCCGGGCCCCGGCTGCGCGCGGGGGGGTCAGAGGGTTTGGTGACCGGGGTGGGGATGAGGGGGGTGCGGGTGCGGGGCCCGG‑CGGTGGCGGTGGGTGGAAAACCTGCCgg                                               >  4:83972/1‑136 (MQ=255)                                                                                    cAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGt                                              <  1:12527/136‑1 (MQ=255)                                                                                               cACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCgtc                                   <  1:110417/136‑1 (MQ=255)                                                                                                 ccGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCgtcgg                                 <  2:11325/136‑1 (MQ=255)                                                                                                    acgaACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCgtcggtcg                              <  1:29608/136‑1 (MQ=255)                                                                                                            ggccggccCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCatg                      <  1:124696/136‑1 (MQ=255)                                                                                                             gccggccCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCatga                     <  4:40059/136‑1 (MQ=255)                                                                                                             gccggccCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCatga                     <  3:214515/136‑1 (MQ=255)                                                                                                              ccggccCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCatgat                    <  5:172835/136‑1 (MQ=255)                                                                                                              ccggccCCGTGCTACGCGCGGGGGTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCatgat                    <  3:23205/136‑1 (MQ=255)                                                                                                               cggccCCCTTCGACGCGCGGGGGTGAGATCGGTTTTTGTCTCCGGGATCGAGTTAAGGGGGATAACGGG‑GACGGCCCCGGCGCTGGCGATGGGTGCAAAACCAGCCCGTGCCCCGGCCGCGGGGGccccccagat                    >  6:325423/1‑130 (MQ=255)                                                                                                                  ccccGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGacc                <  2:313750/136‑1 (MQ=255)                                                                                                                               gctgggTTGACACACGTATTGGTCACCGGGAGCGAGATCAGGTGGAGCAGGGT‑CAGGGCCAGGCCGATGGCCATGGGTGCAAAACCTGCCGGAGCGCGGGCGTCGGTCGCCCCCATGATTACCCCCCGGAACATCg   >  6:229438/4‑136 (MQ=255)                                                                                                                               gcgGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCg   >  1:74723/1‑136 (MQ=255)                                                                                                                               gcgGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGGCGCCCCCATGATGACCACCAGGAACATCg   >  1:248537/1‑136 (MQ=255)                                                                                                                                cgGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGc  <  1:32325/136‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    GGTTCTTTCGCCAGGCCCCGATACAGCGCGCCGCCGATAGCAGCACCGATCAGCGGCGCCACCCAGAACAGCCACAGCTGCTGCAGGGCCCAGCCACCGACGAACAGGGCCGGCCCCGTGCTACGCGCGGGGTTGACAGAGGTATTGGTCACCGGGATCGAGATCAGGTGGATCAGGGT‑CAGGGCCAGGCCGATGGCGATGGGTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGC  >  ACB143/4866247‑4866509

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