Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,685	A→G	intergenic (‑358/‑123)	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,685	0	A	G	100.0%	124.7 / NA	48	intergenic (‑358/‑123)	yieP/rrsC	putative transcriptional regulator YieP/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/24);  total (24/24)

GTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGT  >  NC_000913/3941558‑3941782                                                                                                                                |                                                                                                  gTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTgg                                                                                       >  2:405351/1‑140 (MQ=21)  ttCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTggg                                                                                      >  2:651549/1‑140 (MQ=21)  ttCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTggg                                                                                      <  1:1201498/140‑1 (MQ=25)  ttCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTggg                                                                                      <  1:675095/140‑1 (MQ=25)  ttCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTggg                                                                                      >  2:675095/1‑140 (MQ=21)    ctcCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCa                                                                                    <  2:1219694/140‑1 (MQ=25)        tGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCg                                                                                <  1:1103751/140‑1 (MQ=25)          agaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCtgtgt                                                                                         <  1:379528/129‑1 (MQ=25)          agaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCtgtgt                                                                                         >  2:379528/1‑129 (MQ=25)          agaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACt                                                                                  <  1:495881/136‑1 (MQ=25)          agaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACt                                                                                  >  2:495881/1‑136 (MQ=21)              cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATa                                                                          <  2:867/140‑1 (MQ=21)              cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATa                                                                          >  2:1115141/1‑140 (MQ=21)                    cAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGAtt                                                                    >  1:473345/1‑140 (MQ=21)                      gagaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCt                                                                  >  2:44163/1‑140 (MQ=25)                               aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGc                                                         <  1:1257390/140‑1 (MQ=21)                                      atGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGa                                                  >  2:632061/1‑140 (MQ=25)                                      atGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGa                                                  <  2:1370259/140‑1 (MQ=21)                                       tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGaa                                                 <  1:1115141/140‑1 (MQ=21)                                       tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGaa                                                 <  2:1057339/140‑1 (MQ=21)                                            gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATg                                            >  2:148896/1‑140 (MQ=25)                                               tctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAAt                                         >  1:373925/1‑140 (MQ=25)                                                 tGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATAc                                       <  2:555599/140‑1 (MQ=21)                                                 tGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATAc                                       <  2:485952/140‑1 (MQ=21)                                                   tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCa                                     <  2:862424/140‑1 (MQ=21)                                                   tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCa                                     <  1:632061/140‑1 (MQ=21)                                                   tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCa                                     <  2:536342/140‑1 (MQ=21)                                                        ggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtct                                >  1:875407/1‑140 (MQ=25)                                                        ggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtct                                >  1:364412/1‑140 (MQ=21)                                                        ggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGtct                                <  2:1136504/140‑1 (MQ=21)                                                             gCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAaga                           <  2:735727/140‑1 (MQ=17)                                                                 attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGa                       <  2:364412/140‑1 (MQ=17)                                                                 attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGa                       <  2:924399/140‑1 (MQ=17)                                                                 attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGa                       <  1:1011052/140‑1 (MQ=17)                                                                       cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGt                 >  1:488523/1‑140 (MQ=18)                                                                        acacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTa                >  1:994534/1‑140 (MQ=18)                                                                        acacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTa                >  2:199380/1‑140 (MQ=18)                                                                            cccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTc            >  2:84097/1‑140 (MQ=18)                                                                            cccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTc            >  2:958766/1‑140 (MQ=18)                                                                            cccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTc            >  2:163447/1‑140 (MQ=18)                                                                              cgcgCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCAt          >  2:1281850/1‑140 (MQ=18)                                                                               gcgcCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGa                       >  2:522930/1‑126 (MQ=18)                                                                               gcgcCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGa                       <  1:522930/126‑1 (MQ=17)                                                                                cgccgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa        >  1:151180/1‑140 (MQ=18)                                                                                  ccgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg      <  1:84097/140‑1 (MQ=17)                                                                                  ccgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg      <  1:224210/140‑1 (MQ=17)                                                                                    gcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa    >  2:86741/1‑140 (MQ=17)                                                                                      tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGGCAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt  >  1:848613/1‑140 (MQ=18)                                                                                                                                |                                                                                                  GTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGT  >  NC_000913/3941558‑3941782

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