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%	163.3 / NA	52	*223* (TAG→TAA)	ecpB	ECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/28);  total (24/28)

CTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTAC  >  NC_000913/309228‑309476                                                                                                                                   |                                                                                                                       ctgcGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCAc                                                                                                                >  2:571476/1‑139 (MQ=255)   gcgcACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACgg                                                                                                              >  2:873555/1‑139 (MQ=255)      cACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGaa                                                                                                           >  1:1050090/1‑139 (MQ=255)        cTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATg                                                                                                         <  1:635624/139‑1 (MQ=255)                 gcagcGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCa                                                                                                <  1:931134/139‑1 (MQ=255)                 gcagcGCTGGCATCGGGCTGAACGAACAAAAAGACAATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCa                                                                                                >  1:390792/1‑139 (MQ=255)                     cgcTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTAt                                                                                                  <  2:611686/133‑1 (MQ=255)                     cgcTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTAt                                                                                                  >  1:611686/1‑133 (MQ=255)                        tGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGcc                                                                                         >  1:953972/1‑139 (MQ=255)                              cGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTg                                                                                   <  1:482358/139‑1 (MQ=255)                                  cTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                               <  2:953972/139‑1 (MQ=255)                                  cTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                               <  2:138195/139‑1 (MQ=255)                                  cTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                               <  2:1157722/139‑1 (MQ=255)                                    gaacgaacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAccc                                                                             <  2:616065/139‑1 (MQ=255)                                                   ccATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTCTCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCg                                                              <  1:873555/139‑1 (MQ=255)                                                    cATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATg                                                                                                         >  2:1038484/1‑95 (MQ=255)                                                    cATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATg                                                                                                         <  1:1038484/95‑1 (MQ=255)                                                     aTGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCa                                                            <  1:571476/139‑1 (MQ=255)                                                      tGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                               <  1:1012207/119‑1 (MQ=255)                                                      tGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                               >  2:1012207/1‑119 (MQ=255)                                                        ccGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAAt                                                                                                          >  1:318595/1‑90 (MQ=255)                                                        ccGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAAt                                                                                                          <  2:318595/90‑1 (MQ=255)                                                                                    gggAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGc                             >  1:827529/1‑139 (MQ=255)                                                                                     ggAGAACCGTCGTAAAGTCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                            <  2:1050090/139‑1 (MQ=255)                                                                                         aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCttt                                                                      <  1:1208766/93‑1 (MQ=255)                                                                                         aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCttt                                                                      >  2:1208766/1‑93 (MQ=255)                                                                                         aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGc                        >  1:890105/1‑139 (MQ=255)                                                                                         aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGc                        >  1:599803/1‑139 (MQ=255)                                                                                         aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGc                        >  2:1085696/1‑139 (MQ=255)                                                                                          aCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGa                                                                                                            <  2:996929/54‑1 (MQ=255)                                                                                          aCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGa                                                                                                            >  1:996929/1‑54 (MQ=255)                                                                                          aCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCa                       <  2:876439/139‑1 (MQ=255)                                                                                           ccGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCt                                                                        >  1:608065/1‑89 (MQ=255)                                                                                           ccGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCt                                                                        <  2:608065/89‑1 (MQ=255)                                                                                              tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCAccc                                                                                             <  1:399585/65‑1 (MQ=255)                                                                                              tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCAccc                                                                                             >  2:399585/1‑65 (MQ=255)                                                                                              tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCt                                                                        <  1:87207/86‑1 (MQ=255)                                                                                              tcgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCt                                                                        >  2:87207/1‑86 (MQ=255)                                                                                                gtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCa                       >  1:778285/1‑133 (MQ=255)                                                                                                gtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCa                       <  2:778285/133‑1 (MQ=255)                                                                                                        aTTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCa                                                            >  1:329316/1‑88 (MQ=255)                                                                                                        aTTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCa                                                            <  2:329316/88‑1 (MQ=255)                                                                                                        aTTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtc         <  2:890105/139‑1 (MQ=255)                                                                                                              gTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAG‑AATTCTCTTTac  <  2:599803/139‑1 (MQ=255)                                                                                                                 aTCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                               <  1:1157108/60‑1 (MQ=255)                                                                                                                 aTCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                               >  2:1157108/1‑60 (MQ=255)                                                                                                                 aTCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCTACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtct        >  1:682475/1‑131 (MQ=255)                                                                                                                 aTCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtct        <  2:682475/131‑1 (MQ=255)                                                                                                                      tCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCg                                                              >  1:799842/1‑72 (MQ=255)                                                                                                                      tCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCg                                                              <  2:799842/72‑1 (MQ=255)                                                                                                                            caGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAgg                  >  2:1095829/1‑110 (MQ=255)                                                                                                                            caGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAgg                  <  1:1095829/110‑1 (MQ=255)                                                                                                                                   |                                                                                                                       CTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTAC  >  NC_000913/309228‑309476

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