Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	5,960,590	G→T	9.5%	P313Q (CCG→CAG)	amtB_2 ←	Ammonia channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	5,960,590	0	G	T	9.5%	116.3 / 3.8	42	P313Q (CCG→CAG)	amtB_2	Ammonia channel
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base T (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CACGCCGAAGGCATCCAGCGAATCGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCGGCGAACATCCAGCCCAGGGCGGCGGCGGCGGTAGCGATCTGAGTGACCA  >  KT2440/5960460‑5960716                                                                                                                                   |                                                                                                                               cACGCCGAAGGCATCCAGCGAATCGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCgggg                                                                                                                             >  1:163185/1‑134 (MQ=255)  aCGCCGAAGGCATCCAGCGAATCGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGACACCAGGCCGATGCTCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGa                                                                                                                          >  8:271855/1‑136 (MQ=255)  aCGCCGAAGGCATCCAGCGAATCGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTg                                                                                                                           >  4:182051/1‑135 (MQ=255)  aCGCCGAAGGCATCCAGCGAATCGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGa                                                                                                                          >  7:112624/1‑136 (MQ=255)  aCGCCGAAGGCATCCAGCGAATCGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGa                                                                                                                          >  7:175774/1‑136 (MQ=255)       gAAGGCATCCAGCGAATCGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCg                                                                                                                     <  8:112624/136‑1 (MQ=255)       gAAGGCATCCAGCGAATCGTCATAGCCGAGTTGGCGCTTGAGCGTGGGGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGc                                                                                                                      <  2:211906/135‑1 (MQ=255)            cATCCAGCGAATCGTCATAGCCGAGTTGGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCa                                                                                                                 <  4:172874/135‑1 (MQ=255)                 aGCGAATCGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACCCCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCa                                                                                                            >  1:84609/1‑135 (MQ=255)                   cgatgcGTCATAGCCGAGTTTGCGCGTGAGCGTGGTGGCGCAGAAGTAGCACAGCACACCCTAAACCAGGCCGAGGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGc                                                                                                          <  2:12332/130‑1 (MQ=255)                    gAATCGTCATAGCCGAGTTTGCGCTTGCGCGTGGTGGCGCAGAAGTAGCAGCGCACCCCCGAAACCCGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCa                                                                                                         <  5:233076/135‑1 (MQ=255)                        cGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCa                                                                                                     >  2:211034/1‑135 (MQ=255)                                  gAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAACCCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGc                                                                                           >  8:203223/1‑135 (MQ=255)                                        gcgcTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGccc                                                                                    >  7:2471/1‑136 (MQ=255)                                          gcTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGccc                                                                                    <  4:36097/134‑1 (MQ=255)                                              gAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAACCCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCagag                                                                                <  7:243758/134‑1 (MQ=255)                                               aGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCagag                                                                                <  3:238581/133‑1 (MQ=255)                                               aGCGGGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCagag                                                                                <  3:105211/133‑1 (MQ=255)                                               aGCGGGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCagag                                                                                <  6:104680/133‑1 (MQ=255)                                                 cGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACt                                                                            >  4:19354/1‑135 (MQ=255)                                                  gtggtgGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACtt                                                                           <  3:104820/135‑1 (MQ=255)                                                     gtggCGCAGAAGGAGCAGAGCACACCCCAAACCAGGCCGATGATCAGCCCGCAAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCTGGCGAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGc                                                                        <  5:246494/135‑1 (MQ=255)                                                         cgcAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTAc                                                                    >  3:163819/1‑135 (MQ=255)                                                              aaGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTgac                                                              >  7:268574/1‑136 (MQ=255)                                                                           cacCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCgg                                                  >  2:28352/1‑135 (MQ=255)                                                                           cacCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCgg                                                  >  1:38807/1‑135 (MQ=255)                                                                           cacCCGAAACCAGGCCGATCATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCgg                                                  <  2:3831/135‑1 (MQ=255)                                                                            acCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGCCCCACTCGGc                                                 >  1:139379/1‑135 (MQ=255)                                                                               cGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCTGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCa                                                       >  4:194721/1‑126 (MQ=255)                                                                               cGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCTGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCa                                                       >  4:191144/1‑126 (MQ=255)                                                                               cGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCTGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCa                                                       <  3:194721/126‑1 (MQ=255)                                                                               cGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCTGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCa                                                       <  3:191144/126‑1 (MQ=255)                                                                               cGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCGGCGa                                               >  1:71715/1‑134 (MQ=255)                                                                               cGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCGGCGAAc                                             <  8:268574/136‑1 (MQ=255)                                                                                                tCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCGGCGAACATCCAGCCCAGggcgg                             <  2:163185/135‑1 (MQ=255)                                                                                                       gcCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGAGACCGCCCCACTCGGCGAACATCCAGCCCAGGgcggcggcggcg                     >  7:61077/1‑136 (MQ=255)                                                                                                         cAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGCCCCACTCGGCGAACATCCAGCCCAGGGcggcggcggcgg                    >  3:245988/1‑135 (MQ=255)                                                                                                         cAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACATGGCTTACCGTGACCGACCCACTCGGCGAACATCCAGCCCAGGgcggcggcggcg                     >  4:261587/1‑134 (MQ=255)                                                                                                                   gACGGTGCCCGCAGCCGGGGTGTTGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGTGCACTTGGCTTACCGTGACCGACCCACTCGGCGAACATCCAGCCCAGGGCGGCGGCGGCGGTAGCGATCTGa         <  7:41643/136‑1 (MQ=255)                                                                                                                    aCGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCGGCGAACATCCAGCCCAGGGCGGCGGCGGCGGTAGCGATCTGa         <  2:71715/135‑1 (MQ=255)                                                                                                                      ggTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCGGCGAACATCCAGCCCAGGGCGGCGGCGGCGGTAGCGATCTGAg        >  5:152614/1‑134 (MQ=255)                                                                                                                        tGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCGGCGAACATCCAGCCCAGGGCGGCGGCGGCGGTAGCGATCTGAGTGAc    >  8:100768/1‑136 (MQ=255)                                                                                                                          cccGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCAGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCGGCGAACATCCAGCCCAGGGCGGCGGCGGCGGTAGCGATCTGAGTGACCa  <  7:76549/136‑1 (MQ=255)                                                                                                                                   |                                                                                                                               CACGCCGAAGGCATCCAGCGAATCGTCATAGCCGAGTTTGCGCTTGAGCGTGGTGGCGCAGAAGTAGCAGAGCACACCCGAAACCAGGCCGATGATCAGCGCGCCAATCGGCCCGACGGTGCCCGCAGCCGGGGTGATGGCGACCAGGCCAGCAACCACGCCCGAAGCGATGCCCAGAGCACTTGGCTTACCGTGACCGACCCACTCGGCGAACATCCAGCCCAGGGCGGCGGCGGCGGTAGCGATCTGAGTGACCA  >  KT2440/5960460‑5960716

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: ‑