Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,447,891	G→A	*1260* (TAG→TAA)	tamB →	translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,447,891	0	G	A	100.0%	204.4 / NA	64	*1260* (TAG→TAA)	tamB	translocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/29);  total (35/29)

CGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCG  >  NC_000913/4447768‑4448023                                                                                                                            |                                                                                                                                     cGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGt                                                                                                                      >  1:903989/1‑139 (MQ=255)    gggTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTa                                                                                                                   <  2:721723/139‑1 (MQ=255)                ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                                       >  1:625743/1‑139 (MQ=255)                ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                                       >  1:378415/1‑139 (MQ=255)                   tataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                                                    >  1:380340/1‑139 (MQ=255)                         aaCACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                                                              <  2:903989/139‑1 (MQ=255)                             cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                                          >  1:97044/1‑139 (MQ=255)                             cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGACAGTTTACGCCACAAACAAGGc                                                                                          >  2:689113/1‑139 (MQ=255)                               cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  2:449047/139‑1 (MQ=255)                               cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  1:123540/1‑139 (MQ=255)                               cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  1:272998/1‑139 (MQ=255)                               cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  1:535676/1‑139 (MQ=255)                                acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                       <  2:185592/139‑1 (MQ=255)                                      cgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                 >  2:463611/1‑139 (MQ=255)                                        ttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  1:141767/130‑1 (MQ=255)                                        ttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  2:141767/1‑130 (MQ=255)                                           tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                            <  2:378415/139‑1 (MQ=255)                                                tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTa                                                                                                         <  2:173437/105‑1 (MQ=255)                                                tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTa                                                                                                         >  1:173437/1‑105 (MQ=255)                                                tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                                            <  2:858362/118‑1 (MQ=255)                                                tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                                            >  1:858362/1‑118 (MQ=255)                                                  aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAAt                                                                     >  2:575479/1‑139 (MQ=255)                                                         aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                                       <  1:729797/130‑1 (MQ=255)                                                         aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                                       >  2:729797/1‑130 (MQ=255)                                                              tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  1:615566/108‑1 (MQ=255)                                                              tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGATTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  2:615566/1‑108 (MQ=255)                                                                 cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                                                      <  2:535676/139‑1 (MQ=255)                                                                   ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                                                      >  1:237122/1‑137 (MQ=255)                                                                   ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                                                      <  2:237122/137‑1 (MQ=255)                                                                   ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGa                                                    >  1:580524/1‑139 (MQ=255)                                                                    gacgCCGTTTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                                                   <  2:802182/136‑1 (MQ=255)                                                                        ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                                                >  1:347080/1‑74 (MQ=255)                                                                        ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                                                <  2:347080/74‑1 (MQ=255)                                                                        ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                               <  2:580524/139‑1 (MQ=255)                                                                        ccGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                               <  2:380340/139‑1 (MQ=255)                                                                                 gtgtAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                      <  1:510103/139‑1 (MQ=255)                                                                                       aCCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAcc                                >  2:599507/1‑139 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                                          <  2:2581/61‑1 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                                          >  1:2581/1‑61 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                                                <  2:297563/71‑1 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                                                >  1:297563/1‑71 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  2:832141/1‑79 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  1:257944/79‑1 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  1:832141/79‑1 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                                            >  2:241498/1‑123 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                                            <  1:241498/123‑1 (MQ=255)                                                                                           ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCAACAAACAAGGCaa                                                                                        >  2:257944/1‑79 (MQ=255)                                                                                            gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGc                                                                                                                                >  2:493002/1‑38 (MQ=255)                                                                                            gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGc                                                                                                                                <  1:493002/38‑1 (MQ=255)                                                                                            gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACg                                                                                                                 >  2:667982/1‑53 (MQ=255)                                                                                            gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCac                              <  2:581117/136‑2 (MQ=255)                                                                                            gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTg                           >  2:586353/1‑139 (MQ=255)                                                                                            gCACTGGA‑TTGCTCTATAAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACTAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg                              >  1:581117/1‑136 (MQ=255)                                                                                                  ggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAg                      >  2:498862/4‑139 (MQ=255)                                                                                                             tCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                            >  2:737959/1‑74 (MQ=255)                                                                                                             tCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                            <  1:737959/74‑1 (MQ=255)                                                                                                             tCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                           >  2:35311/1‑107 (MQ=255)                                                                                                             tCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                           <  1:35311/107‑1 (MQ=255)                                                                                                                 ttCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  1:30769/58‑1 (MQ=255)                                                                                                                 ttCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  2:30769/1‑58 (MQ=255)                                                                                                                  tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                   <  1:329036/110‑1 (MQ=255)                                                                                                                  tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                   >  2:329036/1‑110 (MQ=255)                                                                                                                     aGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgc   >  1:461616/1‑139 (MQ=255)                                                                                                                      gTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcg  <  1:586353/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                     CGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCG  >  NC_000913/4447768‑4448023

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