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

CAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTT  >  NC_000913/309245‑309483                                                                                                                  |                                                                                                                              cagcGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCAc                                                                                                      >  1:1282487/1‑139 (MQ=255)         gCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCa                                                                                                                        >  2:1508373/1‑113 (MQ=255)         gCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCa                                                                                                                        <  1:1508373/113‑1 (MQ=255)                        aacAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCt                                                                               <  2:1282487/139‑1 (MQ=255)                         acAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTA‑TTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCttt                                                                             <  1:1616100/139‑1 (MQ=255)                                    aTGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCa                                                                   >  2:104010/1‑139 (MQ=255)                                     tGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAg                                                                  >  2:1554100/1‑139 (MQ=255)                                           aGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCa                                                            >  1:458081/1‑139 (MQ=255)                                             gCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCa                                                                                                                        <  1:917358/77‑1 (MQ=255)                                             gCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCa                                                                                                                        >  2:917358/1‑77 (MQ=255)                                                       cTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAac                                                >  2:305090/1‑139 (MQ=255)                                                       cTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAac                                                >  2:1167161/1‑139 (MQ=255)                                                         ttCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAacga                                              <  1:104010/139‑1 (MQ=255)                                                             gTCCTGGGGAGAACCGTGGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGa                                                                                                               >  2:1285678/1‑70 (MQ=255)                                                             gTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGa                                                                                                               <  1:1285678/70‑1 (MQ=255)                                                             gTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCaa                                                                                       <  2:38592/94‑1 (MQ=255)                                                             gTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCaa                                                                                       >  1:38592/1‑94 (MQ=255)                                                               ccTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGAc                                        <  2:236632/139‑1 (MQ=255)                                                                cTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACt                                       <  1:305090/139‑1 (MQ=255)                                                                   gggAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCAccc                                                                                                    <  2:972876/75‑1 (MQ=255)                                                                   gggAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCAccc                                                                                                    >  1:972876/1‑75 (MQ=255)                                                                        aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAgtg                                                                <  1:1005544/106‑1 (MQ=255)                                                                        aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAgtg                                                                >  2:1005544/1‑106 (MQ=255)                                                                        aaCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGc                               >  2:87845/1‑139 (MQ=255)                                                                         aCCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCa                              <  1:1554100/139‑1 (MQ=255)                                                                              cgtAAAGGCATTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAgg                         <  1:1167161/139‑1 (MQ=255)                                                                                  aaGGCATTATCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATG‑ACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtagta                    >  2:1338990/1‑139 (MQ=255)                                                                                       aTTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                      <  1:469357/69‑1 (MQ=255)                                                                                       aTTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                      >  2:469357/1‑69 (MQ=255)                                                                                       aTTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                      <  1:1199406/69‑1 (MQ=255)                                                                                       aTTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                      >  2:1199406/1‑69 (MQ=255)                                                                                       aTTGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtc                >  1:1010508/1‑139 (MQ=255)                                                                                        ttGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGTAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtct               >  2:188560/1‑139 (MQ=255)                                                                                        ttGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAAt                  <  1:1616294/136‑1 (MQ=255)                                                                                        ttGTCGTTATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAAt                  >  2:1616294/1‑136 (MQ=255)                                                                                              ttATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTac         <  1:87845/139‑1 (MQ=255)                                                                                              ttATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTac         <  2:1010508/139‑1 (MQ=255)                                                                                               tATCCATCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca        >  1:1047337/1‑139 (MQ=255)                                                                                                     tCTGCACAGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCtt  >  1:763165/1‑139 (MQ=255)                                                                                                         cacaGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCa                           <  2:1625523/110‑1 (MQ=255)                                                                                                         cacaGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCa                           >  1:1625523/1‑110 (MQ=255)                                                                                                           caGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                      <  2:1880394/49‑1 (MQ=255)                                                                                                           caGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                      >  1:1880394/1‑49 (MQ=255)                                                                                                           caGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                                   <  2:1775004/100‑1 (MQ=255)                                                                                                           caGCTATTTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGcc                                   >  1:1775004/1‑100 (MQ=255)                                                                                                              ctatttatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                      <  2:1809414/46‑1 (MQ=255)                                                                                                              ctatttatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAAc                                                                                      >  1:1809414/1‑46 (MQ=255)                                                                                                              ctatttatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCa                              >  1:1551787/1‑102 (MQ=255)                                                                                                              ctatttatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCa                              <  2:1551787/102‑1 (MQ=255)                                                                                                                  |                                                                                                                              CAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTT  >  NC_000913/309245‑309483

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