Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,316,789	C→T	R145H (CGT→CAT)	trpA ←	tryptophan synthase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,316,789	0	C	T	100.0%	111.8 / NA	35	R145H (CGT→CAT)	trpA	tryptophan synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/19);  total (16/19)

GTTTTCTGCGCCGGTCACGCCTGCTCGTGACAGCAAATAGGTGTAACCACGACCGTAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGACGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCAATGCCTTTGTTAAACACCAGATTGG  >  NC_000913/1316662‑1316915                                                                                                                                |                                                                                                                               gTTTTCTGCGCCGGTCACGCCTGCTCGTGACAGCAAATAGGTGTAACCACGACCGCAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGcc                                                                                                                     <  1:1129963/139‑1 (MQ=255)                cACGCCTGCTCGTGACAGCAAATAGGTGTAACCACGACCGTAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGcgcg                                                                                                      <  2:767711/139‑1 (MQ=255)                        cTCGTGACAGCAAATAGGTGTAACCACGACCGTAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTc                                                                                              >  1:1185030/1‑139 (MQ=255)                                aGCAAATAGGTGTAACCACGACCGTAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGca                                                                                      <  1:310984/139‑1 (MQ=255)                                  cAAATAGGTGTAACCACGACCGTAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGcaca                                                                                    <  2:418477/139‑1 (MQ=255)                                     aTAGGTGTAACCACGACCGTAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCg                                                                                 <  2:1185030/139‑1 (MQ=255)                                              aCCACGACCGTAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGc                                                                        >  2:122657/1‑139 (MQ=255)                                                        tAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGgcg                                                                                                        <  1:1569618/97‑1 (MQ=255)                                                        tAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGgcg                                                                                                        >  2:1569618/1‑97 (MQ=255)                                                        tAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGActct                                                                                                >  1:1563906/1‑105 (MQ=255)                                                        tAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGActct                                                                                                <  2:1563906/105‑1 (MQ=255)                                                        tAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGa                                                              >  1:1443985/1‑139 (MQ=255)                                                        tAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGa                                                              >  2:1474156/1‑139 (MQ=255)                                                              gCTATCTGGCGCAGCAGGTCGTCATCGGCATTAGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGa                                                        <  1:213607/139‑1 (MQ=255)                                                                tATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGGCGCCGACt                                                      >  1:596447/1‑139 (MQ=255)                                                                   cTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAACACCGAATCGACGCCGACtttt                                                   >  2:106783/1‑139 (MQ=255)                                                                     ggCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGca                                                                                      <  1:265496/102‑1 (MQ=255)                                                                     ggCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGca                                                                                      >  2:265496/1‑102 (MQ=255)                                                                       cgcAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGc                                               >  1:424824/1‑139 (MQ=255)                                                                          agcagGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACt                                            <  1:1474156/139‑1 (MQ=255)                                                                            cagGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTgg                                          >  1:732302/1‑139 (MQ=255)                                                                                 cgtcATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATa                                     <  2:424824/139‑1 (MQ=255)                                                                                     aTCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGGCGCCGACTTTTTCGCACTGGGCATAAAAc                                 <  2:596447/139‑1 (MQ=255)                                                                                     aTCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAAc                                 <  2:1443985/139‑1 (MQ=255)                                                                                              tGGCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCAATg                        >  1:1001937/1‑139 (MQ=255)                                                                                               ggCGGGCAGATGAAGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCAATGc                       <  1:122657/139‑1 (MQ=255)                                                                                                           aaGATAGGTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCAATGCCTTTGTTAAaca           <  1:607937/139‑1 (MQ=255)                                                                                                                  gTGCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCAATGCCTTTGTTAAACACCAGAtt    >  2:246509/1‑139 (MQ=255)                                                                                                                    gCGACATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCAATGCCTTTGTTAAACACCAGATTgg  <  1:246509/139‑1 (MQ=255)                                                                                                                       aCATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACtttt                                                   <  1:304682/87‑1 (MQ=255)                                                                                                                       aCATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACtttt                                                   >  2:304682/1‑87 (MQ=255)                                                                                                                       aCATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCaa                          <  1:834701/112‑1 (MQ=255)                                                                                                                       aCATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCaa                          <  1:398679/112‑1 (MQ=255)                                                                                                                       aCATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCaa                          >  2:398679/1‑112 (MQ=255)                                                                                                                       aCATTATGATGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCaa                          >  2:834701/1‑112 (MQ=255)                                                                                                                                |                                                                                                                               GTTTTCTGCGCCGGTCACGCCTGCTCGTGACAGCAAATAGGTGTAACCACGACCGTAAGAGGCTATCTGGCGCAGCAGGTCGTCATCGGCATTTGGCGGGCAGATGAAGATAGGTGCGACATTATGACGCAACGCGGCCTGGCGGAAGGGCGCGGACTCTTCAACTGGCACATCGGCAACCAGCACCGAATCGACGCCGACTTTTTCGCACTGGGCATAAAACTCATCAATGCCTTTGTTAAACACCAGATTGG  >  NC_000913/1316662‑1316915

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