Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	309,358	C→T	*223* (TAG→TAA)	ecpB ←	ECP production pilus chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	309,358	0	C	T	100.0%	185.4 / NA	57	*223* (TAG→TAA)	ecpB	ECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/27);  total (30/27)

CGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTA  >  NC_000913/309231‑309484                                                                                                                                |                                                                                                                               cgcACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACggg                                                                                                                     >  2:1513276/1‑139 (MQ=255)            cagcagCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTAt                                                                                                          <  1:1000773/139‑1 (MQ=255)             agcagcGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGaa                                                                                                                   <  1:1175780/129‑1 (MQ=255)             agcagcGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGaa                                                                                                                   >  2:1175780/1‑129 (MQ=255)                agctcTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCAcc                                                                                                      >  2:436162/5‑139 (MQ=255)                agcGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCAcc                                                                                                      >  2:1286861/1‑139 (MQ=255)                agcGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCAcc                                                                                                      >  1:889062/1‑139 (MQ=255)                              gCTGAACGAACAAAAAGACCATGTCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCaa                                                                                        <  1:407670/139‑1 (MQ=255)                              gCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCaa                                                                                        <  2:1651892/139‑1 (MQ=255)                               cTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                       <  1:1513276/139‑1 (MQ=255)                                 gaacgaacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAccc                                                                                     >  1:1659015/1‑139 (MQ=255)                                    cgaacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTc                                                                                  <  2:889062/139‑1 (MQ=255)                                                ccATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCg                                                                      <  1:436162/139‑1 (MQ=255)                                                  aTGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCa                                                                    <  1:1286861/139‑1 (MQ=255)                                                   tGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAg                                                                   >  2:225839/1‑139 (MQ=255)                                                             aaaCTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCAcgcg                                                         >  1:511344/1‑139 (MQ=255)                                                               aCTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGa                                                                                                                >  2:962694/1‑82 (MQ=255)                                                               aCTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGa                                                                                                                <  1:962694/82‑1 (MQ=255)                                                                   ggCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTaaa                                                   >  1:533346/1‑139 (MQ=255)                                                                             ccTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGAc                                         >  2:508172/1‑139 (MQ=255)                                                                                  ggAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                                    >  2:757567/1‑139 (MQ=255)                                                                                    agaACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATc                                                                        <  2:869058/101‑1 (MQ=255)                                                                                    agaACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATc                                                                        >  1:869058/1‑101 (MQ=255)                                                                                      aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAAt                                                                                                                  <  2:1658724/57‑1 (MQ=255)                                                                                      aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAAt                                                                                                                  >  1:1658724/1‑57 (MQ=255)                                                                                      aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCg                                                                      >  2:867417/1‑101 (MQ=255)                                                                                      aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCg                                                                      <  1:867417/101‑1 (MQ=255)                                                                                      aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                                    <  1:1066090/135‑1 (MQ=255)                                                                                      aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                                    >  2:1066090/1‑135 (MQ=255)                                                                                       aCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCa                               <  1:1021198/139‑1 (MQ=255)                                                                                        gcgtcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGAc                                         >  1:1173441/2‑128 (MQ=255)                                                                                        ccGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTGATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGTGGGTAAAACGACGCGAc                                         <  2:1173441/128‑1 (MQ=255)                                                                                           tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAacga                                               >  1:215521/1‑119 (MQ=255)                                                                                           tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAacga                                               <  2:215521/119‑1 (MQ=255)                                                                                           tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCa                            <  2:892960/138‑1 (MQ=255)                                                                                           tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCa                            >  1:892960/1‑138 (MQ=255)                                                                                           tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAg                           >  2:885111/1‑139 (MQ=255)                                                                                           tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAg                           >  2:306975/1‑139 (MQ=255)                                                                                               aaaGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtag                       <  1:885111/139‑1 (MQ=255)                                                                                                   gCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAAt                   <  2:861553/139‑1 (MQ=255)                                                                                                      ttGTCGTTATTCATCTGCACAGCTATTTATTTCACGTGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCTGGCATCAGGTAGTAATgct                <  1:1657569/139‑1 (MQ=255)                                                                                                      ttGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtct                >  1:1311014/1‑139 (MQ=255)                                                                                                             tATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca         >  1:1864666/1‑139 (MQ=255)                                                                                                              aTCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTacac        >  1:1664938/1‑139 (MQ=255)                                                                                                                   tCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCa                            >  2:697666/1‑114 (MQ=255)                                                                                                                   tCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCa                            <  1:697666/114‑1 (MQ=255)                                                                                                                   tCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAAt                   <  1:824536/123‑1 (MQ=255)                                                                                                                   tCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAAt                   >  2:824536/1‑123 (MQ=255)                                                                                                                    cTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTa  <  1:306975/139‑1 (MQ=255)                                                                                                                       cacaGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCa                                                                    <  1:1037555/70‑1 (MQ=255)                                                                                                                       cacaGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCa                                                                    >  2:1037555/1‑70 (MQ=255)                                                                                                                          aGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTa           <  1:1596430/124‑1 (MQ=255)                                                                                                                          aGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTCCCCGGCATCAGGTAGTAATTCTCTTTa           >  2:1596430/1‑124 (MQ=255)                                                                                                                           gCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                       >  2:1833035/1‑47 (MQ=255)                                                                                                                           gCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                       <  1:1833035/47‑1 (MQ=255)                                                                                                                           gCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                                    <  2:1802854/98‑1 (MQ=255)                                                                                                                           gCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                                    >  1:1802854/1‑98 (MQ=255)                                                                                                                                |                                                                                                                               CGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTA  >  NC_000913/309231‑309484

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