Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,986	Δ1 bp	coding (347/2526 nt)	ecpC ←	putative fimbrial usher protein EcpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	308,983	0	G	.	100.0%	159.9 / NA	37	coding (350/2526 nt)	ecpC	putative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (16/21);  total (16/21)

CAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCT  >  NC_000913/308910‑309058                                                                          |                                                                            cagcagcagTTGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACgggcgtt                                                                       <  1:69655/80‑5 (MQ=255)  agcagcagTTGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTgg                                                                    <  1:69656/81‑1 (MQ=255)       cagTTGGCGGAAGCTGAGATCCAGCTGGGAGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCaga                                                               <  1:69657/81‑1 (MQ=255)          ttGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGcc                                                            <  1:69659/81‑1 (MQ=255)          ttGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGcc                                                            <  2:69658/81‑1 (MQ=255)            ggCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCa                                                           <  1:69660/80‑1 (MQ=255)            ggCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCAt                                                          >  1:69722/1‑81 (MQ=255)            ggCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCAt                                                          >  1:69716/1‑81 (MQ=255)                 aaGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAgct                                                     <  2:69661/81‑1 (MQ=255)                 aaGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAgct                                                     >  1:69709/1‑81 (MQ=255)                    cTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAgctgct                                                  <  2:69662/81‑1 (MQ=255)                      gagaTCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGt                                                >  1:69741/1‑81 (MQ=255)                       agaTCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGATC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTc                                               >  2:69740/1‑81 (MQ=255)                       agaTCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGt                                                <  2:69665/80‑1 (MQ=255)                       agaTCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTc                                               <  2:69664/81‑1 (MQ=255)                       agaTCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTc                                               >  2:69726/1‑81 (MQ=255)                       agaTCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTc                                               <  1:69663/81‑1 (MQ=255)                        gaTCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTc                                               >  2:69744/1‑80 (MQ=255)                            cAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGtt                                          <  1:69667/81‑1 (MQ=255)                            cAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGt                                           <  1:69666/80‑1 (MQ=255)                               cTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGt                                       <  2:69668/81‑1 (MQ=255)                                  cgcgTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGtt                                      <  2:69669/79‑1 (MQ=255)                                    cgTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCt                                  <  2:69670/81‑1 (MQ=255)                                     gTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTg                                 <  1:69671/81‑1 (MQ=255)                                        gtcagtcaGGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGaca                              >  2:69735/1‑81 (MQ=255)                                              cagGGGGATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCg                        >  1:69731/1‑81 (MQ=255)                                                  gggATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTAt                    <  2:69673/81‑1 (MQ=255)                                                   ggATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATc                   >  2:69736/1‑81 (MQ=255)                                                    gATGGTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCt                  <  1:69674/81‑1 (MQ=255)                                                        gTAAGGGCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCAc              >  2:69721/1‑81 (MQ=255)                                                             ggCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCActctct         <  1:69675/81‑1 (MQ=255)                                                              gCCTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCtt        >  1:69747/1‑81 (MQ=255)                                                                cTCATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTcc      <  1:69676/81‑1 (MQ=255)                                                                  cATTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAg    >  2:69733/1‑81 (MQ=255)                                                                   aTTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGc   >  1:69720/1‑81 (MQ=255)                                                                   aTTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGc   >  1:69742/1‑81 (MQ=255)                                                                   aTTGAAC‑GGGCGTTGGCCAGAGCCATCAGCTACTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGc   <  2:69677/81‑1 (MQ=255)                                                                    ttGAAC‑GGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCt  >  1:69723/1‑81 (MQ=255)                                                                          |                                                                            CAGCAGCAGTTGGCGCAAGCTGAGATCCAGCTGCGCGTTGTCAGTCAGGGGGATGGTAAGGGCCTCATTGAACGGGGCGTTGGCCAGAGCCATCAGCTGCTGTCGAGTTTGTTCGCTGACACTGGCGTTATCTTCACTCTCTTCCAGCT  >  NC_000913/308910‑309058

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