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%	182.5 / NA	50	E212K (GAA→AAA)	nusA	transcription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/22);  total (28/22)

CCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGC  >  NC_000913_3_bme_pgi/3316746‑3317023                                                                                                                                                    |                                                                                                                                   ccTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTaa                                                                                                                                  >  4:401456/1‑150 (MQ=255) ccTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTaa                                                                                                                                  <  1:369723/150‑1 (MQ=255)        acGCGCGCCACGCATACCTACGCAANCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGaa                                                                                                                           >  3:360533/1‑150 (MQ=255)        acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGaa                                                                                                                           >  1:629314/1‑150 (MQ=255)          gcgcgcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACa                                                                                                                         >  1:639666/1‑150 (MQ=255)             cgcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGtt                                                                                                                      <  4:392509/150‑1 (MQ=255)               ccACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACCGTTCg                                                                                                                    >  4:419535/1‑150 (MQ=255)                     aTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                              <  1:681030/150‑1 (MQ=255)                           aCGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTcc                                                                                                        >  1:680928/1‑150 (MQ=255)                            cGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCg                                                                                                       >  3:222259/1‑150 (MQ=255)                             gCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCgg                                                                                                      <  2:657194/150‑1 (MQ=255)                                     cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAAc                                                                                              <  2:673009/150‑1 (MQ=255)                                     cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAAc                                                                                              >  1:366463/1‑150 (MQ=255)                                       aCCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGa                                                                                            <  4:360533/150‑1 (MQ=255)                                        ccGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                           <  2:366463/150‑1 (MQ=255)                                         cGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCCAGCCGCTGCTGTAATTTCAATCACTTCTTCGCCGACTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAgg                                                                                          >  1:590649/1‑149 (MQ=255)                                                                ttttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcgcg                                                                   <  4:195511/150‑1 (MQ=255)                                                                 tttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACgcgc                                                                  >  3:39373/1‑150 (MQ=255)                                                                     aCCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTcc                                                              >  3:30554/1‑150 (MQ=255)                                                                      ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCg                                                             <  1:460967/150‑1 (MQ=255)                                                                       cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCCCGCGCTTCCgg                                                            >  4:268130/1‑150 (MQ=255)                                                                       cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                            <  4:373375/150‑1 (MQ=255)                                                                       cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                            >  1:635969/1‑150 (MQ=255)                                                                       cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                            >  3:114828/1‑150 (MQ=255)                                                                        gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                           >  2:159714/1‑150 (MQ=255)                                                                           aTTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGa                                                        <  4:169930/150‑1 (MQ=255)                                                                            ttttCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGaa                                                       >  2:15889/1‑150 (MQ=255)                                                                                cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGa                                                   <  1:15830/150‑1 (MQ=255)                                                                                cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGa                                                   <  2:429998/150‑1 (MQ=255)                                                                                cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAACATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGa                                                   <  1:126927/150‑1 (MQ=255)                                                                                  cACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAAt                                                 >  2:510870/1‑150 (MQ=255)                                                                                    cGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCt                                                                                                    <  2:313899/97‑1 (MQ=255)                                                                                    cGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCt                                                                                                    >  1:313899/1‑97 (MQ=255)                                                                                     gagaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCtt                                                                                                   <  3:245171/97‑1 (MQ=255)                                                                                     gagaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCtt                                                                                                   >  4:245171/1‑97 (MQ=255)                                                                                                   cgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgcgcg                                                                          >  1:216319/1‑108 (MQ=255)                                                                                                   cgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgcgcg                                                                          <  2:216319/108‑1 (MQ=255)                                                                                                    gAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcgc                               >  1:304302/1‑150 (MQ=255)                                                                                                              tttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAgg                     >  2:390179/1‑150 (MQ=255)                                                                                                               ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAggg                    <  3:401456/150‑1 (MQ=255)                                                                                                                  aTTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAAc                                                      >  4:88860/1‑113 (MQ=255)                                                                                                                  aTTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAAc                                                      <  3:88860/113‑1 (MQ=255)                                                                                                                  aTTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCgg                 >  4:400001/1‑150 (MQ=255)                                                                                                                            cttcttCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGaa                                                                                    >  4:117506/1‑73 (MQ=255)                                                                                                                            cttcttCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGaa                                                                                    <  3:117506/73‑1 (MQ=255)                                                                                                                               cttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAcgcg    >  1:341295/1‑150 (MQ=255)                                                                                                                                 tcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgc  >  1:55140/1‑150 (MQ=255)                                                                                                                                 tcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgc  <  2:27949/150‑1 (MQ=255)                                                                                                                                       aTTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATagag                                             <  4:216391/101‑1 (MQ=255)                                                                                                                                       aTTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATagag                                             >  3:216391/1‑101 (MQ=255)                                                                                                                                                    |                                                                                                                                   CCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGC  >  NC_000913_3_bme_pgi/3316746‑3317023

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