Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,001,023	(A)8→7	coding (1808/2601 nt)	elfC →	putative fimbrial usher protein ElfC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,001,016	0	A	.	100.0%	227.5 / NA	54	coding (1801/2601 nt)	elfC	putative fimbrial usher protein ElfC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (28/26);  total (28/26)

TCGCTGGGGATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGGAAAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGC  >  NC_000913/1000890‑1001150                                                                                                                               |                                                                                                                                       tCGCTGGGGATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGcc                                                                                                                          >  1:691802/1‑140 (MQ=255) tCGCTGGGGATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGcc                                                                                                                          >  2:777858/1‑140 (MQ=255) tCGCTGGGGATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGcc                                                                                                                          <  1:68579/140‑1 (MQ=255)     tGGGGATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGa                                                                                                                      <  2:462935/140‑1 (MQ=255)        ggATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGa                                                                                                                   >  1:819906/1‑140 (MQ=255)        ggATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGa                                                                                                                   <  1:99921/140‑1 (MQ=255)        ggATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGa                                                                                                                   <  1:431550/140‑1 (MQ=255)        ggATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGa                                                                                                                   <  2:691802/140‑1 (MQ=255)         gATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGAc                                                                                                                  >  2:89006/1‑140 (MQ=255)          aTTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGAcc                                                                                                                 <  1:833321/140‑1 (MQ=255)             tCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACctct                                                                                                              <  1:841870/140‑1 (MQ=255)                        cAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACctct                                                                                                              >  2:218880/1‑129 (MQ=255)                        cAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACctct                                                                                                              <  1:218880/129‑1 (MQ=255)                             aGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATa                                                                                              >  1:722598/1‑140 (MQ=255)                                cccGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACg                                                                                           >  1:684239/1‑140 (MQ=255)                                cccGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACg                                                                                           <  1:851863/140‑1 (MQ=255)                                  cGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGaa                                                                                         >  2:573687/1‑140 (MQ=255)                                     ggTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAggg                                                                                      >  1:379586/1‑140 (MQ=255)                                     ggTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAggg                                                                                      <  1:796049/140‑1 (MQ=255)                                          aCCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCacac                                                                                 <  2:819906/140‑1 (MQ=255)                                                 gCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACaa                                                                          >  2:392959/1‑140 (MQ=255)                                                    gATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTACTGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAAc                                                                       <  2:753690/140‑1 (MQ=255)                                                     aTCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAAcc                                                                      >  2:336860/1‑140 (MQ=255)                                                     aTCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAAcc                                                                      >  1:796236/1‑140 (MQ=255)                                                     aTCAGGAGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGACACACTACACAAAAcc                                                                      >  2:413462/1‑140 (MQ=255)                                                         ggTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTggg                                                                  >  2:825313/1‑140 (MQ=255)                                                         ggTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTggg                                                                  <  1:89006/140‑1 (MQ=255)                                                          gtgtTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGt                                                                 >  1:861021/1‑140 (MQ=255)                                                           tgtTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGtt                                                                <  2:722598/140‑1 (MQ=255)                                                            gtTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGttt                                                               >  1:785678/1‑140 (MQ=255)                                                            gtTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGttt                                                               <  2:684239/140‑1 (MQ=255)                                                                gCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACg                                                           <  1:573687/140‑1 (MQ=255)                                                                    tAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAAc                                                                       <  1:36590/124‑1 (MQ=255)                                                                    tAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAAc                                                                       >  2:36590/1‑124 (MQ=255)                                                                    tAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGtt                                                                <  1:782430/131‑1 (MQ=255)                                                                    tAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATGACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGtt                                                                >  2:782430/1‑131 (MQ=255)                                                                            tCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTc                                               >  1:34290/1‑140 (MQ=255)                                                                               cTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGAt                                            <  2:785678/140‑1 (MQ=255)                                                                                        aaTCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAAc                                   <  1:587208/140‑1 (MQ=255)                                                                                        aaTCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAAc                                   <  1:765384/140‑1 (MQ=255)                                                                                          tCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCt                                 >  2:380610/1‑140 (MQ=255)                                                                                                             aaTGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAg              <  1:375983/140‑1 (MQ=255)                                                                                                              aTGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTa                           <  1:73000/126‑1 (MQ=255)                                                                                                              aTGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTa                           >  2:73000/1‑126 (MQ=255)                                                                                                              aTGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCg                      <  2:543430/131‑1 (MQ=255)                                                                                                              aTGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCg                      >  1:543430/1‑131 (MQ=255)                                                                                                              aTGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAg              >  2:481937/1‑139 (MQ=255)                                                                                                              aTGATAGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAg              <  1:481937/139‑1 (MQ=255)                                                                                                                   aGCTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGatat        >  1:831967/1‑140 (MQ=255)                                                                                                                     cTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAAc                >  2:281667/1‑130 (MQ=255)                                                                                                                     cTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAAc                <  1:281667/130‑1 (MQ=255)                                                                                                                     cTATACCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATaa      >  1:58249/1‑140 (MQ=255)                                                                                                                      tataCCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAAc     >  1:86951/1‑140 (MQ=255)                                                                                                                         aCCAGG‑AAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGc  >  2:33141/1‑140 (MQ=255)                                                                                                                               |                                                                                                                                       TCGCTGGGGATTTCCTGGAACTACAGTAAGTCCCGTGGTCAACCTGATGCTGATCAGGTGTTTGCTCTAAATTTTTCCCTGCCGCTCAATCTGTTGCTCCCCAGAAGTAATGATAGCTATACCAGGAAAAAAAATTACGCCTGGATGACCTCTAACACCAGTATCGATAACGAAGGGCACACTACACAAAACCTGGGTTTAACGGAGACATTACTCGATGACGGTAACCTGAGCTACAGCGTGCAACAGGGATATAACAGC  >  NC_000913/1000890‑1001150

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