Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,316,893	C→T	E212K (GAA→AAA)	nusA ←	transcription termination/antitermination L factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	3,316,893	0	C	T	100.0%	190.6 / NA	53	E212K (GAA→AAA)	nusA	transcription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/25);  total (28/25)

ACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGC  >  NC_000913_3_bme_pgi/3316753‑3317036                                                                                                                                             |                                                                                                                                                acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGaa                                                                                                                                        >  1:399259/1‑150 (MQ=255)           cGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATc                                                                                                                              <  2:107194/150‑1 (MQ=255)            gCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCa                                                                                                                             >  4:83628/1‑150 (MQ=255)                         aGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCtt                                                                                                                <  3:265309/150‑1 (MQ=255)                         aGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCtt                                                                                                                <  4:33226/150‑1 (MQ=255)                            aCCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCtt                                                                                                                >  4:110503/1‑147 (MQ=255)                            aCCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCtt                                                                                                                <  3:110503/147‑1 (MQ=255)                              cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGa                                                                                                                                <  1:123977/129‑1 (MQ=255)                              cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGa                                                                                                                                >  2:123977/1‑129 (MQ=255)                              cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAAc                                                                                                           >  1:100565/1‑150 (MQ=255)                                      tCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAGCCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACg                                                                                                   <  1:402798/150‑1 (MQ=255)                                      tCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACg                                                                                                   >  4:210271/1‑150 (MQ=255)                                                       ggTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcg                                                                                  >  4:157588/1‑150 (MQ=255)                                                              aCCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTcc                                                                           >  1:209399/1‑150 (MQ=255)                                                                 gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                        >  1:226848/1‑150 (MQ=255)                                                                 gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                        >  3:115119/1‑150 (MQ=255)                                                                      tttCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAAc                                                                   <  4:87974/150‑1 (MQ=255)                                                                          gCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGaa                                                               <  3:210271/150‑1 (MQ=255)                                                                            aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcg                                                                                  <  2:177539/129‑1 (MQ=255)                                                                            aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcg                                                                                  >  1:177539/1‑129 (MQ=255)                                                                                gaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                        <  2:438638/103‑1 (MQ=255)                                                                                gaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                        <  2:438701/103‑1 (MQ=255)                                                                                gaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                        >  1:438701/1‑103 (MQ=255)                                                                                gaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                        >  1:438638/1‑103 (MQ=255)                                                                                    ccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGAc                                                                                                    <  3:90652/103‑1 (MQ=255)                                                                                    ccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGAc                                                                                                    >  4:90652/1‑103 (MQ=255)                                                                                    ccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGc                                                     >  3:115564/1‑150 (MQ=255)                                                                                    ccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGc                                                     >  3:145325/1‑150 (MQ=255)                                                                                                     gctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCa                                    <  4:115119/150‑1 (MQ=255)                                                                                                        ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                        >  1:200314/1‑79 (MQ=255)                                                                                                        ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                        <  2:200314/79‑1 (MQ=255)                                                                                                        ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAggg                                 <  2:209399/150‑1 (MQ=255)                                                                                                        ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAggg                                 <  3:83628/150‑1 (MQ=255)                                                                                                        ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAggg                                 <  1:421460/150‑1 (MQ=255)                                                                                                         tAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGc                                >  4:277137/1‑150 (MQ=255)                                                                                                          aaTTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCg                               >  2:209636/1‑150 (MQ=255)                                                                                                           aTTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCgg                              <  2:226848/150‑1 (MQ=255)                                                                                                            tttCAATCACTTCTTCCCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGTACGCCACGAACGCGGTCGCCAGGGCGGa                             >  3:257312/1‑150 (MQ=255)                                                                                                             ttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGaa                            >  3:17599/1‑150 (MQ=255)                                                                                                             ttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGaa                            <  4:64704/150‑1 (MQ=255)                                                                                                              tCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAg                           >  2:220150/1‑150 (MQ=255)                                                                                                               cAATCACTTCGTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGt                          >  4:218/1‑150 (MQ=255)                                                                                                                 aTCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGttt                        >  3:59533/1‑150 (MQ=255)                                                                                                                       tcttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAcgc                  <  2:407444/150‑1 (MQ=255)                                                                                                                        cttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAcgcg                 <  1:209636/150‑1 (MQ=255)                                                                                                                            gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCa                                                                                    <  4:131711/79‑1 (MQ=255)                                                                                                                            gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCa                                                                                    >  3:131711/1‑79 (MQ=255)                                                                                                                            gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACg                                                                  <  2:255842/97‑1 (MQ=255)                                                                                                                            gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACg                                                                  >  1:255842/1‑97 (MQ=255)                                                                                                                            gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                             <  2:367547/102‑1 (MQ=255)                                                                                                                            gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                             >  1:367547/1‑102 (MQ=255)                                                                                                                             ccGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgcagc            >  2:234484/1‑150 (MQ=255)                                                                                                                                       gCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTcgc  <  1:139586/150‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                ACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGC  >  NC_000913_3_bme_pgi/3316753‑3317036

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