Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	181,996	T→C	G371G (GGT→GGC)	degP →	periplasmic serine endoprotease DegP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	181,996	0	T	C	100.0%	128.2 / NA	38	G371G (GGT→GGC)	degP	periplasmic serine endoprotease DegP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/20);  total (18/20)

GCAAAAGCGGGCATTAAAGCGGGTGATGTGATCACCTCACTGAACGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGTAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTG  >  NC_000913/181859‑182122                                                                                                                                          |                                                                                                                               gCAAAAGCGGGCATTAAAGCGGGTGATGTGATCACCTCACTGAACGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCaa                                                                                                                              >  1:450050/1‑140 (MQ=255)  cAAAAGCGGGCATTAAAGCGGGTGATGTGATCACCTCACTGAACGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACGATGCCGGTCGGCGGCAGACTGACCCAGGGCTTACTGCGCGACGGCAAg                                                                                                                             >  2:226451/1‑140 (MQ=255)                            gtgatCACCTCACTGAACGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTgca                                                                                                   <  2:450050/140‑1 (MQ=255)                                      cACTGAACGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGagcagc                                                                                         >  2:710667/1‑140 (MQ=255)                                         tGAACGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGTGACGGCAAACAGTTTAGTGTGAACCGGGAACTGCAGCAGAGCAGCCAg                                                                                      >  2:341172/1‑140 (MQ=255)                                         tGAACGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGTGACGGCAAACAGGTTAACGTGAACCTGGAAATGCAGCAGAGCAGCCAg                                                                                      >  2:255131/1‑140 (MQ=255)                                           aaCGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGaa                                                                                    >  2:734491/1‑140 (MQ=255)                                           aaCGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGaa                                                                                    >  1:384934/1‑140 (MQ=255)                                                  aGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGtt                                                                             <  2:38546/140‑1 (MQ=255)                                                     cGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCGGAGCGGCCAGAATCAGGTTGAt                                                                          >  2:230666/1‑140 (MQ=255)                                                       aTCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGCTAACGTGAACCTAGGACTGCAGCAGGGCAGGCAGAATCAGGTTGATTc                                                                        >  2:323779/1‑140 (MQ=255)                                                        tCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTGAACGTGAATCTGTAACCGCAGCAGAGCAGCCAGAATCAGGTTGATTcc                                                                       >  2:216172/1‑140 (MQ=255)                                                         cagcagCTTTGCCGCTCTGCGTGCTTAGGTGGGTACTTTGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCa                                                                      <  1:255131/140‑1 (MQ=255)                                                         cagcagCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCa                                                                      <  2:384934/140‑1 (MQ=255)                                                                      gCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCaa                                                         <  1:734491/140‑1 (MQ=255)                                                                            cGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCAt                                                   >  2:238095/1‑140 (MQ=255)                                                                            cGTGCTCAGGAGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCAt                                                   >  1:383321/1‑140 (MQ=255)                                                                              tGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTg                                                 >  1:501163/1‑140 (MQ=255)                                                                               gCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGa                                                <  2:671551/140‑1 (MQ=255)                                                                                      tGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCt                                         >  1:652029/1‑140 (MQ=255)                                                                                        ggTACTATGCCGGTAGGCAGCAAACTGACCCTGGGC‑TACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTgag                                      <  2:47603/140‑1 (MQ=255)                                                                                        ggTACTATGCCGGTAGGCAGCAAACTGACCCTGGGC‑TACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTgag                                      <  2:365403/140‑1 (MQ=255)                                                                                               tGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGc                                >  2:518227/1‑140 (MQ=255)                                                                                               tGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGc                                <  1:518227/140‑1 (MQ=255)                                                                                                  cGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGcaac                             <  2:383321/140‑1 (MQ=255)                                                                                                    gTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGcaacaa                           >  2:15759/1‑140 (MQ=255)                                                                                                    gTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGcaacaa                           <  1:401386/140‑1 (MQ=255)                                                                                                      aGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAg                         >  1:67549/1‑140 (MQ=255)                                                                                                         cagcaAACTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCa                      <  1:563528/140‑1 (MQ=255)                                                                                                              aaCTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGAt                 <  1:238095/140‑1 (MQ=255)                                                                                                               aCTGACCCTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATc                >  2:515682/1‑140 (MQ=255)                                                                                                                  gACCCTGGGCTTACTTCGCGACCGCAAGCCGGTTAACATGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAgg             <  1:226451/140‑1 (MQ=255)                                                                                                                     ccTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCg          <  1:515682/140‑1 (MQ=255)                                                                                                                     ccTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCg          <  2:501163/140‑1 (MQ=255)                                                                                                                     ccTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCg          <  1:710667/140‑1 (MQ=255)                                                                                                                     ccTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCg          <  2:699463/140‑1 (MQ=255)                                                                                                                     ccTGGGCTTACTGCGCGACGGCAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCg          <  1:15759/140‑1 (MQ=255)                                                                                                                             tACTGCGCGACGGCACGCAGGTCAACGTGCACCTGGACCTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTg  <  1:216172/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                               GCAAAAGCGGGCATTAAAGCGGGTGATGTGATCACCTCACTGAACGGTAAGCCGATCAGCAGCTTTGCCGCACTGCGTGCTCAGGTGGGTACTATGCCGGTAGGCAGCAAACTGACCCTGGGCTTACTGCGCGACGGTAAGCAGGTTAACGTGAACCTGGAACTGCAGCAGAGCAGCCAGAATCAGGTTGATTCCAGCTCCATCTTCAACGGCATTGAAGGCGCTGAGATGAGCAACAAAGGCAAAGATCAGGGCGTGGTAGTG  >  NC_000913/181859‑182122

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