Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	933,022	A→G	D143G (GAC→GGC)	lrp →	leucine‑responsive global transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	933,022	0	A	G	100.0%	111.5 / NA	30	D143G (GAC→GGC)	lrp	leucine‑responsive global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/20);  total (10/20)

ACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTA  >  NC_000913/932885‑933169                                                                                                                                          |                                                                                                                                                    aCAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATAc                                                                                                                                        <  1:340152/151‑1 (MQ=255)   aaaaaCTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACGCGGACATACgt                                                                                                                                      <  1:627744/151‑1 (MQ=255)   aaaaaCTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACgt                                                                                                                                      <  1:498726/151‑1 (MQ=255)     aaaCTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACgttg                                                                                                                                    <  1:530434/151‑1 (MQ=255)     aaaCTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACgttg                                                                                                                                    <  1:382351/151‑1 (MQ=255)               aaaTTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGgaag                                                                                                                           <  1:365700/150‑1 (MQ=255)                 aTTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGaagaag                                                                                                                        >  1:118284/1‑151 (MQ=255)                                          gTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGtta                                                                                                <  1:6274/150‑1 (MQ=255)                                                          gTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAaca                                                                                >  1:642066/1‑150 (MQ=255)                                                           ttGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAacac                                                                               <  2:542514/150‑1 (MQ=255)                                                                                          acagTAAGTTGCTGGGGGAAACCCTGCTGCGCCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGACGTCAAGGAGAGTAATCGTCTGGTTATTACGAGGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTa                                               >  1:559430/4‑151 (MQ=255)                                                                                          aCCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTa                                               >  1:284669/1‑151 (MQ=255)                                                                                                        ggggAAACCCTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAAt                                  <  1:106303/150‑1 (MQ=255)                                                                                                               cccTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAg                          >  1:17984/1‑151 (MQ=255)                                                                                                               cccTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAg                          >  2:377411/1‑151 (MQ=255)                                                                                                                ccTGCTGCGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGc                         <  1:418049/151‑1 (MQ=255)                                                                                                                    ctgcGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAAc                      <  2:42356/150‑1 (MQ=255)                                                                                                                    ctgcGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAAc                      <  1:297799/150‑1 (MQ=255)                                                                                                                       cGTCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCATGTTAATCCATACAGCAACAGTa                  >  1:80091/1‑151 (MQ=255)                                                                                                                         tCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACGCGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTAc                 <  2:232462/150‑1 (MQ=255)                                                                                                                         tCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTAc                 <  1:352867/150‑1 (MQ=255)                                                                                                                         tCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTAc                 <  1:273701/150‑1 (MQ=255)                                                                                                                         tCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTAc                 <  1:684100/150‑1 (MQ=255)                                                                                                                         tCTGCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTAc                 <  1:689149/150‑1 (MQ=255)                                                                                                                            gCCTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTggg             >  1:452704/1‑151 (MQ=255)                                                                                                                             ccTGGCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTgggg            <  1:379149/151‑1 (MQ=255)                                                                                                                                ggCGTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTaa         >  1:134658/1‑151 (MQ=255)                                                                                                                                   gTCAATGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAAc        >  2:437765/1‑149 (MQ=255)                                                                                                                                    tCAATGGCACACGGACATACGTTGTTTTGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCt      <  2:599064/150‑1 (MQ=255)                                                                                                                                       aTGGCACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTa  <  1:377410/151‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                    ACAAAAACTTGAAGAAATTCAGGAGTGTCATTTAGTATCCGGTGATTTCGACTACCTGTTGAAAACACGCGTGCCGGATATGTCAGCCTACCGTAAGTTGCTGGGGGAAACCCTGCTGCGTCTGCCTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTA  >  NC_000913/932885‑933169

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