Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,139,583	G→A	R286H (CGT→CAT)	yghU →	putative S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,139,583	0	G	A	97.3%	120.8 / ‑6.6	37	R286H (CGT→CAT)	yghU	putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (17/19);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTAC  >  NC_000913/3139447‑3139703                                                                                                                                         |                                                                                                                         cGCTGGGCGAAAGGAGTAGGCTAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATc                                                                                                                        <  2:48966/139‑1 (MQ=255) cGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATc                                                                                                                        <  1:755384/139‑1 (MQ=255) cGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATc                                                                                                                        >  2:900095/1‑139 (MQ=255)   cTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTGACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAg                                                                                                                      >  1:908556/1‑139 (MQ=255)      ggCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAgggg                                                                                                                   >  1:925516/1‑139 (MQ=255)       gCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGt                                                                                                                  <  2:475935/139‑1 (MQ=255)       gCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGt                                                                                                                  >  1:65778/1‑139 (MQ=255)       gCGAAAGAAGTAGGCGAACGGCCGGCGGTGAAACGTTGGCGTATTGTTAACCGCACCAACGGACCGCTGAATTAGCAGTTGCATGAGCGCCAGGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGt                                                                                                                  <  1:91044/139‑1 (MQ=255)        cGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGAATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGGGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTa                                                                                                                 >  2:684118/1‑139 (MQ=255)             gaagTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGt                                                                                                            <  1:796194/139‑1 (MQ=255)                    gCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCTAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGtt                                                                                                     <  1:388530/139‑1 (MQ=255)                           tCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCgca                                                                                              >  2:678020/1‑139 (MQ=255)                                                           gCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcg                                                              >  1:918628/1‑139 (MQ=255)                                                           gCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcg                                                              >  1:507888/1‑139 (MQ=255)                                                            cACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgc                                                             >  2:1051805/1‑139 (MQ=255)                                                             aCCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgcc                                                            >  2:818490/1‑139 (MQ=255)                                                                    gACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATc                                                                                    <  1:655210/108‑1 (MQ=255)                                                                    gACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATc                                                                                    >  2:655210/1‑108 (MQ=255)                                                                     aCCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                    <  2:767393/139‑1 (MQ=255)                                                                       cGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                  >  2:874635/1‑139 (MQ=255)                                                                                 gCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGcaca                                        <  2:972142/139‑1 (MQ=255)                                                                                    gTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATg                                     >  1:869891/1‑139 (MQ=255)                                                                                      tGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGcc                                   <  1:874635/139‑1 (MQ=255)                                                                                      tGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGcc                                   <  2:339201/139‑1 (MQ=255)                                                                                                ccATGACGCCAGTGGTTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACg                         >  2:979879/1‑139 (MQ=255)                                                                                                    gACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGc                     <  2:1093013/139‑1 (MQ=255)                                                                                                      cGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCcg                   <  2:4204/139‑1 (MQ=255)                                                                                                          aGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTc               >  1:798824/1‑139 (MQ=255)                                                                                                              aTTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt           <  2:798824/139‑1 (MQ=255)                                                                                                              aTTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCAT‑CAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATc          <  2:822186/139‑1 (MQ=255)                                                                                                                 tCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAg                                                                                         <  2:837710/58‑1 (MQ=255)                                                                                                                 tCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAg                                                                                         >  1:837710/1‑58 (MQ=255)                                                                                                                 tCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAg        <  1:869380/139‑1 (MQ=255)                                                                                                                  cGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAgg       >  1:989879/1‑139 (MQ=255)                                                                                                                   gagaCGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcg                                                              >  2:38621/1‑83 (MQ=255)                                                                                                                   gagaCGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcg                                                              <  1:38621/83‑1 (MQ=255)                                                                                                                       cGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTAc  <  2:715062/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                         CGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTAC  >  NC_000913/3139447‑3139703

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