Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,177,755	C→T	P5S (CCT→TCT)	nusG →	transcription termination factor NusG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,177,755	0	C	T	100.0%	148.4 / NA	41	P5S (CCT→TCT)	nusG	transcription termination factor NusG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/19);  total (22/19)

CCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATC  >  NC_000913/4177678‑4177832                                                                              |                                                                              ccTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCt                                                                             >  1:828424/1‑80 (MQ=255)  cTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTa                                                                            <  2:828378/80‑1 (MQ=255)     tgGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTaaa                                                                          <  2:828379/79‑1 (MQ=255)       gggACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAgc                                                                      <  2:828380/81‑1 (MQ=255)               aTGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTAc                                                              <  2:828381/81‑1 (MQ=255)                 ggTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTAcgt                                                            >  1:828429/1‑81 (MQ=255)                  gTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTAcgtc                                                           >  2:828421/1‑81 (MQ=255)                     ttCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTAcgtcgt                                                         <  2:828383/80‑1 (MQ=255)                     ttCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGtt                                                        <  1:828382/81‑1 (MQ=255)                       cTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTc                                                       <  1:828384/80‑1 (MQ=255)                        tGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAg                                                     >  1:828423/1‑81 (MQ=255)                          gTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGc                                                   >  1:828419/1‑81 (MQ=255)                              gCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGttt                                               >  1:828422/1‑81 (MQ=255)                                cTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTc                                             <  1:828385/81‑1 (MQ=255)                                cTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTc                                             >  1:828448/1‑81 (MQ=255)                                 tGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTcc                                            >  1:828455/1‑81 (MQ=255)                                      tCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGttt                                       <  2:828386/81‑1 (MQ=255)                                                cTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCcgc                             <  1:828387/81‑1 (MQ=255)                                                     cTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAg                         >  1:828427/1‑80 (MQ=255)                                                        aGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCaa                      <  2:828390/80‑1 (MQ=255)                                                        aGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCa                       <  2:828388/79‑1 (MQ=255)                                                        aGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAAc                     >  1:828445/1‑81 (MQ=255)                                                        aGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAAc                     <  2:828389/81‑1 (MQ=255)                                                            tCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTCTTGAAGGCCGCGTAGCAACGTCg                 >  2:828433/1‑81 (MQ=255)                                                             cTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGc                >  1:828446/1‑81 (MQ=255)                                                              tGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCt               >  2:828432/1‑81 (MQ=255)                                                               gagaTGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTg              >  1:828425/1‑81 (MQ=255)                                                                agaTGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCt               <  1:828391/79‑1 (MQ=255)                                                                  aTGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGt           >  2:828426/1‑81 (MQ=255)                                                                     tCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGa         <  2:828392/80‑1 (MQ=255)                                                                      cTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAg        >  2:828443/1‑80 (MQ=255)                                                                      cTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAg        >  1:828393/1‑80 (MQ=255)                                                                       tGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCa      >  2:828437/1‑81 (MQ=255)                                                                        gAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCa      >  1:828449/1‑80 (MQ=255)                                                                         aaGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCata    <  2:828395/81‑1 (MQ=255)                                                                         aaGCTTCTAAAAAGCGCTGGTACGGCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCa      <  1:828394/79‑1 (MQ=255)                                                                          aGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCatat   >  2:828463/1‑81 (MQ=255)                                                                           gCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCatat   >  1:828428/1‑80 (MQ=255)                                                                           gCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCatat   >  2:828434/1‑80 (MQ=255)                                                                           gCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATc  <  1:828397/81‑1 (MQ=255)                                                                           gCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATc  <  1:828396/81‑1 (MQ=255)                                                                              |                                                                              CCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCATATC  >  NC_000913/4177678‑4177832

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: ‑