Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	KT2440	6,167,267	T→C	12.8%	N403D (AAC→GAC)	glmU ←	Bifunctional protein GlmU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	KT2440	6,167,267	0	T	C	12.8%	76.2 / 5.5	32	N403D (AAC→GAC)	glmU	Bifunctional protein GlmU
Reads supporting (aligned to +/- strand):  ref base T (7/20);  new base C (2/2);  total (9/23)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

CAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACC  >  KT2440/6167137‑6167399                                                                                                                                   |                                                                                                                                     cAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCg                                                                                                                                    <  4:44610/133‑1 (MQ=255)      cGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATg                                                                                                                            <  7:53679/136‑1 (MQ=255)       gAAATGTTGCGCTGGCGCGCACGTGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGTCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGa                                                                                                                           <  7:18183/136‑1 (MQ=255)        aactGTTGCGCTGGCGCGCACGGGCTACAGCCAGGGCACCCGCTTCTACAGCCTGGGTGATGGGCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGa                                                                                                                           <  1:218769/132‑1 (MQ=255)                cgcTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACTGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCt                                                                                                                   <  2:263679/135‑1 (MQ=255)                        cgcACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCAt                                                                                                           <  1:193886/135‑1 (MQ=255)                        cgcACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATc                                                                                                          <  7:214136/136‑1 (MQ=255)                             gggCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCGTGATTTCCCCAGGCGCCCCCAGCGCGTTGTTGGCACCGCTGAACCCGTCTTCCCCCCTCACgg                                                                                                      >  2:193886/1‑135 (MQ=255)                                 tACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCtt                                                                                                  <  2:114547/135‑1 (MQ=255)                                            cACTCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGgc                                                                                      >  7:212186/1‑136 (MQ=255)                                            cACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCCGCTGCGGGGGTCGCGCCGGCCGTGATTTCCACAGGCGCCACCCGCCAGTTGGTAGAACCGCTGACCACGGGTTCACCCATCCCGGGCTTGAACatgtttg                                                                                       >  1:233488/1‑133 (MQ=255)                                                 gCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCa                                                                                  <  2:76596/135‑1 (MQ=255)                                                                  gATGGTCGAGCCGGCTGCGGTGGTAGCGCCGGCCGTGATTTCCACAGGCGGCACCAGCGAGTTTTTCGAACCGATGAACACGTCCTCACCCATCACGGGCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATg                                                                >  7:22166/1‑136 (MQ=255)                                                                    tGGGCGAGCCGGCTGCGGTGGGCGCGCCGGCCTTGATTTCCACCGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGGCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGt                                                              <  7:71270/136‑1 (MQ=255)                                                                      gTCGAGCCGGCTGCGGTGGGCGCGCCGGCCTTGTTTTCCAGAGGCGCCACCAGCGAGTTGTTCGACCCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATGGTAGTTGCAGGTGATGGt                                                              <  5:67846/134‑1 (MQ=255)                                                                        cGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAAACGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGcc                                                           >  4:161498/1‑135 (MQ=255)                                                                         gAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCgg                                                         <  7:79775/136‑1 (MQ=255)                                                                          aGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCgg                                                         <  1:152454/135‑1 (MQ=255)                                                                          aGCCGGCTGCGGTGGGCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCgg                                                         <  2:200901/135‑1 (MQ=255)                                                                            ccGGCAGCGGTGGTAGCCCCCGCCTTGATTTGCACACGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTGATCCCCCATCACCGTCTTGAACATGTTGGCGCCAAAGGAGTTGCAGGTTATGGTGCCGgcc                                                       >  2:38983/1‑134 (MQ=255)                                                                                                 ggcTTGATTTCCACTGGCGCCACCAGGGTGTTGTGGGAACCGATGAACACGTCTTCAGCCTTCCGGGGCTTGAACTTGTTGGGGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCg                                  <  4:234164/133‑1 (MQ=255)                                                                                                      gATTGCCACCGGCGCCACCAGCGAGTTGTTCGACCCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGGGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTc                             <  2:261586/135‑1 (MQ=255)                                                                                                       aTTTCCACCGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGGGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTc                             <  1:166134/134‑1 (MQ=255)                                                                                                       aTTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTc                             <  3:11500/134‑1 (MQ=255)                                                                                                       aTTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCg                            <  8:212186/135‑1 (MQ=255)                                                                                                          tCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCAt                        >  7:277940/1‑136 (MQ=255)                                                                                                              caGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATccc                     <  2:193350/135‑1 (MQ=255)                                                                                                                    ccaccaGCGAGTTGTCCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAg                                                                            <  6:171608/74‑1 (MQ=255)                                                                                                                    ccaccaGCGAGTTGTCCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAg                                                                            <  6:171606/74‑1 (MQ=255)                                                                                                                    ccaccaGCGAGTTGTCCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAg                                                                            >  5:171608/1‑74 (MQ=255)                                                                                                                    ccaccaGCGAGTTGTCCGAACCGATGAACACGTCTTCACCCATCACGGTCATGAACTTGTTGGCGCCATCGTAg                                                                            >  5:171606/1‑74 (MQ=255)                                                                                                                             aggtgttCGAACCGATGAACACGGCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGt      <  3:250090/132‑1 (MQ=255)                                                                                                                             aGTTGTTCGATCCGATGAACACGTCTTCACCGATGACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGAGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGt      <  1:38983/135‑1 (MQ=255)                                                                                                                                  ttCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGATCTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGAcc  <  6:215606/134‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     CAGCCCGAAATGTTGCGCTGGCGCGCACGGGCTACAGCCAGGTCACCCGCTTCTACAGCCTGGGTGATGGTCGAGCCGGCTGCGGTGGTCGCGCCGGCCTTGATTTCCACAGGCGCCACCAGCGAGTTGTTCGAACCGATGAACACGTCTTCACCCATCACGGTCTTGAACTTGTTGGCGCCATCGTAGTTGCAGGTGATGGTGCCGGCGCCGATATTGGTGCGCGCACCGATTTCGGCATCCCCCAGGTAGGTCAGGTGACC  >  KT2440/6167137‑6167399

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