Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,327,353	A→G	6.7%	F5L (TTC→CTC)	mbtH ←	Protein MbtH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,327,353	0	A	G	6.7%	191.7 / 4.0	60	F5L (TTC→CTC)	mbtH	Protein MbtH
Reads supporting (aligned to +/- strand):  ref base A (34/22);  new base G (2/2);  total (36/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 = 7.88e-01

GCCAGGCAGTCCTTCTTCAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGAT  >  ACB143/4327220‑4327484                                                                                                                                      |                                                                                                                                    gCCAGGCAGTCCTTCTTCAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGaaaa                                                                                                                                   <  7:274146/136‑1 (MQ=255)    aGGCAGTCCTTCTTCAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACgg                                                                                                                                >  8:207495/1‑136 (MQ=255)       cAGTCCTTCTTCAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAg                                                                                                                             <  5:223127/136‑1 (MQ=255)        aGTCCTTCTTCAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGt                                                                                                                            >  5:346471/1‑136 (MQ=255)         gTCCTTCTTCAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGt                                                                                                                            >  2:226154/1‑135 (MQ=255)         gTCCTTCTTCAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGt                                                                                                                            >  6:123026/1‑135 (MQ=255)         gTCCTTCTTCAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGATCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGt                                                                                                                            >  2:220488/1‑135 (MQ=255)                  cAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGcc                                                                                                                     <  2:127758/133‑1 (MQ=255)                  cAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGcc                                                                                                                     <  3:298700/133‑1 (MQ=255)                   aGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCTGTCGAAAACGGAAGTCATTGCCTGtc                                                                                                                 <  8:123909/136‑1 (MQ=255)                       cGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCt                                                                                                             >  6:206680/1‑136 (MQ=255)                       cGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGACGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCt                                                                                                             >  6:390353/1‑136 (MQ=255)                         cTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCtt                                                                                                            >  7:332622/1‑135 (MQ=255)                            ttGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAat                                                                                                        >  3:380734/1‑136 (MQ=255)                                       tCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAatg                                                                                                       <  5:23441/126‑1 (MQ=255)                                       tCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAatg                                                                                                       >  6:23441/1‑126 (MQ=255)                                       tCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGc                                                                                             >  5:137225/1‑136 (MQ=255)                                         gccagccaTTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGc                                                                                             <  4:319515/134‑1 (MQ=255)                                         gccagccaTTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGc                                                                                             <  8:332622/134‑1 (MQ=255)                                            agccaTTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTgcgc                                                                                        >  2:251843/1‑136 (MQ=255)                                               caTTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGcc                                                                                     <  3:312991/136‑1 (MQ=255)                                                 ttGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTg                                                                                   <  3:336825/136‑1 (MQ=255)                                                   ggggATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTgg                                                                                  >  4:242844/1‑135 (MQ=255)                                                      gATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCa                                                                              <  5:35357/136‑1 (MQ=255)                                                           cTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTTCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGAGCAgggtg                                                                         >  1:200237/1‑136 (MQ=255)                                                              tGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAgggtgggt                                                                      >  8:115485/1‑136 (MQ=255)                                                                  gTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTc                                                                   <  1:114540/135‑1 (MQ=255)                                                                        ccAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGa                                                            >  2:130411/1‑136 (MQ=255)                                                                        ccAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGATTCAAGTAGa                                                            >  6:200772/1‑136 (MQ=255)                                                                        ccAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGACTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGCCGCAGGGTGGGTTTCAAGTAGa                                                            >  8:319223/1‑136 (MQ=255)                                                                             tCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAg                                                       >  4:210026/1‑136 (MQ=255)                                                                             tCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAg                                                       >  7:278137/1‑136 (MQ=255)                                                                               gAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGAAATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTTCGCGCCTGGCCGAGGGAGGGTTTCAAGTAGACGTAGc                                                      >  4:134428/1‑135 (MQ=255)                                                                                         ctTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGCGTGGGTTTCAAGTAGACGTAGCACGCGCTGTa                                            >  4:54715/1‑135 (MQ=255)                                                                                          ttCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTaa                                           >  3:208556/1‑135 (MQ=255)                                                                                          ttCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCCCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTaa                                           >  6:94991/1‑135 (MQ=255)                                                                                             gttgTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAAt                                       <  8:238903/133‑1 (MQ=255)                                                                                             gTGGTTGAGCACTACCTGGAACTGGATATCGTCGAGGTATAAAAAGGAAGTCATTTCCTGTCTCCTTTAATGATGATGCCGCTGCGGGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAAt                                       >  8:359658/1‑136 (MQ=255)                                                                                             gTGGTTGACCACTACCTGGAACTGGATATCGTCTCGTTACAAAAAGGAAGTCATTTCCTGTGTCCTTTAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTattaat                                       >  8:354409/1‑131 (MQ=255)                                                                                             gTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAAt                                       <  8:139939/136‑1 (MQ=255)                                                                                             gTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTAGAAAAAGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAAt                                       >  4:22365/1‑136 (MQ=255)                                                                                               ggTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCt                                                                                                                    >  8:73948/1‑57 (MQ=255)                                                                                               ggTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCt                                                                                                                    <  7:73948/57‑1 (MQ=255)                                                                                                           ccTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCACAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGAc                          <  3:242844/135‑1 (MQ=255)                                                                                                                aaCTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGGCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTACCGGAATTGACGCGccc                    <  4:69921/136‑1 (MQ=255)                                                                                                                   tGGATATCGTCGCGGTCGAGAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAAt                                       <  7:359719/114‑1 (MQ=255)                                                                                                                   tGGATATCGTCGCGGTCGAGAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAAt                                       <  7:365025/114‑1 (MQ=255)                                                                                                                   tGGATATCGTCGCGGTCGAGAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAAt                                       >  8:365025/1‑114 (MQ=255)                                                                                                                   tGGATATCGTCGCGGTCGAGAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAAt                                       >  8:359719/1‑114 (MQ=255)                                                                                                                   tGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCAt                 >  2:121177/1‑136 (MQ=255)                                                                                                                        atCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGc            >  3:214231/1‑136 (MQ=255)                                                                                                                            tcgCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACg        <  6:137225/136‑1 (MQ=255)                                                                                                                             cgcgGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACgg       >  1:103122/1‑136 (MQ=255)                                                                                                                                 gTCTAAAACGTAAGTCATTGCCTGTCTCCTTAAATGATGCTGCCGCTGCGCGCCTGGGGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGa   <  4:214231/136‑1 (MQ=255)                                                                                                                                 gTCTAAAACGGAAGTCATTCCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGCCGCGCCCCATGTGGCCACGGCCGa   <  4:218196/136‑1 (MQ=255)                                                                                                                                 gTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGTGGAGGCCGCGGCGCGCGTGGCGCAGGGTGGGTTTAACGACGAGTAAGCAAGCGAGGTAAAAATTAACGGAATTCACGCGCCCCATGTGGCCACGGCCGa   <  2:398332/136‑1 (MQ=255)                                                                                                                                 gTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGa   <  7:186433/136‑1 (MQ=255)                                                                                                                                  tCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCg                                                                                              <  8:115283/44‑1 (MQ=255)                                                                                                                                  tCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCg                                                                                              >  7:115283/1‑44 (MQ=255)                                                                                                                                  tCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGAt  >  2:29526/1‑136 (MQ=255)                                                                                                                                  tCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGACACGGCCGAt  >  6:198012/1‑136 (MQ=255)                                                                                                                                      |                                                                                                                                    GCCAGGCAGTCCTTCTTCAGGCCGCTTTTGCCCACGGCTCGCCAGCCATTGGGGATGGCTTTGTAGTCGGGCCAGATCGAGTATTGCTCTTCGTGGTTGACCACTACCTGGAACTGGATATCGTCGCGGTCGAAAACGGAAGTCATTGCCTGTCTCCTTAAATGATGATGCCGCTGCGCGCCTGGCGCAGGGTGGGTTTCAAGTAGACGTAGCAAGCGCTGTAAAAATTAACGGAATTGACGCGCCCCATGTGGCCACGGCCGAT  >  ACB143/4327220‑4327484

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