Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,691,784	A→G	S356S (AGT→AGC)	uidC ←	putative outer membrane porin for beta‑glucuronides porin protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,691,784	0	A	G	100.0%	187.1 / NA	50	S356S (AGT→AGC)	uidC	putative outer membrane porin for beta‑glucuronides porin protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/25);  total (25/25)

CGCCCAGCGTACTGCCATCGCGCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTACTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACTTAATGAAAAGGGATAGCCAATGTCGGTAT  >  NC_000913/1691645‑1691929                                                                                                                                            |                                                                                                                                                  cGCCCAGCGTACTGCCATCGCGCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATGCCCGTGCCTTCAATCCGCGCATCTCGTTTGCTGGATTCAt                                                                                                                                         <  2:500193/150‑1 (MQ=255)  gCCCAGCGTACTGCCATCGCGCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTc                                                                                                                                           >  1:821401/1‑147 (MQ=255)   cccAGCGTACTGCCATCGCGCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGc                                                                                                                                      <  1:84010/151‑1 (MQ=255)      aGCGTACTGCCATCGCGCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCcgc                                                                                                                                   >  2:456403/1‑151 (MQ=255)        cGTACTGCCATCGCGCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGTcgcg                                                                                                                                  >  1:269329/1‑150 (MQ=255)              gCCATCGCGCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTcaca                                                                                                                           <  1:521414/151‑1 (MQ=255)                 aTCGCGCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATc                                                                                                                        >  2:959785/1‑151 (MQ=255)                   cgcgCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCgg                                                                                                                      >  2:585180/1‑151 (MQ=255)                        ttttCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCa                                                                                                                  >  2:67672/1‑150 (MQ=255)                          ttCCCGCCCTTTATCGGCAGCAAGAAAGAAGTTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTATGCTGTCTTCTGATTGACAAATCACAATCGGTGCaaa                                                                                                                >  1:605662/1‑150 (MQ=255)                           tCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAAt                                                                                                               <  1:736323/150‑1 (MQ=255)                              cGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTc                                                                                                            >  1:231455/1‑150 (MQ=255)                                 tccgatttggCAGCAAGAAAGAAATTCATCCCTTGTAACGGCTCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAg                                                                                                          >  2:732548/9‑149 (MQ=255)                                    ttATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTa                                                                                                      >  2:452007/1‑150 (MQ=255)                                     tATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTaa                                                                                                     <  1:694920/150‑1 (MQ=255)                                      aTCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAAt                                                                                                   >  1:776043/1‑151 (MQ=255)                                      aTCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCCGGCTTCATAGTCGAAAAACACAATCGGTGCAAATGTCAGCGTAAAt                                                                                                   >  1:59158/1‑151 (MQ=255)                                      aTCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCCTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGCGCAAATGTCAGCGTaaa                                                                                                    >  2:978732/1‑150 (MQ=255)                                         ggCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTgc                                                                                                <  1:330146/151‑1 (MQ=255)                                           cagcaAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTgcgg                                                                                              >  1:180792/1‑151 (MQ=255)                                             gcaAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTgcggcg                                                                                            >  1:289956/1‑151 (MQ=255)                                               aagaaagaAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGtt                                                                                          <  2:118385/151‑1 (MQ=255)                                                   aagaaaTTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAAc                                                                                      >  2:401445/1‑151 (MQ=255)                                                    agaaaTTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACg                                                                                     <  2:703359/151‑1 (MQ=255)                                                     gaaaTTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGa                                                                                    <  2:93450/151‑1 (MQ=255)                                                      aaaTTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGAt                                                                                   <  1:394596/151‑1 (MQ=255)                                                       aaTTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGAt                                                                                   >  2:98775/1‑150 (MQ=255)                                                        aTTCATCCCTTGTGACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAg                                                                                 <  2:158604/151‑1 (MQ=255)                                                           cATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGtt                                                                               <  2:151320/150‑1 (MQ=255)                                                           cATCCCTTGTAACGGCCCTTCGCCGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGtt                                                                               <  1:589344/150‑1 (MQ=255)                                                             tCCCTTGCAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATAGCGGCGTTAAACGATAGTTGAc                                                                            >  2:298405/1‑151 (MQ=255)                                                              cccTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCCGGCTTCATCGCCGAGAGACACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACg                                                                           >  1:545624/1‑151 (MQ=255)                                                                 ttGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGccc                                                                        >  2:9177/1‑151 (MQ=255)                                                                   gTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCa                                                                       <  1:61094/150‑1 (MQ=255)                                                                     aaCGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAAtt                                                                    >  2:833312/1‑151 (MQ=255)                                                                               cGCTGACTCGATAGGTCATCCCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCTGCGTTAAACGATAGTTGCCGTCCAATTGCCAGGACt                                                           <  2:613137/150‑1 (MQ=255)                                                                               cGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTg                                                          <  1:111561/151‑1 (MQ=255)                                                                               cGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTg                                                          >  2:934242/1‑151 (MQ=255)                                                                                gCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGc                                                         <  1:824345/151‑1 (MQ=255)                                                                                   gACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCAt                                                       <  1:158551/150‑1 (MQ=255)                                                                                          tAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGGTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCAc                                               >  2:648989/1‑151 (MQ=255)                                                                                                    caccTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATc                                     <  2:321953/151‑1 (MQ=255)                                                                                                           gATACCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACtt                              >  2:398067/1‑151 (MQ=255)                                                                                                              aCCCGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACTTaa                            <  2:977521/150‑1 (MQ=255)                                                                                                                 cGTGCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACTTAATGaa                        <  2:289923/151‑1 (MQ=255)                                                                                                                    gCCTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACTTAATGAAAAg                     <  1:70847/151‑1 (MQ=255)                                                                                                                     ccTTCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACTTAATGAAAAgg                    >  1:785521/1‑151 (MQ=255)                                                                                                                        tCAATCCGCACATCTCGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATCGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACTTTATGAAAAGGGa                  <  2:49252/150‑1 (MQ=255)                                                                                                                                       cGTTTGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACTTAATGAAAAGGGATAGCCAATGTCGGTAt  <  1:456455/151‑1 (MQ=255)                                                                                                                                         tttGCTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACTTAATGAAAAGGGATAGCCAATGTCGGTAt  <  2:384065/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                  CGCCCAGCGTACTGCCATCGCGCTTTTCCCGCCCTTTATCGGCAGCAAGAAAGAAATTCATCCCTTGTAACGGCCCTTCGCTGACCCGATAGTTCATACCACCTAAGATACCCGTGCCTTCAATCCGCACATCTCGTTTACTGGATTCATAGCCGCGAGTCACAATCGGTGCAAATGTCAGCGTAAATTGCGGCGTTAAACGATAGTTGACGCCCAATTGCCAGGACTGCATTCCTTCACGGTTACGATCAAGACTTAATGAAAAGGGATAGCCAATGTCGGTAT  >  NC_000913/1691645‑1691929

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