Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,341,617	A→T	D4E (GAT→GAA)	melR ←	melibiose operon transcriptional regulator; autoregulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,341,617	0	A	T	100.0%	118.9 / NA	31	D4E (GAT→GAA)	melR	melibiose operon transcriptional regulator; autoregulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (16/15);  total (16/15)

GAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTC  >  NC_000913/4341470‑4341760                                                                                                                                                    |                                                                                                                                                gAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCt                                                                                                                                              >  1:536997/1‑151 (MQ=255)       cacaTCCCCTCCGCCCTGCCCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAACCTCAATCTCCATCCGCTGGTATTCCGAATCCAGCGCCAGCGGGTTGCGGGTCTG‑TTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTc                                                                                                                                        <  2:532962/150‑1 (MQ=255)        acaTTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCa                                                                                                                                       <  2:198868/151‑1 (MQ=255)         caTTCACTTCGACCTGACCATGCCAGCGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCat                                                                                                                                      >  1:370609/1‑151 (MQ=255)          aTTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCatg                                                                                                                                     <  2:680147/151‑1 (MQ=255)                                 aGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAatat                                                                                                              <  1:214259/151‑1 (MQ=255)                                  gTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATc                                                                                                             <  1:499877/151‑1 (MQ=255)                                      cTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGa                                                                                                          >  1:509632/1‑150 (MQ=255)                                         gtggGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTa                                                                                                      <  1:434147/151‑1 (MQ=255)                                                         gTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTgaggg                                                                                       >  1:262648/1‑150 (MQ=255)                                                         gTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTgaggg                                                                                       >  2:61476/1‑150 (MQ=255)                                                              cgAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATg                                                                                 >  1:346415/1‑151 (MQ=255)                                                                    tCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCg                                                                           >  1:219280/1‑151 (MQ=255)                                                                     cAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGa                                                                          <  2:142196/151‑1 (MQ=255)                                                                               ccGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGAACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCAc                                                                >  1:286073/1‑151 (MQ=255)                                                                                 gCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACgg                                                              >  2:54104/1‑151 (MQ=255)                                                                                    ggTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGttt                                                           <  1:419667/151‑1 (MQ=255)                                                                                      tATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGtttt                                                          >  2:539967/1‑150 (MQ=255)                                                                                             aaTACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCctct                                                  <  2:260010/151‑1 (MQ=255)                                                                                                 cagcgacagcgGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTccc                                              <  1:189849/151‑1 (MQ=255)                                                                                                    cgacagcgGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTccc                                              >  1:329128/1‑148 (MQ=255)                                                                                                       cagcgGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAggg                                        <  1:121019/151‑1 (MQ=255)                                                                                                                gCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAgggg                                >  2:642676/1‑150 (MQ=255)                                                                                                                gCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGa                               <  1:287109/151‑1 (MQ=255)                                                                                                                  gggTCTGTTTTTCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGAct                             <  1:630002/151‑1 (MQ=255)                                                                                                                          ttttCGTCGCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGt                     >  1:44151/1‑151 (MQ=255)                                                                                                                               gtcgCTGCTGCACATAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCAt                <  2:90591/151‑1 (MQ=255)                                                                                                                                         cacaTAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCTAGGGGACTCTCCGAGTATCATGAGGCCGaa       >  1:201421/1‑150 (MQ=255)                                                                                                                                          acaTAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGaaaa     <  1:624148/151‑1 (MQ=255)                                                                                                                                          acaTAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGaaa      >  2:478712/1‑150 (MQ=255)                                                                                                                                             tAAACGTTTCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAActc  >  2:366989/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                GAAAGGCACATTCACTTCGACCTGACCATGCCAGTGGCTGGTGGGCATGATATGCGGTGCGCGAAACTCAATCTCCATCCGCTGGTATTCCGAATACAGCGACAGCGGGCTGCGGGTCTGTTTTTCGTCGCTGCTGCACATAAACGTATCTGTATTCATGGATGGCTCTCTTTCCTGGAATATCAGAATTATGGCAGGAGTGAGGGAGGATGACTGCGAGTGGGAGCACGGTTTTCACCCTCTTCCCAGAGGGGCGAGGGGACTCTCCGAGTATCATGAGGCCGAAAACTC  >  NC_000913/4341470‑4341760

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