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%	102.0 / NA	33	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 (18/15);  total (18/15)

GGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCG  >  NC_000913/3941575‑3941726                                                                            |                                                                             ggCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAagagcaa                                                                         <  2:785517/81‑1 (MQ=37)   cAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAagagcaaag                                                                       <  1:785518/81‑1 (MQ=255)    agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAagagcaaagc                                                                      <  1:785520/81‑1 (MQ=255)    agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAagagcaaagc                                                                      <  2:785519/81‑1 (MQ=255)     gagaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCg                                                                     >  1:785571/1‑81 (MQ=255)       gaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCgg                                                                    >  1:785549/1‑80 (MQ=255)         aaGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCAc                                                                 >  2:785576/1‑81 (MQ=255)          aGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACt                                                                <  1:785521/81‑1 (MQ=255)          aGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACt                                                                >  1:785554/1‑81 (MQ=255)            cAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGc                                                              <  1:785522/81‑1 (MQ=255)            cAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGc                                                              >  1:785564/1‑81 (MQ=255)             aaaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGct                                                             >  1:785558/1‑81 (MQ=255)              aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGctc                                                            >  1:785575/1‑81 (MQ=255)                   aaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTaa                                                       >  1:785574/1‑81 (MQ=255)                      tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACaa                                                    <  1:785523/81‑1 (MQ=255)                           gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTAt                                               >  2:785569/1‑81 (MQ=255)                            aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTAt                                               <  1:785524/80‑1 (MQ=255)                             ctctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTAt                                               >  1:785568/1‑79 (MQ=255)                                 gTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAATGCGGCACTGCTCTTTAACAATTTATCAGAc                                          >  2:785556/1‑80 (MQ=255)                                    gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGGGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAAt                                       <  1:785525/80‑1 (MQ=255)                                    gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAAt                                       <  1:785526/80‑1 (MQ=255)                                      gggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtg                                    >  2:785580/1‑81 (MQ=255)                                      gggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtg                                    <  2:785527/81‑1 (MQ=255)                                         aaaGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgtgt                                 >  1:785559/1‑81 (MQ=255)                                          aaGCGTATTATGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgtgtg                                >  2:785578/1‑81 (MQ=37)                                               tattatGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCAc                           <  2:785528/81‑1 (MQ=18)                                                 ttatGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                         <  1:785529/81‑1 (MQ=18)                                                    tGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGaa                      <  2:785531/81‑1 (MQ=18)                                                    tGCACACCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGa                       <  1:785530/80‑1 (MQ=21)                                                         acCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATAc                 <  2:785532/81‑1 (MQ=18)                                                         acCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATAc                 <  2:785533/81‑1 (MQ=18)                                                         acCCCGCGCCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATAc                 >  2:785563/1‑81 (MQ=18)                                                           cccGCGCCGCTGAGAAAGAGCAAAGCGGCACAGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACgg               <  2:785534/81‑1 (MQ=255)                                                              gcgcCGCTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAtt            <  1:785535/81‑1 (MQ=18)                                                                gccgcTGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCt          >  2:785566/1‑81 (MQ=18)                                                                      gagaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGt    >  1:785570/1‑81 (MQ=18)                                                                       agaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTc   >  1:785567/1‑81 (MQ=18)                                                                         aaagaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCg  >  2:785560/1‑80 (MQ=21)                                                                            |                                                                             GGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCG  >  NC_000913/3941575‑3941726

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