Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	189,506	T→C	15.2%	T1455T (ACT→ACC)	CJ019_00185 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	189,506	0	T	C	15.2%	60.0 / 6.3	34	T1455T (ACT→ACC)	CJ019_00185	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/9);  new base C (2/3);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

CCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCG  >  KT2440/189380‑189639                                                                                                                               |                                                                                                                                      ccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCCGCGTCACCATCGAGAACGCCACTGGCGGCAAc                                                                                                                              >  4:165998/1‑136 (MQ=25)        aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTgg                                                                                                                                     <  6:66067/122‑1 (MQ=32)        aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTgg                                                                                                                                     >  5:66067/1‑122 (MQ=32)        aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGaa                                                                                                                        >  8:358212/1‑135 (MQ=38)        aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGa                                                                                                                         >  2:303791/1‑134 (MQ=21)         aCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAAc                                                                                                                       <  1:303791/135‑1 (MQ=37)            gCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGCCCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGc                                                                                                                    >  7:176617/1‑135 (MQ=255)             ccAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGCCGTCTACAACAACGGTTCGACCGTCAGCGTCCCCATCGAGAACTCCACTGGCGGCAACTTCGAACAGCt                                                                                                                   >  8:365545/1‑135 (MQ=255)             ccAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCt                                                                                                                   >  7:124876/1‑135 (MQ=255)                aaCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGAcc                                                                                                               <  4:315734/136‑1 (MQ=255)                         cggtGACGCCGGCACGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACTTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGACACTGGGGGCAACTTCGAACAGCTGACCCCGAACCCg                                                                                                      <  4:10973/133‑1 (MQ=11)                          gTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGa                                                                                                     <  4:187811/136‑1 (MQ=37)                          gTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGa                                                                                                     >  4:37704/1‑136 (MQ=37)                          gTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAACTTCGAACAGCTGACCCCGAACCCg                                                                                                      <  8:162861/135‑1 (MQ=17)                              aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGCCCCCGGAGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCACCATAGAACAGCTGCCCCCGAACCCGACGc                                                                                                  >  5:105391/1‑135 (MQ=11)                                gcgggCAAAACCATGGGCAGCGTCGTTTTCCAGACCCCGGCGAATGACGTCTACAACAAGGGTTCGACCGTCAGGGTCAGCATCGAGAACGCCACTGGCGGCAACTTCTAACAGCTGACCCCGAACCCGACGCCg                                                                                                <  8:5349/132‑1 (MQ=11)                                gCCGGCATGACCCAGGGCAGCGTCTTTATCCAGACCCCGGCGAATGACGTCTAGAACAACGGTTCGACGGTCTGCGTCACCTTCGAGAACGCCACTGGCGGCTAGTTCGAACAGCTGACCCCGAACCCGACGCCg                                                                                                <  7:211052/135‑1 (MQ=11)                                gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGCCCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGGAAGTTCGAACAGCTGACCCCGAACCCGACGCCg                                                                                                >  7:104480/1‑135 (MQ=255)                                gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGGAAGTTCGACCACCTGACCCCGAACCCGACGCCg                                                                                                >  6:158727/1‑135 (MQ=11)                                      aaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTggcgg                                                                                                                                  <  1:75757/92‑1 (MQ=18)                                      aaGACCGAGGGCAGCGGCTATGTCCACACCCCGGCGAATGACGGCTACAACATCGGTGCGACCGTCAGCGTCACCATCGAGAACCCCACCGGCGGAAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGa                                                                                         <  3:255122/136‑1 (MQ=11)                                            ctgggCTGCGTGGATTTCCAGACCCCGGCGAATGACGTATACAACAACGGGTCGACCGTCTGCGTCACCATCTAGAGCGCCAGTGGGGTAAATTTCGATCTGCTGACCCCTAACCCGACGCCGGGTCTGACCACg                                                                                    <  8:42662/133‑1 (MQ=2)                                                       cGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGACCCCGACGCCGGCTCAGACCACGCTCAACGACTc                                                                         >  8:54668/1‑135 (MQ=21)                                                       cGATTTCCAGA‑CCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGGAAATTCGAACAGCTGACCCCGAACCCCACCCCCGCACAGACCACGCTCAACGACTCg                                                                        >  8:350238/1‑135 (MQ=11)                                                            tCCAGACCCCGGCGAATGACGTCTACAACAAAGGTTCGCCCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGCCGCCGGCTCAGACCACGATCAACGACTCGGTCg                                                                    >  7:336239/1‑135 (MQ=255)                                                                 aCCCCGGCGAATGACGTCTACCACAACGGTTCGACCGTCAGCGTCACCATCGAGACCGCCACTGGCGGCCACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCCCGATCACCGCCTCGGTCCATGcc                                                               >  2:296775/1‑135 (MQ=11)                                                                           aTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGa                                                     >  7:110425/1‑135 (MQ=32)                                                                                              ttCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGCCCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTc                                  >  7:134006/1‑135 (MQ=17)                                                                                                cGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGCACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGCTGCCACCACCGCGACCCTGCCGGCGAGCCCGTCgg                                >  5:356309/1‑135 (MQ=12)                                                                                                cGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCgg                                >  5:112408/1‑135 (MQ=32)                                                                                                           gTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggc                     >  1:293959/1‑135 (MQ=32)                                                                                                            gcACCAGCGAGAACCCCACTGGCGCAAACTGCGAACCGCGGACCCCGAACCCGCCGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                    <  8:64777/134‑1 (MQ=11)                                                                                                                             aCTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACAcc   >  5:16753/1‑135 (MQ=32)                                                                                                                             aCTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCg  >  3:341938/1‑136 (MQ=32)                                                                                                                               |                                                                                                                                      CCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCG  >  KT2440/189380‑189639

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