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%	107.9 / NA	32	P5S (CCT→TCT)	nusG	transcription termination factor NusG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/16);  total (16/16)

CCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCAT  >  NC_000913/4177678‑4177829                                                                              |                                                                           ccTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCt                                                                          >  2:761584/1‑80 (MQ=255)         gACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAgcgc                                                                 <  2:761547/81‑1 (MQ=255)         gACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAgcgc                                                                 >  1:761609/1‑81 (MQ=255)          aCTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCt                                                                <  2:767363/81‑1 (MQ=255)           cTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTg                                                               >  2:761577/1‑81 (MQ=255)           cTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTg                                                               >  1:761587/1‑81 (MQ=255)              gATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTa                                                            >  2:761581/1‑81 (MQ=255)                     ttCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGtt                                                     <  1:761548/81‑1 (MQ=255)                         ggTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGATGCTTCTAAAAAGCGCTGGTACGTCGTTCAgg                                                 >  2:761595/1‑81 (MQ=255)                          gTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAgg                                                 <  1:761549/80‑1 (MQ=255)                          gTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGc                                                >  1:761606/1‑81 (MQ=255)                              gCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGttt                                            >  2:761590/1‑81 (MQ=255)                              gCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGtt                                             <  2:761550/80‑1 (MQ=255)                                 tGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTcc                                         >  2:761600/1‑81 (MQ=255)                                   gTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTcc                                         <  2:761551/79‑1 (MQ=255)                                    tATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGt                                      >  1:761611/1‑81 (MQ=255)                                     aTCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGtt                                     >  2:761599/1‑81 (MQ=255)                                            aTCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAg                               <  1:761552/80‑1 (MQ=255)                                                    ccTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAg                      <  2:761553/81‑1 (MQ=255)                                                    ccTGAGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAg                      >  2:761619/1‑81 (MQ=255)                                                        aGGTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCaa                   <  1:761554/80‑1 (MQ=255)                                                          gTTCTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGt                <  1:761555/81‑1 (MQ=255)                                                             cTGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCg              <  1:761556/80‑1 (MQ=255)                                                              tGAGATGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCt            >  2:761593/1‑81 (MQ=255)                                                               gagaTGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTg           >  2:761602/1‑81 (MQ=255)                                                                agaTGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGc          <  1:761557/81‑1 (MQ=255)                                                                 gaTGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCg         <  2:761558/81‑1 (MQ=255)                                                                  aTGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGt        <  2:761559/81‑1 (MQ=255)                                                                   tGTCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGt        >  1:761605/1‑80 (MQ=255)                                                                     tCTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGa      <  2:761560/80‑1 (MQ=255)                                                                      cTGAAGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAg     >  2:761596/1‑80 (MQ=255)                                                                         aaGCTTCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCat  <  2:761561/80‑1 (MQ=255)                                                                              |                                                                           CCTGTGGGGACTGGATGGTATTCTGGTTCGCCTGGTATCCTTTATCACTGGCCTGAGGTTCTGAGATGTCTGAAGCTCCTAAAAAGCGCTGGTACGTCGTTCAGGCGTTTTCCGGTTTTGAAGGCCGCGTAGCAACGTCGCTGCGTGAGCAT  >  NC_000913/4177678‑4177829

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