Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	744,025	C→T	Q145* (CAA→TAA)	pxpB →	5‑oxoprolinase component B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	744,025	0	C	T	100.0%	110.1 / NA	31	Q145* (CAA→TAA)	pxpB	5‑oxoprolinase component B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/18);  total (13/18)

TGGAATACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAACAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCTGTCGGGATCGGCGGGCCGCAGA  >  NC_000913/743960‑744100                                                                  |                                                                            tGGAATACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCAc                                                              <  1:156980/81‑1 (MQ=255) tGGAATACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCAc                                                              <  2:156982/81‑1 (MQ=255) tGGAATACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCAc                                                              <  2:156981/81‑1 (MQ=255) tGGAATACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCAc                                                              <  1:156979/81‑1 (MQ=255)  ggAATACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACg                                                             >  2:157020/1‑81 (MQ=255)   gAATACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACg                                                             <  2:156983/80‑1 (MQ=255)    aaTACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGgc                                                           <  1:156984/81‑1 (MQ=255)      tACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGgcgc                                                         >  2:157045/1‑81 (MQ=255)       aCGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGgcgcg                                                        >  1:157012/1‑81 (MQ=255)          tGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTg                                                     <  2:156985/81‑1 (MQ=255)            gTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGaa                                                   <  1:156986/81‑1 (MQ=255)                ggTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACcgc                                               >  2:157015/1‑81 (MQ=255)                 gTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAAcc                                                 <  1:156987/78‑1 (MQ=255)                 gTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACcgcg                                              >  1:157042/1‑81 (MQ=255)                     ttAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCtt                                           <  1:156988/80‑1 (MQ=255)                        ggTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACt                                        <  2:156989/80‑1 (MQ=255)                        ggTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACt                                        <  2:156990/80‑1 (MQ=255)                         gTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCg                                      >  1:157038/1‑81 (MQ=255)                          ttttCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGt                                     <  2:156991/81‑1 (MQ=255)                           tttCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACATGCCACGGCGCGCTGAACCGCGCTTACTCGtt                                    >  2:157033/1‑81 (MQ=255)                            ttCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGtt                                    <  2:156992/80‑1 (MQ=255)                                aCCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCgg                                >  2:157019/1‑80 (MQ=255)                                    ggCTTCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGt                           >  2:157035/1‑81 (MQ=255)                                       ttCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCt                        >  2:157050/1‑81 (MQ=255)                                        tCCCGTATCTCGGGAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCTg                       >  1:157055/1‑81 (MQ=255)                                                   gggAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCTGTCGGGATcggc            >  1:157044/1‑81 (MQ=255)                                                    ggAGTTTGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCTGTCGGGATcggcg           >  1:157041/1‑81 (MQ=255)                                                          tGCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCTGTCGGGATCGGCGGGCCg      <  1:156993/80‑1 (MQ=255)                                                           gCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCTGTCGGGATCGGCGGGCCGc     <  1:156994/80‑1 (MQ=255)                                                           gCCGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCTGTCGGGATCGGCGGGCCGc     <  2:156995/80‑1 (MQ=255)                                                             cGGAATAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCTGTCGGGATCGGCGGGCCGCAGa  <  2:156996/81‑1 (MQ=255)                                                                  |                                                                            TGGAATACGTGGTCTGGTTTTTAGGTTTTCAACCGGGCTTCCCGTATCTCGGGAGTTTGCCGGAACAACTACACACGCCACGGCGCGCTGAACCGCGCTTACTCGTTCCGGCAGGTTCTGTCGGGATCGGCGGGCCGCAGA  >  NC_000913/743960‑744100

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