Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	5,486,639	A→C	26.8%	A396A (GCT→GCG)	cysG_2 ←	Siroheme synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	5,486,639	0	A	C	26.8%	59.9 / 3.0	50	A396A (GCT→GCG)	cysG_2	Siroheme synthase
Reads supporting (aligned to +/- strand):  ref base A (26/9);  new base C (7/7);  total (33/17)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCCCAGCGGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCG  >  KT2440/5486516‑5486774                                                                                                                            |                                                                                                                                        aTGCCCCAGCGGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCCGCGGCCATGGc                                                                                                                              >  2:380884/1‑135 (MQ=255)    ccccAGCGGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCCGCGGCCCTGGCGaa                                                                                                                           >  7:175262/1‑135 (MQ=255)    ccccAGCGGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGa                                                                                                                            >  1:259779/1‑134 (MQ=255)    ccccAGCGGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCGTCGCGCACCGCAGCGGCCATGGCGa                                                                                                                            >  5:364226/1‑134 (MQ=255)     cccAGCGGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCCCCGCAGCGGCCATGGCGaaa                                                                                                                          >  6:214857/1‑135 (MQ=255)       cAGCGGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCAACGCAGCGGCCATGGCGAAAAc                                                                                                                        >  4:444452/1‑135 (MQ=255)         gCGGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGCGCACCGCAGCGGCCATGGCGAAAACGc                                                                                                                      >  1:82451/1‑135 (MQ=255)           ggTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGAGTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGccc                                                                                                                    >  5:343505/1‑135 (MQ=255)            gTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCg                                                                                                                   >  6:423657/1‑135 (MQ=255)                  ctgcctgGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCGGGATTCGGGGCATGGGTCGCGCCCCCCCGCGGGCCTGGCGCAAACGCCCcgggggc                                                                                                             >  5:276047/1‑131 (MQ=255)                       tgGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTgcc                                                                                                        >  1:321385/1‑133 (MQ=255)                          gcggcgGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCCGCGGCCATGGCGAAAACGCCCGGGTCGCCTgacgac                                                                                                     >  4:40003/1‑135 (MQ=255)                          gcggcgGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCGGGGTTCGTGCAATGCTTCGAGCCCCGCAGCGGCCCTGGCGACACCGCCCGGGTCGCCGgacgcc                                                                                                     >  4:334120/1‑133 (MQ=255)                           cggcggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCAAGGCGAAAACGCCCCGGTCGCCTGACGacac                                                                                                   >  7:127467/1‑136 (MQ=255)                            ggcggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGCTTCGAGCAATTGGTCGAGCACCGCCTCGGCCATGGCGAAACCCCCCGGGGCGCCGGACGaaac                                                                                                   >  6:33732/1‑135 (MQ=255)                            ggcggTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCCCCGCCGCGGCCATGGCGAAAACGCCCGGGTCGCCTgacgcccc                                                                                                   >  1:256021/1‑131 (MQ=255)                                  tGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAAAGCTTCGAGCCCCCCCCCGGCCCTGGCGAAAACGCCCGGGTCGCCGGACGCCCcccccgc                                                                                             >  6:400214/1‑135 (MQ=255)                                      cgcgcGGCGGAACCCCCCGGCAATAGCTGCCGGGCAACGCGAGGCCCGGGGGGGGCGCGGGATTCGTGCAAGGCGGCGCCCACCCCCGCGGCCAGGGCAAAAACGCCCGGGTCGCCTGACGACAccaccgccacc                                                                                         <  1:285840/134‑1 (MQ=255)                                      aGCGAGGCGTAAACCCCCGCAATCATCTGCCGGTAACGGCGAGGCCAGGGGGGGGCGGTGTATTCGGCAAATGCTTCGCGCCCCGCCGCGGCCATGGCTAAAACGCCCGGGTCGCCTGACGACAccaccgccacc                                                                                         <  4:242590/135‑1 (MQ=255)                                                   cccccGCCAATAGCTGCTGGTCAACGCGAGGCCTGTGGGGTTGGGTGGTTTGGTGAAATGCTGGGACCACCGCAGCGGCCATGGCTAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                           <  8:37194/136‑1 (MQ=255)                                                    ccccGGCAATAGCGGCCGGCCACCGCGAGCCGGGGGGGGGGCGCGGGATGCGTCCAATGCTGCCACCCCCGCCGCGGCCTTGGCTAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                           <  2:346760/135‑1 (MQ=255)                                                    ccccGGCAAGATCTGCTGGTCAACGCGATGCCAGGTGGGGTCGCTGTATGCGTGCAATGCTGCGCGCACCGCAGGGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGgc                                                                           <  2:350136/135‑1 (MQ=255)                                                          cAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGCCGACACCACCGCCACCGAACGGCCTTGGGCCGCCAg                                                                     >  3:210678/1‑135 (MQ=255)                                                            aGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGtt                                                                   >  1:10173/1‑135 (MQ=255)                                                                 tGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGCAgg                                                             >  8:14735/1‑136 (MQ=255)                                                                           cgcgATGCCAGGTGGGGACGGGGGATTCGTGCAATGGTACGAGCACCCCACCGGCCTTGCCGAAAACGCCCGGGTCGCCGGACTCCACCCCCGCCCCCGACCGGCGTTGGGCCGCCAGTGGGAAGGCATGGCggc                                                    >  4:354774/1‑134 (MQ=255)                                                                            gctaTGCCAGGGGGGGTCGCTGGATTCGTGCAATGCTTCTCGCACCGCCGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgc                                                   <  6:171006/132‑1 (MQ=255)                                                                            gcgATGCCAGGTGGGGTCGCTGGTTTCGTGCAATGCTTCGAGCCCCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGCCGACACCCCCGCCCCCCAACGGCCGTGGGCCGCCAGTTCGCAGGGCTGGCGGgc                                                   >  1:347385/1‑135 (MQ=255)                                                                              gATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGACCGGCCTTGGGCCGCCAGTTCGAAGGCATGTCCGgcgcg                                                >  8:197304/1‑136 (MQ=255)                                                                              gATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGgcgc                                                 >  3:366938/1‑135 (MQ=255)                                                                                  ccAGGGGGGGGCGCGGGAGGCGTGCAACGCTCCGCCCACCCCCCCGCCCAGGGCCAACACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTg                                             <  5:197017/135‑1 (MQ=255)                                                                                   cAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACCCCACCGCCCCCCAACCGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTg                                             >  1:125172/1‑134 (MQ=255)                                                                                           ggcGCTGGATGCGTGCAATGCTTCGAGCCCCGCCGCGGCCTTGGCTAACACGCCCGGGGCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTc                                    <  4:74194/133‑1 (MQ=255)                                                                                              gCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCgg                                >  8:88205/1‑136 (MQ=255)                                                                                               ttGGATTCGGGAATTGCTGCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGGCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGt                               <  8:127467/135‑1 (MQ=255)                                                                                               cTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGGATGTCGGGCGCGCTGCATCTCTTCGCGGt                               >  8:411437/1‑136 (MQ=255)                                                                                                tGGATTCGGGCAAGGCTGCGAGCACCGCGGCGGCCATGGCGAACACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGt                               <  6:189450/135‑1 (MQ=255)                                                                                                gggATTCGTGCAAGGCTGCGCGCCCCGCAGCGGCCATGGCGAACACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGt                               <  3:414308/134‑1 (MQ=255)                                                                                                gggATTCGTGAAATGCTTCGCGCCCCGCCGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGt                               <  1:22705/134‑1 (MQ=255)                                                                                                  gATTAGTTCAATGCTTCGAGCACCTCAGCGGCCATGGCGAAAAAGCCCGGGTCGCCTGCCGACACCACCGCCACCGAACGGCCTTGGGCCGACCGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGGTg                             >  4:342350/1‑135 (MQ=255)                                                                                                    ttCGTGCAATGCGTCGAGCACCGCAGCGGCCATGGCGAAAACGCTCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                          <  7:286024/136‑1 (MQ=255)                                                                                                     tCGTGCAATGCTCCGAGCACCGCAGCGGCCATGGCTAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                          <  1:22697/135‑1 (MQ=255)                                                                                                      cGTGCAATGGTTGGAGCAGCGCAGCGGCCATGGCGAAAACGCCCGTGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCg                          <  3:444452/134‑1 (MQ=255)                                                                                                         gaaaggCTTCTAGAACCGCAGCGGCAATGCGTAAAACGCCCTGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgca                     <  8:63654/131‑1 (MQ=255)                                                                                                           aaTGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGTCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGgtgcag                    >  1:204986/1‑135 (MQ=255)                                                                                                               cTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCCCCGAACTGCCTTGGGCAGCCAGTTCGAAGGAATGTCGGGCCGGCTGCATCTCTTAGCCGTTGTCGgtgcagtgca                >  5:51008/1‑135 (MQ=255)                                                                                                                     gCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCCCCGACCGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTg         >  7:296532/1‑136 (MQ=255)                                                                                                                     gCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCt          >  2:271427/1‑135 (MQ=255)                                                                                                                     gCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCt          >  5:431628/1‑135 (MQ=255)                                                                                                                          gCAGCGGCCATGGCGAAAACGCCCGGGACGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcg     <  4:69938/135‑1 (MQ=255)                                                                                                                            aGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGt                               <  3:156241/107‑1 (MQ=255)                                                                                                                            aGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGt                               >  4:156241/1‑107 (MQ=255)                                                                                                                            aGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgccg  >  7:87999/1‑136 (MQ=255)                                                                                                                            aGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTcgcc   >  2:354122/1‑135 (MQ=255)                                                                                                                            |                                                                                                                                        ATGCCCCAGCGGTGCACCTGCCTGGGCGGCGGTTGCCAGCGAGGCGGATACCCCCGGCAAGATCTGCAGGTCAACGCGATGCCAGGTGGGGTCGCTGGATTCGTGCAATGCTTCGAGCACCGCAGCGGCCATGGCGAAAACGCCCGGGTCGCCTGACGACACCACCGCCACCGAACGGCCTTGGGCCGCCAGTTCGAAGGCATGTCGGGCGCGCTGCATCTCTTCGCGGTTGTCGGTGCAGTGCAGCACCTGGTCGCCG  >  KT2440/5486516‑5486774

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