Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,081,989	G→A	G26E (GGA→GAA)	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,081,989	0	G	A	100.0%	210.1 / NA	64	G26E (GGA→GAA)	loiP	Phe‑Phe periplasmic metalloprotease, OM lipoprotein, low salt‑inducible, Era‑binding heat shock protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/34);  total (30/34)

ATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTA  >  NC_000913/3081855‑3082123                                                                                                                                       |                                                                                                                                       aTTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTg                                                                                                                                    <  2:305530/139‑1 (MQ=255)   tCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGct                                                                                                                                  <  1:1077048/139‑1 (MQ=255)   tCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGct                                                                                                                                  <  2:641095/139‑1 (MQ=255)   tCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGct                                                                                                                                  <  1:552774/139‑1 (MQ=255)      ccTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCtctc                                                                                                                               >  2:1127903/1‑139 (MQ=255)       cTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCtca                                                                                                                              <  2:787384/139‑1 (MQ=255)            aaCGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAgg                                                                                                                         >  1:296258/1‑139 (MQ=255)                    aTTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAg                                                                                                                 >  1:1042152/1‑139 (MQ=255)                       tttAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCtt                                                                                                              >  2:156080/1‑139 (MQ=255)                         tAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCtttt                                                                                                            <  2:296258/139‑1 (MQ=255)                         tAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCtttt                                                                                                            <  1:156080/139‑1 (MQ=255)                           gAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCa                                                                                                          <  1:866539/139‑1 (MQ=255)                               aatCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCt                                                                                                      <  2:1019081/139‑1 (MQ=255)                                atCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCtt                                                                                                     <  2:123376/139‑1 (MQ=255)                                     gACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGt                                                                                                <  2:722666/139‑1 (MQ=255)                                          tgtgCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAg                                                                                           >  2:1083419/1‑139 (MQ=255)                                          tgtgCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAg                                                                                           >  1:434759/1‑139 (MQ=255)                                           gtgCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCtcatc                                                                                                                            >  1:467033/1‑105 (MQ=255)                                           gtgCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCtcatc                                                                                                                            <  2:467033/105‑1 (MQ=255)                                           gtgCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAgg                                                                                                        >  1:500108/1‑125 (MQ=255)                                           gtgCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAgg                                                                                                        <  2:500108/125‑1 (MQ=255)                                           gtacGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCtcatc                                                                                                                            >  2:811974/4‑105 (MQ=255)                                           gtacGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCtcatc                                                                                                                            <  1:811974/102‑1 (MQ=255)                                            tgCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTg                                                                                                                                    >  2:825070/1‑96 (MQ=255)                                            tgCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTg                                                                                                                                    <  1:825070/96‑1 (MQ=255)                                            tgCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTg                                                                                         >  2:214923/1‑139 (MQ=255)                                               ggAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGacg                                                                                      <  1:1127436/139‑3 (MQ=255)                                                 aaGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGAcgcg                                                                                    >  2:415712/1‑139 (MQ=255)                                                 aaGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGAcgcg                                                                                    <  1:415712/139‑1 (MQ=255)                                                         aCATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGaa                                                                            <  2:1042152/139‑1 (MQ=255)                                                           aTGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACT‑GTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAg                                                                                                                          <  1:844377/90‑1 (MQ=255)                                                           aTGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACT‑GTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAg                                                                                                                          >  2:844377/1‑90 (MQ=255)                                                             gAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAAcc                                                                        >  2:1132371/1‑139 (MQ=255)                                                                           tATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAgg                                                           <  2:457413/138‑1 (MQ=255)                                                                           tATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAgg                                                           >  1:457413/1‑138 (MQ=255)                                                                                gTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTc                                                     <  1:1127903/139‑1 (MQ=255)                                                                                   gCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGacg                                                                                      <  1:822823/103‑3 (MQ=255)                                                                                   gCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGacg                                                                                      >  2:822823/1‑101 (MQ=255)                                                                                     aaTGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGacg                                                                                      >  2:37400/1‑99 (MQ=255)                                                                                     aaTGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGacg                                                                                      <  1:37400/101‑3 (MQ=255)                                                                                     aaTGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGgag                                                >  1:797977/1‑139 (MQ=255)                                                                                       tGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCg                                                                                                                     >  2:95321/1‑68 (MQ=255)                                                                                       tGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCg                                                                                                                     <  1:95321/68‑1 (MQ=255)                                                                                         aGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGc                                                          <  1:1021908/125‑1 (MQ=255)                                                                                         aGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGc                                                          >  2:1021908/1‑125 (MQ=255)                                                                                           cGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTAACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGAc                                          >  2:687916/1‑139 (MQ=255)                                                                                              ggCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCtcatc                                                                                                                            >  1:492660/1‑54 (MQ=255)                                                                                              ggCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCtcatc                                                                                                                            <  2:492660/54‑1 (MQ=255)                                                                                                   cGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTg                                                                     <  2:877285/104‑1 (MQ=255)                                                                                                   cGGTACTGACTGGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTg                                                                     >  1:877285/1‑104 (MQ=255)                                                                                                        ctgactgGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCa                                                                                                          >  1:792449/1‑62 (MQ=255)                                                                                                        ctgactgGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCa                                                                                                          <  2:792449/62‑1 (MQ=255)                                                                                                          gactgGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAAcc                                                                        <  2:243677/94‑1 (MQ=255)                                                                                                          gactgGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAAcc                                                                        >  1:243677/1‑94 (MQ=255)                                                                                                             tgGTTGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTgcg                        <  1:1083419/139‑1 (MQ=255)                                                                                                                ttGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAgg                                                           >  2:919633/1‑101 (MQ=255)                                                                                                                ttGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAgg                                                           <  1:919633/101‑1 (MQ=255)                                                                                                                ttGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATg                                             <  2:153083/115‑1 (MQ=255)                                                                                                                ttGCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATg                                             >  1:153083/1‑115 (MQ=255)                                                                                                                  gCCAGAATATGGACTCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagc                   <  1:263253/139‑1 (MQ=255)                                                                                                                                tCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAg                                                   <  1:737749/93‑1 (MQ=255)                                                                                                                                tCCAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAg                                                   >  2:737749/1‑93 (MQ=255)                                                                                                                                 ccAACGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGc    >  1:8810/1‑139 (MQ=255)                                                                                                                                   aaCGAACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGACGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTa  >  2:207589/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                       ATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCGGAAGAGAAAACATGAAAATTCGCGCCTTATTGGTAGCAATGAGCGTGGCAACGGTACTGACTGGTTGCCAGAATATGGACTCCAACGGACTGCTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTA  >  NC_000913/3081855‑3082123

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: ‑