Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,186,737	T→A	28.0%	I378F (ATC→TTC)	U712_16400 ←	3‑ketoacyl‑CoA thiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,186,737	0	T	A	28.0%	51.6 / 6.1	45	I378F (ATC→TTC)	U712_16400	3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (23/9);  new base A (2/11);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 8.06e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCGATTGCGCCCCCGTTGACATTCACTTTTT  >  CP006881/3186599‑3186868                                                                                                                                                |                                                                                                                                    ggTGGAAACAATTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGAt                                                                                                                                     <  8:275362/139‑1 (MQ=255)     tataCATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGC‑TTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATAcaca                                                                                                                                >  3:279916/1‑139 (MQ=255)          aTTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTACCAATCCGCCACCGATACACCTGGt                                                                                                                            >  6:15339/1‑139 (MQ=255)          aTTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCGTTTTTTTCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAAAA‑CA‑CCGGCCGCGCCCATTCCGCCGCCGATACACATgg                                                                                                                             >  8:36725/1‑138 (MQ=255)            ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTACCATTCCGCCGCCGAAACACAAGGTCa                                                                                                                          >  3:162331/1‑139 (MQ=255)              gATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCCGCCGCGTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCGCCCATTCCGCCGCCGCTCCACCTGGTCACg                                                                                                                        >  5:121585/1‑139 (MQ=255)               aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTGCCCC‑TTTAGC‑ATAATTCAAAAA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGa                                                                                                                       >  5:375964/1‑139 (MQ=255)               aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGa                                                                                                                       >  7:21077/1‑139 (MQ=255)               aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGa                                                                                                                       >  1:268165/1‑139 (MQ=255)               aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTAACgc                                                                                                                       >  4:49659/1‑138 (MQ=255)               aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCAATCCGCCGCCGATACACATGGTCActa                                                                                                                       >  6:117249/1‑137 (MQ=255)               aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTACCATTCCGCCGCCGATACAAATTGTAACgg                                                                                                                       >  5:60843/1‑138 (MQ=255)                 cGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTCCAA‑TTTTTTTTTCCCCA‑TTTAGA‑TTATTTAAAAAA‑CA‑CCGCCCGATCCAATGCCCCCGCCGAACCACAGGTGCACGaga                                                                                                                     >  7:255707/1‑137 (MQ=255)                           aTCCGCGATTAAAAATCCGCCGCGTTTTTGAACGTCAGCCGCTTTTTTTGCCA‑TTTTTTTTCCCCCC‑TTTACC‑AAAATTAAAAAA‑CC‑CCGGCCGCTCCAATTCCGCCGCCGAAACACATGGTCACGACTCCAAATTGCt                                                                                                           <  8:21077/139‑1 (MQ=255)                             ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCCTTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCg                                                                                                         >  2:60443/1‑139 (MQ=255)                             ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTTCtag                                                                                                         >  7:203903/1‑137 (MQ=255)                              cTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTACGCCCCCCATACAAATTGTCACCACACCAAAATgatagt                                                                                                        >  4:92927/1‑134 (MQ=255)                               tCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCGACCATTCCGCCGCCGATACACCTGGTCACGACGACAAATTTCTCGtt                                                                                                       >  7:12569/1‑139 (MQ=255)                                 gATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCA‑TTTAGC‑AAAATTAAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCTATACACATGGGCACGACTCCAAATTGCTCGTTtc                                                                                                     <  2:47286/139‑1 (MQ=255)                                         aTCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTACCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCGGCCAATAAGCCCCCCATACACAAGGTGAAGACTCCAAATTTCTAGTTTGTCCCtttc                                                                                             >  1:184467/1‑139 (MQ=255)                                          cccgccgCCTTTTGGTCCGTCCCCCGCTTTTTTTGCCA‑TTTGTTTTCCCCCC‑TTAAGC‑AAAATTCAAACA‑CC‑CCGGCCGCTCCCATTCCGCGGCCTAAACACTGGGCCACGACTCCAAATTGCTCGTTTCTCCGtttca                                                                                            <  4:192010/138‑1 (MQ=255)                                              ccggatttttGTCCGGCCGCCGATTTTTTTGCCT‑GTTGTTTTACCCCC‑TTTAGC‑ATAATAAAAAAA‑CC‑CCGCCCCCACCCATTCCGCCGCCGTAACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGtttcatttc                                                                                        <  2:282277/134‑1 (MQ=255)                                                  cTTTTTGTAGGCCAGCCTGTTTTTTTGCCT‑TTTGTTTTACCCCA‑TATAG‑‑AAAATTAAAAAT‑CACCCGGCCGCCCCCTAGCCGCCGCCGAAACACTGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATgg                                                                                    <  6:228861/139‑1 (MQ=255)                                                     tttGTCCGCCAGCCGCTTTTTTTGCCA‑TTTGTTTTACCCCC‑TTTACA‑ATAATTCAAAAC‑CC‑CCGGCCGCACCCATACCCCCGCCGATACACATGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAt                                                                                  <  5:380650/138‑1 (MQ=255)                                                     tttGTACGTCCGCCGCTTTTTTTGCCA‑TTTGTTTTCCCCCA‑TTAAGA‑AAAATAAAAATA‑CC‑CCGGCCCCGCCCATCCCGCCGCCTAAACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAt                                                                                  <  3:49659/138‑1 (MQ=255)                                                         tACGTCCGCCGCTTTTTTTGCCA‑TTTGTTTTCCCCCC‑TTTACC‑ATAATTAAAAAC‑CC‑CCGGCCGACCCCTTGCCGCCGCCGAACCACTTGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                                                             <  8:12569/139‑1 (MQ=255)                                                          aCGTCCGCCGTTTTTTTTG‑CA‑TTTTTTTTTCCCCCATTAAGA‑AAAATTAAAATA‑CA‑CCGGCCGCCCCAATGCCGCCGCCTATACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                                                             <  5:15339/138‑1 (MQ=255)                                                           cctCAGCCGC‑TTTTTTGCCATTTTGTTTTCCCCCC‑TTTAGC‑ATAATTCAAAAC‑CC‑CGGGCCGCGCCCATGCCGCCGCCGAAACACAGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                                                            <  2:284865/136‑1 (MQ=255)                                                                  cctgttttttTCCA‑TTTTGTTTACCCCC‑ATTAAC‑AAAATTAAAAAACCC‑CCTGCCGCACCATTCCCGCCGCTGAAACACTTGTTCACGACTCCAAATTGCTCGGTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTa                                                                     <  8:331607/135‑1 (MQ=255)                                                                       tttttGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCtt                                                                <  4:143428/138‑1 (MQ=255)                                                                          ttGCAA‑TTTGTTTTTCCCCA‑TTACGC‑AAAATACAATAA‑CC‑CCGGCCGCTCCCATTCCGCCGCCGATCCACTTGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTtgt                                                             <  3:407227/138‑1 (MQ=255)                                                                               a‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGTATGAGTGAGAGTGTTAGCGTTGTGcgcg                                                        >  1:279793/1‑138 (MQ=255)                                                                                   tGTTTTCCCCCC‑TTTAGC‑ATAATTCAAAAA‑CA‑CCGGCCGCCCCCATTCCGCCGCCGAAACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                                                     <  5:299837/138‑1 (MQ=255)                                                                                   tGTTTTACCCCC‑TTTATA‑ATAATTAAAATA‑CC‑CCGCCCGCTCCCTTCCGGCCGCCGAAACCCATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                                                     <  4:162331/138‑1 (MQ=255)                                                                                   tGTTTTACCCAC‑TTAACC‑ATAATACAATAC‑CC‑CCGGCCGCACCCATTCCGCCTCCGATAAACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                                                     <  8:203903/138‑1 (MQ=255)                                                                                    ttttttCCCCC‑CTTAACAATATAACAAACA‑CC‑CCGCCCGCGCCAATCCCCCCGCCGAAACACAGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGGGTAAGCTTTGTGCCCGTAc                                                     <  1:171732/137‑1 (MQ=255)                                                                                            ccc‑TTTCGC‑ATAATTAAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACAAATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCgg                                           <  8:255707/139‑1 (MQ=255)                                                                                            ccc‑TTAACC‑AAAATTCAAAAA‑CA‑CCGGCCGCCCCCATACCGCCCCCTATACACAGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCgg                                           <  1:60443/139‑1 (MQ=255)                                                                                             ccatttACC‑ATAATTAAAATA‑CA‑CCGGCCGCCCCCATTCCGCCGCCTATACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCgg                                           <  7:366168/136‑1 (MQ=255)                                                                                                tttAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGgtg                                        >  7:106987/1‑139 (MQ=255)                                                                                                           ttCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCCCGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGcgcg                              >  6:68199/1‑139 (MQ=255)                                                                                                           ttCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGGTTCATTTCATGGATGAGTGAGAGTGTTAGCTTTGTGCCCGTACCGCCAAGCGGGTGTCCAAgcggg                              >  6:87465/1‑137 (MQ=255)                                                                                                                aTA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTGTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCGATTgc                         >  3:207686/1‑139 (MQ=255)                                                                                                                          gCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCGATTGCGCCCCCGt                 >  7:151298/1‑139 (MQ=255)                                                                                                                                         ccgccgATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCGATTGCGCCCCCGTTGACATTCACttttt  >  2:365719/1‑139 (MQ=255)                                                                                                                                         ccgccgATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCGATTGCGCCCCCGTTGACATTCACttttt  >  4:249706/1‑139 (MQ=255)                                                                                                                                                |                                                                                                                                    GGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCC‑TTTAGC‑ATAATTCAAATA‑CA‑CCGGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCGATTGCGCCCCCGTTGACATTCACTTTTT  >  CP006881/3186599‑3186868

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