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%	230.4 / NA	71	*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/36);  total (35/36)

ACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTG  >  NC_000913/4447764‑4448007                                                                                                                                |                                                                                                                     aCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                          <  1:805006/139‑1 (MQ=255)          gTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACg                                                                                                 <  2:1234890/139‑1 (MQ=255)                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                          >  2:895289/1‑139 (MQ=255)                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                          >  1:437310/1‑139 (MQ=255)                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGACACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGt                                                                                            >  2:1183722/1‑121 (MQ=255)                                 cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGACACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGt                                                                                            <  1:1183722/121‑1 (MQ=255)                                    acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                       >  2:414999/1‑139 (MQ=255)                                              aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAGTATTTGTCTACGGCAGTTTa                                                                                         <  2:887509/111‑1 (MQ=255)                                              aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTa                                                                                         >  1:887509/1‑111 (MQ=255)                                                cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGTAGTTTACGCCACAAACAAGGCAACAGTCACTgg                                                              >  2:1409001/1‑136 (MQ=255)                                                cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGTAGTTTACGCCACAAACAAGGCAACAGTCACTgg                                                              <  1:1409001/136‑1 (MQ=255)                                                   cTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGAcc                                                        <  2:437310/139‑1 (MQ=255)                                                    tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                       >  1:1037219/1‑139 (MQ=255)                                                    tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                       <  1:414999/139‑1 (MQ=255)                                                            aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACg                                                                                                 <  1:1400288/89‑1 (MQ=255)                                                            aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACg                                                                                                 >  2:1400288/1‑89 (MQ=255)                                                             aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGt                                              >  2:1275888/1‑139 (MQ=255)                                                                  tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                       <  2:434609/93‑1 (MQ=255)                                                                  tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                       >  1:434609/1‑93 (MQ=255)                                                                  tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                    >  2:635581/1‑112 (MQ=255)                                                                  tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                    <  1:635581/112‑1 (MQ=255)                                                                        gAAGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                                   <  2:416863/139‑1 (MQ=255)                                                                         aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                        >  1:1309909/1‑101 (MQ=255)                                                                         aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                        <  2:1309909/101‑1 (MQ=255)                                                                         aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                    >  2:345472/1‑105 (MQ=255)                                                                         aaGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                    <  1:345472/105‑1 (MQ=255)                                                                          aGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                                <  1:493531/76‑1 (MQ=255)                                                                          aGCCGTGTCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                                >  2:493531/1‑76 (MQ=255)                                                                          aGCCGTGTATGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGAAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAttt                                 <  1:895289/139‑1 (MQ=255)                                                                              gtgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGagc                                                        >  1:843568/1‑110 (MQ=255)                                                                              gtgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGagc                                                        <  2:843568/112‑3 (MQ=255)                                                                               tgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGTCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                            >  2:1263399/1‑139 (MQ=255)                                                                               tgtCTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                            <  1:1275888/139‑1 (MQ=255)                                                                                  cTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                        >  2:1300789/1‑92 (MQ=255)                                                                                  cTGGTGTAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                        <  1:1300789/92‑1 (MQ=255)                                                                                       gtAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                          >  1:792957/1‑69 (MQ=255)                                                                                       gtAGACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                          <  2:792957/69‑1 (MQ=255)                                                                                          gACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                                      <  2:114463/70‑1 (MQ=255)                                                                                          gACCAGGCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                                      >  1:114463/1‑70 (MQ=255)                                                                                               ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGaa                                                                                                             <  1:665794/42‑1 (MQ=255)                                                                                               ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGaa                                                                                                             >  2:665794/1‑42 (MQ=255)                                                                                               ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                                    <  1:667754/67‑1 (MQ=255)                                                                                               ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCca                                                                                    >  2:667754/1‑67 (MQ=255)                                                                                               ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                                >  1:751318/1‑71 (MQ=255)                                                                                               ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                                <  2:751318/71‑1 (MQ=255)                                                                                               ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                            <  1:358564/75‑1 (MQ=255)                                                                                               ggCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                                            >  2:358564/1‑75 (MQ=255)                                                                                                gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                        <  2:51129/78‑1 (MQ=255)                                                                                                gCACTGGA‑TTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                        >  1:51129/1‑78 (MQ=255)                                                                                                     tggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAttt                                 >  2:29645/5‑113 (MQ=255)                                                                                                     tggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAttt                                 <  1:29645/109‑1 (MQ=255)                                                                                                     tggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGtata       >  1:880182/5‑139 (MQ=255)                                                                                                          tGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTg  <  2:880182/139‑1 (MQ=255)                                                                                                            ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                            >  1:296340/1‑79 (MQ=255)                                                                                                            ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                            <  2:296340/79‑1 (MQ=255)                                                                                                            ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                       >  2:434892/1‑84 (MQ=255)                                                                                                            ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                       <  1:434892/84‑1 (MQ=255)                                                                                                               tATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                   >  1:330501/1‑117 (MQ=255)                                                                                                               tATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                   <  2:330501/117‑1 (MQ=255)                                                                                                                     ttCGAGTTTTAACAATGCGAATACTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                               >  1:310207/1‑99 (MQ=255)                                                                                                                     ttCGAGTTTTAACAATGCGAATACTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                               <  2:310207/99‑1 (MQ=255)                                                                                                                      tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                                     >  2:1506926/1‑44 (MQ=255)                                                                                                                      tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                                     <  1:1506926/44‑1 (MQ=255)                                                                                                                      tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                        <  2:522660/57‑1 (MQ=255)                                                                                                                      tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                        >  1:522660/1‑57 (MQ=255)                                                                                                                      tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                  <  1:891915/79‑1 (MQ=255)                                                                                                                      tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                  >  2:891915/1‑79 (MQ=255)                                                                                                                             ttAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAttt                                 >  1:163811/1‑89 (MQ=255)                                                                                                                             ttAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAttt                                 <  2:163811/89‑1 (MQ=255)                                                                                                                             ttAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                            <  2:1394380/94‑1 (MQ=255)                                                                                                                             ttAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                            >  1:1394380/1‑94 (MQ=255)                                                                                                                                |                                                                                                                     ACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTG  >  NC_000913/4447764‑4448007

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