Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,650	A→G	intergenic (‑323/‑158)	yieP ← / → rrsC	putative transcriptional regulator YieP/16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,650	0	A	G	100.0%	85.2 / NA	29	intergenic (‑323/‑158)	yieP/rrsC	putative transcriptional regulator YieP/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/16);  total (13/16)

AGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTC  >  NC_000913/3941580‑3941725                                                                       |                                                                            agaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCgg                                                                   <  2:718333/81‑1 (MQ=255)  gaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCgg                                                                   >  2:718400/1‑80 (MQ=255)     aGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGTACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACt                                                               >  2:718393/1‑81 (MQ=255)     aGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACt                                                               <  2:718334/81‑1 (MQ=255)          aaataaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGctct                                                          >  2:718391/1‑81 (MQ=255)                   cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAAtt                                                 <  1:718336/81‑1 (MQ=255)                   cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAAtt                                                 >  2:718385/1‑81 (MQ=255)                      gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTAt                                              >  2:718395/1‑81 (MQ=255)                       aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATc                                             <  2:718336/81‑1 (MQ=255)                          ctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGa                                          >  1:718397/1‑81 (MQ=255)                            gTAGCGGGAAAGCGTATTATGCACACTCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                                         >  1:718383/1‑80 (MQ=255)                            gTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                                         >  1:718385/1‑80 (MQ=255)                             tAGCGGGAAAGCGTATTATGCACGCCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                                         <  2:718338/79‑1 (MQ=255)                             tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                                         <  1:718338/79‑1 (MQ=255)                                cGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCt                                    <  2:718339/81‑1 (MQ=255)                                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgtg                                 <  1:718341/81‑1 (MQ=255)                                       gCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTggg                             <  1:718342/81‑1 (MQ=18)                                          tattatGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCAc                          <  1:718343/81‑1 (MQ=18)                                           attatGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACt                         >  1:718395/1‑81 (MQ=18)                                           attatGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACt                         >  1:718388/1‑81 (MQ=18)                                              atGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCg                       <  2:718343/80‑1 (MQ=21)                                                  acacCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGa                   >  2:718398/1‑80 (MQ=21)                                                  acacCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGAt                  <  2:718344/81‑1 (MQ=18)                                                      cccGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg               >  2:718401/1‑80 (MQ=21)                                                          cgccgcTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTc          <  2:718345/81‑1 (MQ=18)                                                            ccgcTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCtt        <  2:718346/81‑1 (MQ=18)                                                                tGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACg    >  1:718414/1‑81 (MQ=18)                                                                  agaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGt   <  2:718348/80‑1 (MQ=21)                                                                  agaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTc  <  1:718348/81‑1 (MQ=18)                                                                       |                                                                            AGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTC  >  NC_000913/3941580‑3941725

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