Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,541	G→T	intergenic (‑214/‑267)	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,541	0	G	T	100.0%	121.1 / NA	44	intergenic (‑214/‑267)	yieP/rrsC	putative transcriptional regulator YieP/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (25/19);  total (25/19)

TATGGCACATTAACGGGGCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAA  >  NC_000913/3941415‑3941673                                                                                                                               |                                                                                                                                     tATGGCACATTAACGGGGCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGc                                                                                                                         >  2:359143/1‑140 (MQ=255)     gCACATTAACGGGGCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCgggg                                                                                                                     <  1:336125/140‑1 (MQ=255)       acaTTAACGGGGCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGtt                                                                                                                   <  2:743015/140‑1 (MQ=255)       acaTTAACGGGGCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGtt                                                                                                                   <  2:1118136/140‑1 (MQ=255)         aTTAACGGGGCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTtct                                                                                                                 >  2:127919/1‑140 (MQ=255)            aaCGGGGCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCt                                                                                                              <  1:247499/140‑1 (MQ=255)              cGGGGCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTga                                                                                                            >  1:657593/1‑140 (MQ=255)               ggggCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTgag                                                                                                           <  2:125085/140‑1 (MQ=255)                 ggCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGaa                                                                                                         <  2:1071474/140‑1 (MQ=255)                 ggCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGaa                                                                                                         >  1:559973/1‑140 (MQ=255)                   cTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCTGGGTCAGCGGGGTTCTCCTGAGAACt                                                                                                       >  2:109928/1‑140 (MQ=255)                    ttttGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTgag                                                                                                           <  2:758938/135‑1 (MQ=37)                    ttttGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTgag                                                                                                           >  1:758938/1‑135 (MQ=37)                                                                    gTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGt                                                      >  2:1079015/1‑140 (MQ=17)                                                                       aGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtatt                                                   >  1:53531/1‑140 (MQ=17)                                                                       aGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtatt                                                   >  1:968772/1‑140 (MQ=17)                                                                       aGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGCTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGtatt                                                   >  2:60722/1‑140 (MQ=17)                                                                            ggAATAACTCCCTATAATGCGCCCCCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGca                                              <  1:175308/140‑1 (MQ=17)                                                                              aaTAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGcaca                                            >  2:62885/1‑140 (MQ=17)                                                                                       cTATAATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGcgcc                                   >  2:249864/1‑140 (MQ=17)                                                                                        tataATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGcgccg                                  >  2:45205/1‑140 (MQ=17)                                                                                        tataATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGcgccg                                  >  2:734771/1‑140 (MQ=17)                                                                                         ataATGCGCCACCACTGACACGGAACAACGGCAAACGCTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGcgccgc                                 <  2:1021772/140‑1 (MQ=17)                                                                                          taATGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCt                                >  2:351651/1‑140 (MQ=17)                                                                                           aaTGCGCCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTg                               <  1:62885/140‑1 (MQ=17)                                                                                              gcgcCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTgaga                            <  1:45205/140‑1 (MQ=17)                                                                                               cgcCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGaa                           >  1:455148/1‑140 (MQ=17)                                                                                                gcCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGaaa                          >  1:876686/1‑140 (MQ=17)                                                                                                gcCACCACTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGaaa                          >  2:818504/1‑140 (MQ=17)                                                                                                  caccacTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGaaaaa                        >  1:196068/1‑140 (MQ=17)                                                                                                      acTGACACGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAaaagcaa                    <  1:1059482/140‑1 (MQ=18)                                                                                                          acacGGAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCg                >  1:1088195/1‑140 (MQ=21)                                                                                                              ggAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCAc            <  2:968772/140‑1 (MQ=25)                                                                                                              ggAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCAc            <  1:333441/140‑1 (MQ=25)                                                                                                              ggAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCAc            <  1:1116813/140‑1 (MQ=25)                                                                                                              ggAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCAc            <  2:876686/140‑1 (MQ=25)                                                                                                               gAACAACGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt           >  2:1112719/1‑140 (MQ=21)                                                                                                                aacaacGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCAc            >  1:591704/1‑138 (MQ=21)                                                                                                                aacaacGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCAc            <  2:591704/138‑1 (MQ=25)                                                                                                                aacaacGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTg          >  1:372415/1‑140 (MQ=21)                                                                                                                aacaacGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTg          <  1:1079015/140‑1 (MQ=25)                                                                                                                  caacGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGct        >  1:1072255/1‑140 (MQ=21)                                                                                                                     cGGCAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCtt     >  2:988089/1‑140 (MQ=21)                                                                                                                        cAAACACTCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTaa  <  2:53531/140‑1 (MQ=25)                                                                                                                               |                                                                                                                                     TATGGCACATTAACGGGGCTTTTGCTGAAAAAATGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAA  >  NC_000913/3941415‑3941673

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