Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,537,286	T→C	*159Q (TAA→CAA)	greB →	transcription elongation factor GreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,537,286	0	T	C	100.0%	164.0 / NA	47	*159Q (TAA→CAA)	greB	transcription elongation factor GreB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/20);  total (27/20)

TAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGG  >  NC_000913/3537153‑3537419                                                                                                                                      |                                                                                                                                      tAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGaa                                                                                                                                 <  2:6989/140‑1 (MQ=255)              tctATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAAcctc                                                                                                                    <  2:650153/140‑1 (MQ=255)               ctATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAAcctcc                                                                                                                   >  1:585448/1‑140 (MQ=255)                tATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGaa                                                                                                                                 <  2:974449/125‑1 (MQ=255)                tATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGaa                                                                                                                                 >  1:974449/1‑125 (MQ=255)                tATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAAcctcct                                                                                                                  >  1:732010/1‑140 (MQ=255)                  tCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGc                                                                                                                >  1:847051/1‑140 (MQ=255)                        ccccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGAtgg                                                                                                          >  2:101595/1‑140 (MQ=255)                          ccGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGAtggct                                                                                                        >  2:151805/1‑140 (MQ=255)                               ggCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCAt                                                                                                   <  2:847051/140‑1 (MQ=255)                                          tGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTgccagcc                                                                                        >  1:433316/1‑140 (MQ=255)                                           gCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTgccagcca                                                                                       >  1:980341/1‑140 (MQ=255)                                            cTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGccagccag                                                                                      <  2:433316/140‑1 (MQ=255)                                             tGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGt                                                                                     <  1:151805/140‑1 (MQ=255)                                             tGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGt                                                                                     <  2:980341/140‑1 (MQ=255)                                             tGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGt                                                                                     <  2:385854/140‑1 (MQ=255)                                             tGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGt                                                                                     <  1:101595/140‑1 (MQ=255)                                             tGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGt                                                                                     <  1:290552/140‑1 (MQ=255)                                               aaaaaaGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTcc                                                                                   >  1:918454/1‑140 (MQ=255)                                                    agaagGCGGCGATCTGGCGGTGGAGAATACCCCTGCCGGGGAACCGAGCTGGTATGCCAATGCTATCGAGTACGAGAAACCGCAAAGAAGAGTCTTAACCTCCTGCCGATTGCGGGCCTTTTTGCCAGCCAGTCCGtata                                                                              >  2:891951/1‑140 (MQ=255)                                                       agTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACt                                                                           <  2:787073/140‑1 (MQ=255)                                                              gATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCt                                                                    <  2:27471/140‑1 (MQ=255)                                                                       gtggtgAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGa                                                           >  2:451669/1‑140 (MQ=255)                                                                          gtgAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATcc                                                        >  2:843405/1‑140 (MQ=255)                                                                                cccccTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAg                                                  <  2:881799/139‑1 (MQ=255)                                                                                  cccTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGAt                                                >  2:236511/1‑140 (MQ=255)                                                                                    cTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGa                                              >  2:192088/1‑140 (MQ=255)                                                                                       cgggggAAGCGAGTTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACt                                           <  2:731592/138‑1 (MQ=255)                                                                                             aaGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAAcc                                     <  1:843405/140‑1 (MQ=255)                                                                                                 gAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCatga                                 >  1:495342/1‑140 (MQ=255)                                                                                                    cTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCAtgatga                              >  1:564943/1‑140 (MQ=255)                                                                                                    cTGGTATGTTAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCAtgatga                              <  2:564943/140‑1 (MQ=255)                                                                                                             tAATGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGt                     >  1:221988/1‑140 (MQ=255)                                                                                                               aTGCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTc                   >  1:442120/1‑140 (MQ=255)                                                                                                                 gCTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTg                 >  1:870175/1‑140 (MQ=255)                                                                                                                  cTATCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGc                >  2:978156/1‑140 (MQ=255)                                                                                                                  cTATCAAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGc                >  2:329158/1‑140 (MQ=255)                                                                                                                    aTCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCg              <  2:495342/140‑1 (MQ=255)                                                                                                                     tcaaGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGc             <  1:329158/137‑1 (MQ=255)                                                                                                                     tCGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGc             >  1:605484/1‑140 (MQ=255)                                                                                                                      cGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGtt                    <  2:718410/132‑1 (MQ=255)                                                                                                                      cGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGtt                    >  1:718410/1‑132 (MQ=255)                                                                                                                      cGAGTACGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCa            >  2:923841/1‑140 (MQ=255)                                                                                                                            cGTGAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTg      >  2:910211/1‑140 (MQ=255)                                                                                                                               gAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCgg   >  2:634663/1‑140 (MQ=255)                                                                                                                               gAAACCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCgg   >  2:299073/1‑140 (MQ=255)                                                                                                                                  aCCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCggg  >  2:468175/1‑138 (MQ=255)                                                                                                                                  aCCGCAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCggg  <  1:468175/138‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      TAAAGATTACATCTCTATCGATTCCCCGATGGCCCGCGCATTGCTGAAAAAAGAAGTCGGCGATCTGGCGGTGGTGAATACCCCTGCCGGGGAAGCGAGCTGGTATGTTAATGCTATCGAGTACGTGAAACCGTAAGGAAGAGTCTTAACCTCCTGCCGATGGCTGGCATTTTTGCCAGCCAGTCCGTATAACTATCCCCTGATTTTTGATCCGAAAAGATGAACTCAAACCATGATGAATGATTCGTTCTGCCGCATTATTGCGGG  >  NC_000913/3537153‑3537419

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