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	100.0%	110.5 / NA	33	R286H (CGT→CAT)	yghU	putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/16);  total (17/16)

CGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACG  >  NC_000913/3139447‑3139704                                                                                                                                         |                                                                                                                          cGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATc                                                                                                                         <  2:588878/139‑1 (MQ=255) cGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATc                                                                                                                         <  1:360342/139‑1 (MQ=255)              aagTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGtt                                                                                                            <  2:526372/139‑1 (MQ=255)                                     aaaCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGAGTTCGTCGCAGCAAGCCATc                                                                                     >  2:390594/1‑139 (MQ=255)                                        cGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAg                                                                                  <  2:446037/139‑1 (MQ=255)                                          tGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGc                                                                                >  1:200208/1‑139 (MQ=255)                                                            cACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgc                                                              >  1:264694/1‑139 (MQ=255)                                                            cACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgc                                                              >  1:491534/1‑139 (MQ=255)                                                            cACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgc                                                              >  2:617160/1‑139 (MQ=255)                                                             aCCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGc                                                                                <  1:443510/120‑1 (MQ=255)                                                             aCCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGc                                                                                >  2:443510/1‑120 (MQ=255)                                                             aCCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgcc                                                             >  1:27068/1‑139 (MQ=255)                                                                     aCCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                     <  1:617160/139‑1 (MQ=255)                                                                         cTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgccg                                                            >  2:499344/1‑128 (MQ=255)                                                                         cTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgccg                                                            <  1:429302/128‑1 (MQ=255)                                                                         cTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgccg                                                            <  1:499344/128‑1 (MQ=255)                                                                         cTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCAcgccg                                                            >  2:429302/1‑128 (MQ=255)                                                                         cTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt                                                 >  1:151286/1‑139 (MQ=255)                                                                               gAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGca                                           <  2:200208/139‑1 (MQ=255)                                                                                              cgcCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGa                            >  2:10887/1‑139 (MQ=255)                                                                                               gcCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAg                                                                                  >  1:539152/1‑84 (MQ=255)                                                                                               gcCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAg                                                                                  <  2:539152/84‑1 (MQ=255)                                                                                                     aCGCCAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGcc                     <  2:151286/139‑1 (MQ=255)                                                                                                        ccAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGc                      <  2:230243/135‑1 (MQ=255)                                                                                                        ccAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGc                      >  1:230243/1‑135 (MQ=255)                                                                                                        ccAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCcgcg                  <  1:10887/139‑1 (MQ=255)                                                                                                         cAGTGATTTCGAGACGAATACGGAAGATAAGCATCAGGGGTAAGGGTAGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGt                 <  1:379823/139‑1 (MQ=255)                                                                                                                      aCGAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTa    >  1:422485/1‑139 (MQ=255)                                                                                                                        gAATACGGAAGATAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACg  >  1:403289/1‑139 (MQ=255)                                                                                                                                   aTAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCg                                               >  1:8038/1‑83 (MQ=255)                                                                                                                                   aTAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCg                                               <  2:8038/83‑1 (MQ=255)                                                                                                                                   aTAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACg                          >  1:578408/1‑104 (MQ=255)                                                                                                                                   aTAAGCATCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACg                          <  2:578408/104‑1 (MQ=255)                                                                                                                                         |                                                                                                                          CGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACG  >  NC_000913/3139447‑3139704

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