Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	REL606	743,238	G→A	8.8%	G309D (GGC→GAC)	sucB →	dihydrolipoamide acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	743,238	0	G	A	8.8%	135.5 / 2.1	45	G309D (GGC→GAC)	sucB	dihydrolipoamide acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (19/22);  new base A (2/2);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTTACCACAACTATTTCGACGTCAGCATGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGATCGTCCGATGGCGG  >  REL606/743095‑743384                                                                                                                                                |                                                                                                                                                   ggTTTACCACAACTATTTCGACGTCAGCATGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCt                                                                                                                                              >  2:520897/1‑150 (MQ=255)  gTTTACCACAACTATTTCGACGTCAGCATGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTg                                                                                                                                             >  2:367667/1‑150 (MQ=255)          caACTATTTCGACGTCAGCATGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGa                                                                                                                                     >  2:456263/1‑150 (MQ=255)               aTTTCGACGTCAGCATGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATc                                                                                                                                <  1:139141/150‑1 (MQ=255)                 ttCGACGTCAGCATGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTg                                                                                                                              >  2:600145/1‑150 (MQ=255)                             aTGGCGTTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGAGGTCGATACCCGCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAAc                                                                                                                  <  1:348667/150‑1 (MQ=255)                              tGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACt                                                                                                                 <  2:97071/150‑1 (MQ=255)                              tGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACt                                                                                                                 <  1:55077/150‑1 (MQ=255)                              tGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACt                                                                                                                 <  2:474798/150‑1 (MQ=255)                                   gTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGCCAGGCAAAGGCCGTGACGGCAAGCTGGCCGTTGAAGATCTGACCGGTGGTAACTtcccc                                                                                                            >  1:11761/1‑150 (MQ=255)                                                   gCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGgtggtg                                                                                            >  2:74874/1‑150 (MQ=255)                                                       ggTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACGTCGAGAAGAAAATCAAAGCGCTGGCAGTCAAAGGCCGTTACGGCCAGCTGACCGTTGAAGATCTGACCGGTGGTAACTCCACCATCCCCACCGGTGGTGTGtt                                                                                        >  1:502270/1‑150 (MQ=255)                                                           acacCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGt                                                                                    <  2:459102/150‑1 (MQ=255)                                                            cacCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGCAAACGAGAAGAAAATCAAAGAGCTGGCAGACAAAGGCCGTGACGGCAAGCTGACCGTTGACGATCTGACCGGTGGTAACTTCACCACCACGACCGGTGGTGCGTTCGGtg                                                                                   >  1:663549/1‑149 (MQ=255)                                                            cacCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGtt                                                                                   >  2:346120/1‑150 (MQ=255)                                                            cacCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGtt                                                                                   >  2:527333/1‑150 (MQ=255)                                                            cacCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGtt                                                                                   >  2:89187/1‑150 (MQ=255)                                                            cacCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGtt                                                                                   >  2:501325/1‑150 (MQ=255)                                                              ccGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTcc                                                                                 >  2:202343/1‑150 (MQ=255)                                                                 gTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTg                                                                              >  2:575181/1‑150 (MQ=255)                                                                  tactaggTGATGTCGGTACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGa                                                                             <  1:519406/144‑1 (MQ=255)                                                                  tACGGCGTGATGCCGATACCCCCGGCATGGCAGATAGCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGATGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCCGTTCCCTGa                                                                             <  1:361759/150‑1 (MQ=255)                                                                    aggCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATg                                                                           <  2:485525/148‑1 (MQ=255)                                                                      gCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATgct                                                                         >  1:484150/1‑148 (MQ=255)                                                                       cGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCt                                                                        <  2:397076/150‑1 (MQ=255)                                                                           aTGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGc                                                                    <  2:377968/150‑1 (MQ=255)                                                                             gTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCg                                                                  >  2:470956/1‑150 (MQ=255)                                                                             gTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCg                                                                  >  2:539727/1‑150 (MQ=255)                                                                             gTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGACAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCg                                                                  >  2:297124/1‑150 (MQ=255)                                                                             gTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGACAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCg                                                                  >  2:464670/1‑150 (MQ=255)                                                                                     ccTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCaa                                                          >  2:545147/1‑150 (MQ=255)                                                                                          gCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACccgc                                                     >  1:475895/1‑150 (MQ=255)                                                                                               gCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCag                                                <  2:121109/150‑1 (MQ=255)                                                                                                            aaaaacTCAAAGAGCTGGCAGTAAAAGGCCGGGACGGCAAGCTGACCGGTGAGGATCGGACCGGTGGTAACGTAACCATCACCAACGGTGGTGTTTTCGGTTCCCTGATGTTTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTgg                                   <  1:338293/148‑1 (MQ=255)                                                                                                              aaaaTCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGt                                 <  2:428327/150‑1 (MQ=255)                                                                                                              aaaaTCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGt                                 <  2:503702/150‑1 (MQ=255)                                                                                                                            ggCAGTCAAAGGCCGTGACGACAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCaa                   <  1:296278/150‑1 (MQ=255)                                                                                                                            ggCAGTCAAAGGCCGTGACGACAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCaa                   <  1:463352/150‑1 (MQ=255)                                                                                                                                gttaactcgcctGACGTCAAGCTCACCTTTTAAGATCTGACCTTTCGCAACCTCACCATCACCGACGGTAGTGTGTTCCGTCCCCTGATGTCTAGGCCTTCCATCAACCCGCCGCAGGGCGTAATTCTGGGTATGCACGCTATCAAAGAt               <  1:414122/139‑1 (MQ=255)                                                                                                                                gTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGGTGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGAt               <  1:201760/150‑1 (MQ=255)                                                                                                                                gTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGAt               <  1:74655/150‑1 (MQ=255)                                                                                                                                gTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGAt               >  2:550414/1‑150 (MQ=255)                                                                                                                                gTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGAt               <  1:366604/150‑1 (MQ=255)                                                                                                                                gTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGAt               <  2:1058/150‑1 (MQ=255)                                                                                                                                gTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGAt               <  2:65970/150‑1 (MQ=255)                                                                                                                                gTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGAt               <  2:477246/150‑1 (MQ=255)                                                                                                                                gTCAAAGGCCGGGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTGCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGAt               <  1:469616/150‑1 (MQ=255)                                                                                                                                             aCGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGATCGTCCGATGGCgg  <  2:119914/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GGTTTACCACAACTATTTCGACGTCAGCATGGCGGTTTCTACGCCGCGCGGCCTGGTGACACCGGTACTGCGTGATGTCGATACCCTCGGCATGGCAGACATCGAGAAGAAAATCAAAGAGCTGGCAGTCAAAGGCCGTGACGGCAAGCTGACCGTTGAAGATCTGACCGGTGGTAACTTCACCATCACCAACGGTGGTGTGTTCGGTTCCCTGATGTCTACGCCGATCATCAACCCGCCGCAGAGCGCAATTCTGGGTATGCACGCTATCAAAGATCGTCCGATGGCGG  >  REL606/743095‑743384

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