Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	1,196,093	C→T	P324L (CCT→CTT)	icd →	isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	1,196,093	0	C	T	100.0%	122.8 / NA	40	P324L (CCT→CTT)	icd	isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/19);  total (21/19)

CGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAA  >  NC_000913_3_bme_pgi/1195947‑1196235                                                                                                                                                   |                                                                                                                                               cGTCATTAAAGCCGTGATTGCTGATGCCTTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCTCTGTCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTg                                                                                                                                              <  1:237917/149‑1 (MQ=255) cGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGCCTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTg                                                                                                                                              <  2:645496/149‑1 (MQ=255)  gTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTgg                                                                                                                                             <  1:535809/149‑1 (MQ=255)   tCATTAAAGACGTGCTTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGt                                                                                                                                            <  1:202396/149‑1 (MQ=255)    cATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCATGTTGGCGGTATCGGTATCGCCCTTGGTg                                                                                                                                           <  1:542670/149‑1 (MQ=255)         aaGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGAAGTTATCGCCTGGATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTGTCGGTATCGCCCTTGGTGCAAAc                                                                                                                                      >  1:561814/1‑149 (MQ=255)            aCGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCCGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATc                                                                                                                                   >  1:436686/1‑149 (MQ=255)               tGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCGGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGt                                                                                                                                >  1:108232/1‑149 (MQ=255)                     cgtaTGCATTCCTGCAACAGATCCTGCTGCGTCGGGCAGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACGACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGaa                                                                                                                          <  2:417425/146‑1 (MQ=255)                          gCATTACCGCAACAGATTCTGCTGCGCCCGGCTCAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACCACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATgcgc                                                                                                                     <  2:18644/149‑1 (MQ=255)                                  gCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTg                                                                                                             >  1:640185/1‑149 (MQ=255)                                    aaCAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGaa                                                                                                           >  1:210863/1‑149 (MQ=255)                                    aaCAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGaa                                                                                                           <  2:640330/149‑1 (MQ=255)                                     aCAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAg                                                                                                          >  1:520875/1‑149 (MQ=255)                                          tCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGccacc                                                                                                     <  2:21631/149‑1 (MQ=255)                                            ctgccgCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGccacccc                                                                                                   >  1:413112/1‑148 (MQ=255)                                                   gTCCGGCTGAATATTATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACt                                                                                            <  1:263746/149‑1 (MQ=255)                                                            aaTATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGaaa                                                                                   <  1:517840/149‑1 (MQ=255)                                                                 gatgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGc                                                                              >  2:572830/1‑149 (MQ=255)                                                                   tgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCCAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCAAGGTACTGCGCCGTAATATGcgg                                                                            >  2:544853/1‑147 (MQ=255)                                                                   tgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCg                                                                            <  2:436686/149‑1 (MQ=255)                                                                     ttATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGACTGCGCCCTGTTTGAAGCCACCCACGGTACTCCGCCGAAGTATGcgtgt                                                                          >  2:481905/1‑145 (MQ=255)                                                                        tCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCg                                                                                                                                  <  1:562979/90‑1 (MQ=255)                                                                        tCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCg                                                                                                                                  >  2:562979/1‑90 (MQ=255)                                                                          gCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTAAGGa                                                                     >  2:388565/1‑149 (MQ=255)                                                                                 tGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTa                                                              >  1:412651/1‑149 (MQ=255)                                                                                    aCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTa                                                                                              >  2:32607/1‑114 (MQ=255)                                                                                    aCCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTa                                                                                              <  1:32607/114‑1 (MQ=255)                                                                                      cTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATcc                                                         >  1:265911/1‑149 (MQ=255)                                                                                                    cATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACAAGCAGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATtc                                           >  1:159619/1‑149 (MQ=255)                                                                                                     aTTTGTGCCGCCCTGGCAGCGCAGGTTTGCGGTATCGGTATCGCCATTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCAGGGTACTTCGCACGAATATTCCTGCTACGCGAAAGTCACTCCTGGCTCTTTTACtcc                                          >  2:153547/1‑148 (MQ=255)                                                                                                     aTTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTAACGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAATAAATCCCAGCTCTATTATtct                                          >  2:312437/1‑149 (MQ=255)                                                                                                     aTTTCTCACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATtct                                          <  1:87505/149‑1 (MQ=255)                                                                                                                 cTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATAATTCTCTCCGCTGAgat                              >  1:86129/1‑149 (MQ=255)                                                                                                                       gcgcAGGTTGGCGGTATCGGTATCGCCCTTGATGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCCGAGATGAtgct                        >  2:43839/1‑149 (MQ=255)                                                                                                                           aGGTTGGCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTgcgc                    <  1:445017/149‑1 (MQ=255)                                                                                                                                 gCGGTATCGGTATCGCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATg              <  1:572830/149‑1 (MQ=255)                                                                                                                                  cagtagcggtatcgTCGTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGGCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg             <  2:471965/147‑1 (MQ=255)                                                                                                                                     tatcggtatcgCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGAGCCACATGGGtt          >  1:282418/1‑149 (MQ=255)                                                                                                                                             tcgCCCTTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGaa  <  1:630441/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               CGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCAGCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAA  >  NC_000913_3_bme_pgi/1195947‑1196235

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