Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,106,320	G→A	*167* (TAG→TAA)	cpxP →	inhibitor of the cpx response, periplasmic adaptor protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,106,320	0	G	A	100.0%	112.6 / NA	36	*167* (TAG→TAA)	cpxP	inhibitor of the cpx response, periplasmic adaptor protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/18);  total (18/18)

ACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGA  >  NC_000913/4106230‑4106454                                                                                           |                                                                                                                                       acaacGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCcaca                                                                                        <  1:16415/139‑1 (MQ=255)      gAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCtg                                                                                   >  2:728999/1‑139 (MQ=255)             gCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGaca                                                                                                             <  1:136602/106‑1 (MQ=255)             gCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGaca                                                                                                             >  2:136602/1‑106 (MQ=255)             gCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGaca                                                                                                             >  2:411453/1‑106 (MQ=255)             gCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGaca                                                                                                             <  1:411453/106‑1 (MQ=255)               aGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttt                                                                          <  1:600564/139‑1 (MQ=255)                 ttGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                         <  1:2337/90‑1 (MQ=255)                 ttGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                         >  2:2337/1‑90 (MQ=255)                 ttGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGttttt                                                                        >  2:786096/1‑139 (MQ=255)                      tgacgtgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTc                                                                                                  >  1:223180/1‑108 (MQ=255)                      tgacgtgacgCAATGGCAAAAAAGTTCATCGTTCAAGCTATTGAGTAGTAGCATCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTc                                                                                                  <  2:223180/108‑1 (MQ=255)                           tgacgCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCa                                                              >  2:731233/1‑139 (MQ=255)                                  aTGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGcc                                                       <  1:786096/139‑1 (MQ=255)                                    ggCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCctg                                                     <  1:673773/139‑1 (MQ=255)                                        aaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                         <  2:845732/67‑1 (MQ=255)                                        aaaaaGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                         >  1:845732/1‑67 (MQ=255)                                                             aTTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCa                                         <  1:527445/126‑1 (MQ=255)                                                             aTTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCa                                         >  2:527445/1‑126 (MQ=255)                                                             aTTGAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGAttt                            >  1:920850/1‑139 (MQ=255)                                                                gAGTAGTAGCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCt                         <  1:728999/139‑1 (MQ=255)                                                                      tagCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                        >  1:847810/1‑118 (MQ=255)                                                                      tagCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCAt                                        <  2:847810/118‑1 (MQ=255)                                                                        gCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                         >  2:249518/1‑35 (MQ=255)                                                                        gCAACTCACGTTCCCAGTAATAAACCCTGTTTTcc                                                                                                                         <  1:249518/35‑1 (MQ=255)                                                                        gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCt                         >  1:691860/1‑131 (MQ=255)                                                                        gCAACTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCt                         <  2:691860/131‑1 (MQ=255)                                                                           aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTc                                                                                                  <  2:396246/55‑1 (MQ=255)                                                                           aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTc                                                                                                  >  1:396246/1‑55 (MQ=255)                                                                           aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCcaca                                                                                        <  1:779342/65‑1 (MQ=255)                                                                           aCTCACGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCcaca                                                                                        >  2:779342/1‑65 (MQ=255)                                                                                cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGacac                                                                                                            >  2:642713/1‑40 (MQ=39)                                                                                cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGacac                                                                                                            <  1:642713/40‑1 (MQ=255)                                                                                cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGc                          >  2:112405/1‑122 (MQ=255)                                                                                cGTTCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGc                          <  1:112405/122‑1 (MQ=255)                                                                                  ttCCCAGTAATAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCa       >  2:850992/1‑139 (MQ=255)                                                                                       agtaatAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCg                     <  2:42118/120‑1 (MQ=255)                                                                                       agtaatAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCg                     >  1:42118/1‑120 (MQ=255)                                                                                       agtaatAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGa  >  1:667241/1‑139 (MQ=255)                                                                                           |                                                                                                                                       ACAACGAATGGAGCAGTTGCGTGACGTGACGCAATGGCAAAAAAGTTCATCGTTGAAGCTATTGAGTAGTAGCAACTCACGTTCCCAGTAGTAAACCCTGTTTTCCTTGCCATAGACACCATCCCTGTCTTCCCCCACATGCTGTGGGGGTTTTTTTTATCCTCAATTTGCCTGCTGCTTAATGCATTGCAGATGATTTGCTTCCGTTATACTAGCGTCAGTTGA  >  NC_000913/4106230‑4106454

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