Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	4,809,650	T→C	9.1%	Y66H (TAC→CAC)	CJ019_04517 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	4,809,650	0	T	C	9.1%	127.8 / 4.9	44	Y66H (TAC→CAC)	CJ019_04517	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/17);  new base C (2/2);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

ACATCGCTTCTTTGGTCAGCGCTGGCCCCTTGCATACTCTTCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTATTCCTGACTGCCCACACACT  >  KT2440/4809545‑4809777                                                                                                          |                                                                                                                                acaTCGCTTCTTTGGTCAGCGCTGGCCCCTTGCATACTCTTCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTcg                                                                                                               <  4:40057/124‑1 (MQ=255) acaTCGCTTCTTTGGTCAGCGCTGGCCCCTTGCATACTCTTCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTcg                                                                                                               >  3:40057/1‑124 (MQ=255) acaTCGCTTCTTTGGTCAGCGCTGGCCCCTTGCATACTCTTCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCGGGTTCTACTCATACATCTCTGCCACCCTCGCGTTGTTTGGa                                                                                                    >  5:24543/1‑135 (MQ=255)                            ccTTGCATACTCTTCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCt                                                                                                                 >  4:29250/1‑95 (MQ=255)                            ccTTGCATACTCTTCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCt                                                                                                                 <  3:29250/95‑1 (MQ=255)                            ccTTGCATACTCTTCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAggg                                                                        <  8:120955/136‑1 (MQ=255)                                 cATACTCTTCATCATCACCCAGCCAACTCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACGCATACATCTCTGCCACCCTCTCGTTCTTTGGATTAAGCGTCGACCTAATCGGCAGTAGGGAAAg                                                                    >  3:232498/1‑135 (MQ=255)                                      tctTCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCt                                                               >  1:28404/1‑135 (MQ=255)                                        ttCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGcgt                                                              >  1:256698/1‑132 (MQ=255)                                           atcatcaCCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCAcg                                                          >  5:8820/1‑135 (MQ=255)                                           accatcaCCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTTGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCAcg                                                          >  4:240127/3‑135 (MQ=255)                                             catcaCCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCAcgcg                                                        >  2:229837/1‑135 (MQ=255)                                                  cccAGCCAACGCATGCAGCCTCGACCTCTACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   <  6:40648/135‑1 (MQ=255)                                                  cccAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   <  3:119680/135‑1 (MQ=255)                                                  cccAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   <  6:75960/135‑1 (MQ=255)                                                  cccAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   <  1:229837/135‑1 (MQ=255)                                                   ccAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   >  3:40044/1‑134 (MQ=255)                                                     aGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTTGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGGGTCCCGCGGGCCTTCa                                                >  5:28267/1‑135 (MQ=255)                                                                    ccTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAgcgc                                <  7:676/136‑1 (MQ=255)                                                                    ccTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAgcg                                 >  4:227309/1‑135 (MQ=255)                                                                    ccTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCCGTCCTTCATCTACCAACACAgcg                                 >  7:277403/1‑135 (MQ=255)                                                                         aCCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTTGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATc                            <  2:256698/135‑1 (MQ=255)                                                                         aCCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAggg                                                                        >  3:165908/1‑91 (MQ=255)                                                                         aCCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAggg                                                                        <  4:165908/91‑1 (MQ=255)                                                                          ccTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAAt                             >  4:54097/1‑133 (MQ=255)                                                                           cTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCaa                          >  4:122295/1‑135 (MQ=255)                                                                                 cGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   >  4:237232/1‑104 (MQ=255)                                                                                 cGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   <  3:237232/104‑1 (MQ=255)                                                                                 cGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTAtt                    >  4:74810/1‑135 (MQ=255)                                                                                 cGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTAtt                    >  2:257325/1‑135 (MQ=255)                                                                                 cGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTATTc                   <  7:98354/136‑1 (MQ=255)                                                                                    ttCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTATTCCt                 <  2:253953/135‑1 (MQ=255)                                                                                      cAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTATTCCTGa               >  2:37550/1‑135 (MQ=255)                                                                                      cAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTATTCCTGa               >  4:40830/1‑135 (MQ=255)                                                                                      cAGCCGCCTGGTTCTACTCACACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTGATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   <  8:35368/99‑1 (MQ=255)                                                                                      cAGCCGCCTGGTTCTACTCACACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   >  7:35368/1‑99 (MQ=255)                                                                                      cAGCCGCCTGGTTCTACTCACACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   >  7:39782/1‑99 (MQ=255)                                                                                      cAGCCGCCTGGTTCTACTCACACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCt                                                   <  8:39782/99‑1 (MQ=255)                                                                                            ccTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGGCCTTCATCTACCAACACAGCGCAATCAAGCTATTCCTGACTGccc         >  2:34603/1‑135 (MQ=255)                                                                                              tGGTTCTACTCCTACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCTGGAGGAGGGATAGTGGCTTCACGCGGTCCTTCATCTACCAACTCAGCGCAATCAAGCTATTCCTGACTGCCcac       <  5:220827/135‑1 (MQ=255)                                                                                                 ttCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTATTCCTGACTGCCcaca      <  3:122295/133‑1 (MQ=255)                                                                                                    tACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTATTCCTGACTGCCCACACACt  <  4:232086/134‑1 (MQ=255)                                                                                                      cTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAggg                                                                        <  3:996/62‑1 (MQ=255)                                                                                                      cTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAggg                                                                        >  4:996/1‑62 (MQ=255)                                                                                                          |                                                                                                                                ACATCGCTTCTTTGGTCAGCGCTGGCCCCTTGCATACTCTTCATCATCACCCAGCCAACGCATGCAGCCTCGACCTATACCGCTTCAGCCGCCTGGTTCTACTCATACATCTCTGCCACCCTCGCGTTCTTTGGATTCAGCTTCTACCTCATCGGCAGGAGGGAAAGTGGCTTCACGCGGTCCTTCATCTACCAACACAGCGCAATCAAGCTATTCCTGACTGCCCACACACT  >  KT2440/4809545‑4809777

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