Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	4,245,154	T→G	40.1%	D175A (GAT→GCT)	CJ019_03982 ←	Crotonyl‑CoA hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	4,245,154	0	T	G	40.1%	32.1 / 8.6	56	D175A (GAT→GCT)	CJ019_03982	Crotonyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base T (24/6);  new base G (2/22);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 9.32e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGC  >  KT2440/4245032‑4245272                                                                                                                           |                                                                                                                       gCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCCGCCCGGGCGCGAACCGCCCCCcacc                                                                                                            >  2:196564/1‑134 (MQ=255)     ggTCTGCGCATAGGCGAAGGTCGGGCCATTTGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCCCGCACGGCGCCCcccaccac                                                                                                        >  8:82362/1‑131 (MQ=255)     ggTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACtc                                                                                                        >  8:263012/1‑135 (MQ=255)     ggTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCAAcgcg                                                                                                       >  4:89894/1‑136 (MQ=255)     ggTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACCACtgc                                                                                                       >  3:272385/1‑134 (MQ=255)         tGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGGGCGCTCCGCGGCCCCCCACacgtcg                                                                                                    >  7:176086/1‑135 (MQ=255)           cgcATCGGCGATGGTCGGGCCATTGGCCAGGCGCCGGGCCAGTGCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAg                                                                                                  <  8:176086/135‑1 (MQ=255)            gcATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACGTCGCTGACCAACCCGGCGGCCACGGCGCGATCGGCGCCCCCCCACTCGTCGAg                                                                                                  >  2:244592/1‑134 (MQ=255)                cggCGACGGGCGGGCCATGGGCCCGGCGCCCGGCCAGTGCGGGGGCAACGGCAGGAAACTGCTCGTCGGCGCAGACTTCGCTCACCACCCCGGGGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAAca                                                                                             <  6:153150/134‑1 (MQ=255)                cggCGACGGGCGGGCCATGGGCCCGGCGCCCGGCCAGTGCGGCGGCACCGGCCGGCAACTGCTCGGCGGCGCAGACTTCGCTGACCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAAca                                                                                             <  2:80131/134‑1 (MQ=255)                   cGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACGTCGCTGACCAACCCGGCGGCCCCGGCGCGCTCCGCGCCCCCCCACACCTGGCGGAAaagca                                                                                          >  6:304105/1‑135 (MQ=255)                    gAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACCACTCGTCGAGGCAccgccg                                                                                         >  2:130705/1‑133 (MQ=255)                    gAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCCGGCGCCGTCGGGGCCCCACCACACCGCGAGGCAcaccacg                                                                                        >  3:327042/1‑134 (MQ=255)                     atggTGGGGCCTTTGGCCAGGCGCCCGGCCTGGCGGGGGGCAAGGGCAGTAACCTCGTCGGCGGCTCAGACTGGGCTGACAACCCCGGCGCCCCGGGCGGTATCGCCGCCCCACAACTCGTCGAGGAACAGCAgg                                                                                        <  6:307437/133‑1 (MQ=255)                     atggTCGGGCCATTGGCCAGGCGCCCGGCAAGTTGGGCGGCAACGGCAGGCAACTGCTGGTCGGGGCTGACTTCGCTCACCAACCCGGCGGCCAGGGCGCGATGGGCGCCCCACAACTCGTCGAGGAACAGCAgg                                                                                        <  8:267626/133‑1 (MQ=255)                        gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGCCCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACCACTCGTCGAGGCACAGCAAgcgc                                                                                     >  6:339609/1‑135 (MQ=255)                        gTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACGCGTCGCGGCACAGCACgcgcg                                                                                    >  3:165537/1‑135 (MQ=255)                               cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACTCGGCGAGGCACAGCAGGCGCGTGGCCt                                                                              >  6:97286/1‑135 (MQ=255)                               cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCACAACTCGTCGAGGGACAGCAGGCGCTTGGCCt                                                                              >  6:196809/1‑135 (MQ=255)                               cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGCGCGGCGCCCCCCAACACGGCGCGGCACACCAGGCGCGTgggccg                                                                             >  4:337145/1‑132 (MQ=255)                               cATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCACGGCGCCCTAGGCGCCCCCCCACACCTCGCGGAACAGCCGGCGCTTGGccg                                                                              >  5:176027/1‑134 (MQ=255)                                   gcccgggCGCCCGGCCAGTGCGGCGCAAAGGGCAGGAATCTGCGGGGGGGCGACCACTGCGCTGCCCACCCCGCCGGCCTGGGCGCTATCGCCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCg                                                                         <  4:349966/131‑1 (MQ=255)                                    gCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGGCGAGGCACAGCAGGCGCTTGGCCTGCTCg                                                                         >  8:213734/1‑135 (MQ=255)                                     ccAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACCACTCGTCGAGGCACAGCAGGCGCTTGGGCTGCTCGc                                                                        >  8:103932/1‑135 (MQ=255)                                        ggCGCCGGGCCAGGGCGGCGGCAAGGGCAGGCAACGGGGCGGCGGCGCCGACTGCGCTGCCCAACCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                    <  4:227270/136‑1 (MQ=255)                                        ggCGCCCGGCCTGTTCGGCGGCACCGCCAGCCAACTGCGCGGCGGCGCAGACTGCGCTCCCCACCCCGGCGCCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGc                                                                     <  2:70380/135‑1 (MQ=255)                                        ggCGCCCGGCCTGTTCGGCGGCAACGGCAGGCAACTGCTCGGCGGCGCAGACTGCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                    <  3:84815/136‑1 (MQ=255)                                         gcgcCCGGCCTGTGCGGCGCCACCGCCCGCCAAGGGCGCGGCGGCCCAGACTGCGCTCCCCACCCGGGCGGCCGGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGcc                                                                    <  2:340185/135‑1 (MQ=255)                                          cgcCCGGACAGGGCGGCGGCCACGGCAGGCAACTGCGCGGCGGCGCAGACTTCGCTGACCACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCg                                                                   <  8:256997/135‑1 (MQ=255)                                            cccGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACCGCCGGCGCTTGGCCTGCTCGCTGCCGAt                                                                 >  4:111092/1‑135 (MQ=255)                                                 ccAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGcc                                                            >  2:33987/1‑135 (MQ=255)                                                   tgtgcggcggCAACGGCAGGCAACGGCGCGGCGGCGCAGACTGCGCTGACCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGt                                                          <  8:365325/131‑1 (MQ=255)                                                        ggcggcAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCGCGGCGCCCCCCCACTCGTCGCGGCACAGCAGGCGCGTGGCCTGCGCGCTGCCGGTCAGCCGGGGccc                                                     >  1:66326/1‑133 (MQ=255)                                                                   caggcaACTGCGCGGCGGCGCAGACTTCGCGGACCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGgc                                         <  4:165537/136‑1 (MQ=255)                                                                    tgccacctgctgGTCGGCCCTGAGTTCGCTGCCCAACCCGGGGGCCAGGGGGGGTGGGCCGCCCCACAACTCGTCGAGGAACAGCTGGGGGTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGgc                                         <  8:55358/129‑1 (MQ=255)                                                                       caACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCACAACTCGGCGCGGAACAGCCCGCGCTTGGCCTGCTCGCTGCCGCTCAGCCGGGGCAGGTGCCCGCTGGCGCCg                                     >  4:74247/1‑136 (MQ=255)                                                                          cTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACGCGGCGCGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGCTCCGCCGTGGCAGGTGCCAGCTGGCGCCGGc                                   >  8:255579/1‑135 (MQ=255)                                                                          cTGCGCGGCGGCGCAGACTTCGCTGACCACCCCGGGGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                  <  4:175497/136‑1 (MQ=255)                                                                           tGCGCGGCGGCGCAGACTGCGCGGACCACCCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCa                                  <  8:20587/135‑1 (MQ=255)                                                                              tcgtcgGCGCAGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGATCGGCGCCCCCCAACTCGTCGAGGCACAGCAGGCGCTTGGCCTGCTCGCTGGCGATCCGCCGTGGCAGGTGCCCGCTGGCGCCGGCATCg                               >  6:364047/1‑135 (MQ=255)                                                                              tcggcgGCGCAGACTGCGCTGACCACCCCGGCGCCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCAt                                 <  5:196809/133‑1 (MQ=255)                                                                               cgtcgGCGCAGACTGCGCTGACCACCCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCggg                             <  3:89894/136‑1 (MQ=255)                                                                                ttcgGCGCATACATGGCTCCCAACCCCGGCGGCCCGGGCGCGAGCGGCGCACCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                            <  3:161565/135‑1 (MQ=255)                                                                                 gcgGCGCTCACTTCGCGGCCCAACCGGGCGGCCTGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                            <  1:256855/134‑1 (MQ=255)                                                                                 gcgGCGCAGACTGCGCTGACCACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                            <  7:313704/134‑1 (MQ=255)                                                                                 gcgGCGCAGACCTCGCTCCCACCCCCGGCGCCCCGGGCGCGAGCGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCgggg                            <  1:289713/134‑1 (MQ=255)                                                                                  cgGCGCTGAGTTCGCGGCCCAACCCGGCGGCCCGGGCGCGAGGGGCCCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCt                          <  4:327042/136‑1 (MQ=255)                                                                                        aGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACTCGGCGAGGGACAGCAGGCGCGTGGCCTGCTCGCTGCCGCTCAGCCGTGGCAGGTGCCCGCTGGCGCCGGCCTCGGGGCGGTAGcc                    >  3:98900/1‑136 (MQ=255)                                                                                        aGACTTCGCTGACCAACCCGGCGGCCCGGGCGCGCTAGGCGGCCCCCCACTCGGCGAGGAACACCCCGCGCGTTGCCAGCTCGCTGCCGCTCAGCCGGGGCAGGTGGCACCTGGCGCCCGCCTCGGGGCGGTatc                     >  6:291401/1‑133 (MQ=255)                                                                                          aCTTCGCTGACCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCa                  <  4:98900/136‑1 (MQ=255)                                                                                               gCTGACCAACCCGGCGGCCCGGGCGCGCTCGGCGCCCCCCCACTCGGCGCGGCACCGCAGGCCGGTTGCCGGGTCCCAGGCGCTCAGCCGGGGGAGGTGGCAGCGGGCGCCGGCCACGGGGCGGTAGccccagcg              >  6:371144/1‑130 (MQ=255)                                                                                                   accaccCCGGCGGCCCGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt         <  4:47132/136‑1 (MQ=255)                                                                                                   aCCAACCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt         <  4:187436/136‑1 (MQ=255)                                                                                                    ccACCCGGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt         <  6:176027/135‑1 (MQ=255)                                                                                                    ccACCCCGGCGGCCAGGGCGCGAGCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGgtgt         <  8:186297/135‑1 (MQ=255)                                                                                                           ggcggcCAGGGCGCGATCGGCGCCCCCCAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCCGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCCGGGGGGTGTAGCCGGc  >  7:34825/1‑135 (MQ=255)                                                                                                                           |                                                                                                                       GCTTGGTCTGCGCATAGGCGAAGGTCGGGCCATTGGCCAGGCGCCCGGCCAGTTCGGCGGCAACGGCAGGCAACTGCTCGTCGGCGCAGACTTCGCTGACCAACCCGGCGGCCAGGGCGCGATCGGCGCCCCACAACTCGTCGAGGAACAGCAGGCGCTTGGCCTGCTCGCTGCCGATCAGCCGTGGCAGGTGCCAGCTGGCGCCGGCATCGGGGCTGTAGCCCATGCTGGTGTAGCCGGC  >  KT2440/4245032‑4245272

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