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%	77.6 / NA	23	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 (14/9);  total (14/9)

CGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCG  >  NC_000913/3941574‑3941726                                                                             |                                                                             cGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAagagc                                                                           >  1:720776/1‑80 (MQ=255)   gCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAagagcaaa                                                                        <  2:720745/81‑1 (MQ=255)      gagaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCg                                                                     >  2:720775/1‑81 (MQ=255)        gaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCgg                                                                    >  1:720772/1‑80 (MQ=255)           aGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACt                                                                <  1:720746/81‑1 (MQ=255)            gCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTg                                                               <  2:720747/81‑1 (MQ=255)              aaaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGct                                                             >  2:720779/1‑81 (MQ=255)               aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGctc                                                            <  2:720748/81‑1 (MQ=255)                    aaatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTaa                                                       >  2:720767/1‑81 (MQ=255)                     aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTaa                                                       <  2:720749/80‑1 (MQ=255)                         cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAAt                                                   <  2:720750/80‑1 (MQ=255)                                    aGCGGGAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAAt                                       >  1:720783/1‑81 (MQ=255)                                        ggAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgt                                   >  1:720790/1‑81 (MQ=255)                                        ggAAAGCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgt                                   >  1:720769/1‑81 (MQ=255)                                             gCGTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTggg                              >  2:720774/1‑81 (MQ=255)                                               gTATTATGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                            >  1:720773/1‑81 (MQ=255)                                                  ttatGCACACCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                         <  1:720751/81‑1 (MQ=255)                                                       cacaCCCCGCGCCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCg                        >  2:720782/1‑77 (MQ=255)                                                               gcgcCGCGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAtt            <  1:720752/81‑1 (MQ=255)                                                                cgccgcGGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTc           >  1:720778/1‑81 (MQ=255)                                                                   cgcgGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCtt         >  1:720777/1‑80 (MQ=255)                                                                   cgcgGAGAAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTa        <  2:720753/81‑1 (MQ=255)                                                                      ggagaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACg     >  1:720781/2‑81 (MQ=18)                                                                       gagaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGt    <  2:720754/81‑1 (MQ=18)                                                                         gaAAGAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCg  >  1:720786/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 < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑