Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,082,038	T→C	D42D (GAT→GAC)	loiP →	Phe‑Phe periplasmic metalloprotease, OM lipoprotein, low salt‑inducible, Era‑binding heat shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,082,038	0	T	C	100.0%	168.5 / NA	47	D42D (GAT→GAC)	loiP	Phe‑Phe periplasmic metalloprotease, OM lipoprotein, low salt‑inducible, Era‑binding heat shock protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/23);  total (24/23)

AAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTC  >  NC_000913/3081903‑3082171                                                                                                                                        |                                                                                                                                      aaGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGAcgcg                                                                                                                                    <  2:189930/139‑1 (MQ=255)           aTGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGAcgcg                                                                                                                                    >  2:939689/1‑129 (MQ=255)           aTGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGAcgcg                                                                                                                                    <  1:939689/129‑1 (MQ=255)              aaaaTTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAAccc                                                                                                                       >  1:133319/1‑139 (MQ=255)                     gcgcCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGa                                                                                                                <  1:337297/139‑1 (MQ=255)                        ccTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCa                                                                                                             >  2:717990/1‑139 (MQ=255)                              tGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATg                                                                                                       <  2:166667/139‑1 (MQ=255)                                      aTGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGgaga                                                                                               >  2:1171078/1‑139 (MQ=255)                                       tGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGAt                                                                                              >  1:1330480/1‑139 (MQ=255)                                          gCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCa                                                                                                                                  <  1:1213194/100‑1 (MQ=255)                                          gCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCa                                                                                                                                  >  2:1213194/1‑100 (MQ=255)                                           cGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGAc                                                                                          >  2:95816/1‑139 (MQ=255)                                                        ctgactgGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGAcgcg                                                                                                                                    <  1:1208333/84‑1 (MQ=255)                                                        ctgactgGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGAcgcg                                                                                                                                    >  2:1208333/1‑84 (MQ=255)                                                             tgGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTgcg                                                                        <  2:697810/139‑1 (MQ=255)                                                              ggTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTgcgc                                                                       >  2:711488/1‑139 (MQ=255)                                                                    cAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcca                                                                 <  2:1287969/139‑1 (MQ=255)                                                                              aCTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATa                                                       <  1:711488/139‑1 (MQ=255)                                                                              aCTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATa                                                       <  1:717990/139‑1 (MQ=255)                                                                                  cAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGa                                                                                                                             <  2:255839/65‑1 (MQ=255)                                                                                  cAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGa                                                                                                                             >  1:255839/1‑65 (MQ=255)                                                                                   aaCGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTa                                                  >  1:227115/1‑139 (MQ=255)                                                                                        aCTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTg                                                                                                                     <  2:240600/67‑1 (MQ=255)                                                                                        aCTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTg                                                                                                                     >  1:240600/1‑67 (MQ=255)                                                                                          tGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCt                                           >  2:52867/1‑139 (MQ=255)                                                                                           gCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTg                                          <  1:52867/139‑1 (MQ=255)                                                                                            ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCa                                                                                                             >  2:283837/1‑71 (MQ=255)                                                                                            ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCa                                                                                                             <  1:283837/71‑1 (MQ=255)                                                                                               tcatcaGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACaa                                      <  1:1171078/139‑1 (MQ=255)                                                                                                tatcaGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGc                                                                                                          <  2:1250349/69‑1 (MQ=255)                                                                                                tatcaGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGc                                                                                                          >  1:1250349/2‑70 (MQ=255)                                                                                                catcaGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAAccc                                                                                                                       >  1:1085557/1‑57 (MQ=255)                                                                                                catcaGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAAccc                                                                                                                       <  2:1085557/57‑1 (MQ=255)                                                                                                         cGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTc                                                                                                     >  2:934061/1‑66 (MQ=255)                                                                                                         cGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTc                                                                                                     <  1:934061/66‑1 (MQ=255)                                                                                                         cGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTAt                                  <  2:780564/133‑1 (MQ=255)                                                                                                         cGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTAt                                  >  1:780564/1‑133 (MQ=255)                                                                                                         cGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCaa                            <  1:95816/139‑1 (MQ=255)                                                                                                         cGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCaa                            >  2:430667/1‑139 (MQ=255)                                                                                                             aGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATgcg                        <  2:291703/139‑1 (MQ=255)                                                                                                                ttttCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTa                     <  2:133319/139‑1 (MQ=255)                                                                                                                 tttCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAg                    >  1:277121/1‑139 (MQ=255)                                                                                                                  ttCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCa                                                                                                             >  1:275211/1‑49 (MQ=255)                                                                                                                  ttCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCa                                                                                                             <  2:275211/49‑1 (MQ=255)                                                                                                                  ttCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAgg                   <  2:417921/139‑1 (MQ=255)                                                                                                                     aGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGc                                                                                                          >  2:659318/1‑49 (MQ=255)                                                                                                                     aGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGc                                                                                                          <  1:659318/49‑1 (MQ=255)                                                                                                                         ttACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAat            >  1:783112/1‑139 (MQ=255)                                                                                                                          tACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAata           >  2:399671/1‑139 (MQ=255)                                                                                                                          tACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAata           >  2:398829/1‑139 (MQ=255)                                                                                                                                   aGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTc  <  2:227115/139‑1 (MQ=255)                                                                                                                                   aGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTc  >  2:1465984/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                      AAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTATTGCCAATGCGCTAGGCAACAATATCAACGGTC  >  NC_000913/3081903‑3082171

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