Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,186,711	A→C	23.4%	G386G (GGT→GGG)	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,711	0	A	C	23.4%	63.7 / 6.6	45	G386G (GGT→GGG)	U712_16400	3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (24/10);  new base C (2/9);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 3.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACAC‑ATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCGATTGCG  >  CP006881/3186578‑3186846                                                                                                                                        |                                                                                                                                          tGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑Gc                                                                                                                                     <  7:89241/139‑1 (MQ=255) tGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑Gc                                                                                                                                     <  5:71224/139‑1 (MQ=255)    gcGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGc                                                                                                                                  <  6:243242/138‑1 (MQ=255)     cGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGcg                                                                                                                                 >  1:131526/1‑138 (MQ=255)     cGTCAGTGAAATCCTCCGGTGTATACATTAGCTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGGCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCt                                                                                                                                 >  3:134155/1‑139 (MQ=255)     cGTCAGTGAAATCCTCCCGTGAATACATTTGATCGTATGTGACATCCTCGATTAAAAATCAGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGACTTTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCggt                                                                                                                                 >  8:196344/1‑137 (MQ=255)      gTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTACCCCTTT‑AGCATAATTCAAAAA‑CACCG‑GCCGCTc                                                                                                                                >  4:299647/1‑139 (MQ=255)      gTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCGTTTTTTACGTCCGCCGCGGTTTTTGGCGTGTGGTTTTCCCCCGTT‑AGCATAATTCAATTA‑CCCCG‑GGCGCtg                                                                                                                                >  5:184555/1‑138 (MQ=255)      gTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATCAGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCAAAATTAAAATA‑CACCG‑GCCCCTc                                                                                                                                >  5:259772/1‑139 (MQ=255)                tcctccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGGCATTTGTTTTTCCCCCTTT‑ATCATAATTAAAATA‑CCCCG‑GCCGCGGCCACTtcgcc                                                                                                                      >  1:152476/1‑139 (MQ=255)                tcctccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTccgcc                                                                                                                      <  5:287472/139‑1 (MQ=255)                   tccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCAGTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTccgccgcc                                                                                                                   >  2:254211/1‑139 (MQ=255)                     cGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTcagccgccg                                                                                                                  >  7:270490/1‑138 (MQ=255)                            taCATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACAC‑At                                                                                                         >  5:168987/1‑139 (MQ=255)                            taCATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCGTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGAGACAC‑At                                                                                                         >  3:42739/1‑139 (MQ=255)                            taCATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTATGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCACCTTT‑AGCATAATTAAACTA‑CACCC‑GCCCCTCCCATTCCGCCGCCGATacac                                                                                                            >  5:154643/1‑137 (MQ=255)                                 ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGATCCCATTACCCCGCCGATACAC‑ATGGGCa                                                                                                    >  4:54875/1‑139 (MQ=255)                                  tGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACCC‑ATGGTCAc                                                                                                   >  6:133943/1‑139 (MQ=255)                                    aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTTTTTTTCCCCCGTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCCATACCCCCCCGATACAC‑AAGGTCACgc                                                                                                 >  6:217443/1‑138 (MQ=255)                                    aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCt                                                                                                                                 >  2:235279/1‑108 (MQ=255)                                    aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCt                                                                                                                                 <  1:235279/108‑1 (MQ=255)                                    aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACAC‑ATGGTCACGa                                                                                                 >  3:149960/1‑139 (MQ=255)                                    aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACAC‑ATGGTCACGa                                                                                                 >  7:225012/1‑139 (MQ=255)                                    aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCGTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACAC‑ATGGTCACGa                                                                                                 >  4:108712/1‑139 (MQ=255)                                      cGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCCTTCCAC‑ATGGTCACGGCt                                                                                               >  7:162133/1‑139 (MQ=255)                                        tATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTTGCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTGCCATTCCGCCGCCCCTACAC‑AAGGTCCCGACTc                                                                                              >  2:315462/1‑138 (MQ=255)                                                  ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACCC‑CTGGTCACGACTCCAAATTGCTCg                                                                                   >  6:335666/1‑139 (MQ=255)                                                  ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACAC‑ATGGTCACGACTCCAAATTGCTCg                                                                                   >  5:63232/1‑139 (MQ=255)                                                  ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑AACCG‑GCCGCTCCAATTCCGCCGCCGCTACAC‑ATGTTCACGACTCCAAATTGCTCg                                                                                   >  4:271958/1‑139 (MQ=255)                                                   cTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCC‑GCCGCTCCCATTACGCCCCCCATACCC‑ATGGTCACGACTCCAAAATGCTAGt                                                                                  >  4:265927/1‑139 (MQ=255)                                                             aaTCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCCTTCCGCCGCCGATACAC‑ATGGTCACGACTCCAAATTGCTAGTTTCTACGtt                                                                         >  8:15900/1‑138 (MQ=255)                                                                       cTTTTTGTACGCCCGCCGCTTTTTTTGCCTTTTGTTTTACCCCCTTT‑ACAAAAATAAAAAAC‑CCCCC‑CCCGCCCCAATCCCCCCCCCGATACCC‑ATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATgg                                                              <  3:108712/139‑1 (MQ=255)                                                                              tACGTCCGGCGTTTTTTTTGCTTTTTGTTTAGCCCACTTT‑ACAATAATAAAAAAC‑CCCCG‑CCCGCCCCCATCCCGCCCCCTAACCAC‑AGGCGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                                       <  6:34636/139‑1 (MQ=255)                                                                              tACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACAC‑ATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                                       <  6:63232/139‑1 (MQ=255)                                                                               ccgcccgCCGCTTTTTTTCCCATTTG‑TTTTCCCCCCTT‑ACCAAAAAAAAAAAACCCCCG‑CCCGCTCCCATCCCGCCGCCCAAACAC‑TGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                                      <  4:149960/133‑1 (MQ=255)                                                                                cGTCCGCCGCTTTTTTTGCCTTTTGTTTTTCCCCTTTT‑ATCATAATACAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACAC‑ATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                                      <  5:80751/138‑1 (MQ=255)                                                                                cGTCAGCCGCTTTTTTTCCAATTTGTTTTACCCCTTTA‑CGAATAATACAAAAC‑CACCG‑GCCGCCCCCATTCCGCCCCCGAAACAC‑ATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                                      <  8:225012/138‑1 (MQ=255)                                                                                      cggttttttttGACTTTTGTTTTCCCCCATTT‑ACAAAAATAAAAAAC‑CCCCG‑GCCCCACCCATTCCGCCCCAGAAACCC‑TGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTa                                               <  4:31208/135‑1 (MQ=255)                                                                                      ccGCTTTTTTGGCCATTTTTTTTTCCCCATAT‑AGAAAAATTCAAAAA‑CACCG‑GCCCATCCCATTCCGCCGCCTAAACAA‑ATGCTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTa                                               <  7:15900/139‑1 (MQ=255)                                                                                              tttCCCATTT‑TTTTTCCCCCTTTAATAATAAATAAAATC‑CCCCG‑GCCGCTCCATTTCCGCCGCCTATAAAC‑TTGCCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTtg                                        <  5:217443/138‑1 (MQ=255)                                                                                                       tttttttCCCCATTT‑AACAAAAATACAAAA‑CCCCGCGCCCCCCCAATCCCCCCGCCTAAACAC‑ATGCCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTa                                <  1:315462/136‑1 (MQ=255)                                                                                                       ttgtttcCCGCCTTA‑ACCATATATAAAATA‑CCCGG‑GCGGCTCCCATGCCGCCGCGGAAACAG‑AGGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCTTGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                               <  6:259772/135‑1 (MQ=255)                                                                                                       tGTTTTACCCGTTTT‑AGAAAAATTAAAATA‑CCCCG‑GCCGCTCCCTTTGCGCCGCAG‑TAAACGTTGGGCACGACTTCAATTTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTTAGCTTTGTGCGCGttc                               <  6:154643/138‑3 (MQ=255)                                                                                                                cccTAT‑AGAAAAATAAAAATC‑CCCCG‑GCCCGTCCCAATCCGCGCCCTATACAC‑AGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCgg                     <  3:54875/139‑1 (MQ=255)                                                                                                                                    tA‑CACCG‑GCCGCCCCCATTCCGCCGACGATACAC‑AGGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCGATTgcg  <  8:34210/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                          TGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTT‑AGCATAATTCAAATA‑CACCG‑GCCGCTCCCATTCCGCCGCCGATACAC‑ATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCGATTGCG  >  CP006881/3186578‑3186846

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