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%	60.4 / NA	22	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/9);  total (13/9)

CGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCG  >  NC_000913/3941574‑3941726                                                                             |                                                                             cGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAagagc                                                                           <  2:770838/80‑1 (MQ=21)   gCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAagagcaaa                                                                        <  1:770839/81‑1 (MQ=255)        gaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGc                                                                   <  2:770840/81‑1 (MQ=255)         aaaGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCa                                                                  >  2:770874/1‑81 (MQ=255)             cAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGc                                                              <  2:770841/81‑1 (MQ=255)                aaataaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGctct                                                           >  2:770879/1‑81 (MQ=255)                   taaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCttt                                                         >  1:770871/1‑80 (MQ=255)                           tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTa                                                >  1:770882/1‑81 (MQ=255)                                 tGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                            <  2:770842/79‑1 (MQ=255)                                       gggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtg                                    >  1:770883/1‑81 (MQ=255)                                           aaGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgtgtg                                >  1:770881/1‑81 (MQ=37)                                               gTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                            >  1:770886/1‑81 (MQ=18)                                                tattatGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCAc                           >  2:770891/1‑81 (MQ=18)                                                  ttatGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACt                          <  1:770844/80‑1 (MQ=21)                                                  ttatGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                         <  1:770843/81‑1 (MQ=18)                                                    atGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCg                        <  2:770845/80‑1 (MQ=21)                                                     tGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGa                       <  1:770846/80‑1 (MQ=21)                                                         cacCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATa                  >  1:770900/1‑81 (MQ=18)                                                                cgccgcTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTc           >  1:770880/1‑81 (MQ=18)                                                                   cgcTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCtt         >  1:770876/1‑80 (MQ=21)                                                                      tGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACg     >  2:770901/1‑81 (MQ=18)                                                                        agaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTc   <  1:770847/81‑1 (MQ=18)                                                                         gaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCg  <  2:770848/81‑1 (MQ=18)                                                                         gaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCg  >  2:770887/1‑81 (MQ=18)                                                                             |                                                                             CGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCG  >  NC_000913/3941574‑3941726

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