Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,810,695	T→C	M162T (ATG→ACG)	ydjM →	putative inner membrane protein regulated by LexA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,810,695	0	T	C	100.0%	113.1 / NA	31	M162T (ATG→ACG)	ydjM	putative inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/16);  total (15/16)

ATCCTGGTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCA  >  NC_000913/1810622‑1810770                                                                          |                                                                            aTCCTGGTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCcat                                                                      >  2:332904/1‑81 (MQ=255) aTCCTGGTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCcat                                                                      >  1:332902/1‑81 (MQ=255)  tCCTGGTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCcat                                                                      <  1:332866/80‑1 (MQ=255)  tCCTGGTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCca                                                                       <  1:332868/79‑1 (MQ=255)  tCCTGGTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGAATTATTTGTCTGGTCGGTATGGACGCCCca                                                                       <  1:332867/79‑1 (MQ=255)   ccTGGTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCcatt                                                                     >  2:332905/1‑80 (MQ=255)   ccTGGTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCcatt                                                                     >  1:332898/1‑80 (MQ=255)       gTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCcattcatt                                                                 <  1:332869/80‑1 (MQ=255)       gTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCcattcat                                                                  <  1:332870/79‑1 (MQ=255)           cccAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTAccc                                                             >  1:332901/1‑80 (MQ=255)           cccAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCg                                                            <  1:332871/81‑1 (MQ=255)           cccAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCg                                                            <  1:332872/81‑1 (MQ=255)             cAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGCTCGGTATGGACGCCCCATTCATTAccc                                                             >  1:332889/1‑78 (MQ=255)                aaGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTTGGTATGGACGCCCCATTCATTACCCGAGaa                                                        >  2:332907/1‑80 (MQ=255)                 aGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGaa                                                        >  2:332900/1‑79 (MQ=255)                      aaccaacTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCt                                                 <  1:332873/81‑1 (MQ=255)                          aacTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTc                                             <  2:332874/81‑1 (MQ=255)                                   gTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCAt                                    <  2:332875/81‑1 (MQ=255)                                     tttATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCg                                  >  1:332906/1‑81 (MQ=255)                                       tATCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGc                                 <  1:332876/80‑1 (MQ=255)                                        aTCTGCATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCa                                <  1:332877/80‑1 (MQ=255)                                             cATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCaaa                              <  1:332880/77‑1 (MQ=255)                                             cATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGAt                          <  1:332878/81‑1 (MQ=255)                                             cATGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGAt                          <  2:332879/81‑1 (MQ=255)                                              aTGGCATTATTTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGAt                          >  1:332909/1‑80 (MQ=255)                                                   attattTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATAc                    >  2:332903/1‑81 (MQ=255)                                                    ttattTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATAcc                   >  1:332899/1‑81 (MQ=255)                                                     tattTGTCTGGTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCt                  >  1:332917/1‑81 (MQ=255)                                                               gTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATc          >  1:332908/1‑79 (MQ=255)                                                               gTCGGTATGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATc          <  1:332881/79‑1 (MQ=255)                                                                     aTGGACGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCa  >  1:332915/1‑81 (MQ=255)                                                                          |                                                                            ATCCTGGTTCCCCAAAAGGGCAACCAACTGGAACGTTTTATCTGCATGGCATTATTTGTCTGGTCGGTATGGATGCCCCATTCATTACCCGAGAACAGCGCTGTTCGTTGGTCATCGCAAATGATCAATACCTTGCAGATCCAGTTTCA  >  NC_000913/1810622‑1810770

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