Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,594	T→C	T247A (ACC→GCC)	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,594	0	T	C	100.0%	95.3 / NA	30	T247A (ACC→GCC)	ecpC	putative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/20);  total (10/20)

TGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCG  >  NC_000913/308480‑308723                                                                                                                   |                                                                                                                                  tGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGc                                                                                                                 >  1:464551/1‑133 (MQ=255) tGGCTGACCGGCTGCAGTCGAATATGGTGGGGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGc                                                                                                                 <  1:200581/133‑1 (MQ=255)   gCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCaa                                                                                                               >  1:346029/1‑133 (MQ=255)                     aaGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAAt                                                                                             <  2:462368/133‑1 (MQ=255)                            tggAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACAGGCAAATGGGTGACAGGCAAATTCGCGTTc                                                                                      >  1:247254/1‑133 (MQ=255)                                gCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATAc                                                                                  <  2:431570/133‑1 (MQ=255)                                 cTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACa                                                                                 <  2:346029/133‑1 (MQ=255)                                 cTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACa                                                                                 <  1:354850/133‑1 (MQ=255)                                 cTGGCCTGGTATCCCCAGGAAAGGTCGTAAATTTTCCCTGCGGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATata                                                                                 <  2:252554/133‑3 (MQ=255)                                  tGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATcg                                                                                           <  1:149626/122‑1 (MQ=255)                                  tGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATcg                                                                                           >  2:149626/1‑122 (MQ=255)                                    gCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCg                                                                              <  1:323122/133‑1 (MQ=255)                                      cTGGTTTCCCGAGGAAAGGCCGTAGATCGTCCATGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCt                                                                            <  2:230345/133‑1 (MQ=255)                                          tttCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTtat                                                                        <  2:378131/133‑1 (MQ=255)                                              cccAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTTATTCATCGCTTTATATaa                                                                    >  1:230345/1‑133 (MQ=255)                                                          cGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTg                                                        <  1:27081/133‑1 (MQ=255)                                                                tCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGCCTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATAGCTTTATATAATTCACTGTCCTGTTGAcc                                                  <  2:16257/133‑1 (MQ=255)                                                                tCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGAcc                                                  <  2:328119/133‑1 (MQ=255)                                                                  ttCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCCGGCGTCGAGCATTCCCCCCGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCCCTTTATAAAATACAGTGTCCTGTTGACCGc                                                >  1:216388/1‑133 (MQ=255)                                                                    cccTGCTGAAATGGCGGTCATCGGCCCCAAGGGCTGCAAGTTTCCGGCGGCGGGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTGTATGCATCGCTTTTCAGGATTCACTTGCCTGATGACCGcct                                              >  2:315282/1‑131 (MQ=255)                                                                                   ggTCATCGGCCCTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGc                               >  1:334879/1‑133 (MQ=255)                                                                                       aTCGGCCCTAAGGACAGCAAGTTCCAGGCGTCTAGCATTCCTCCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGc                           <  2:269590/133‑1 (MQ=255)                                                                                             ccTAAGGATTGCAAGTTCGAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGa                     <  1:315282/133‑1 (MQ=255)                                                                                              cTAAGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAg                    <  2:463641/133‑1 (MQ=255)                                                                                               tAAGGACTGCAAGTTCCAGGCGTCTAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGc                   <  2:269511/133‑1 (MQ=255)                                                                                                 aGGACTGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGcac                 >  2:129309/1‑133 (MQ=255)                                                                                                      tGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTATAGCTAGTCGTCGATCACCACat            >  1:211936/1‑133 (MQ=255)                                                                                                      tGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACat            <  2:507950/133‑1 (MQ=255)                                                                                                      tGCAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACat            >  1:154892/1‑133 (MQ=255)                                                                                                        cAAGTTCCAGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACatga          >  1:636068/1‑133 (MQ=255)                                                                                                                aGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTcgcg  <  2:334879/133‑1 (MQ=255)                                                                                                                aGGCGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTcgcg  <  1:610081/133‑1 (MQ=255)                                                                                                                aGGCGTCGAGAATTCCACCGGCAACTCGCTGATCGGCAAAAACGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTcgcg  <  2:247254/133‑1 (MQ=255)                                                                                                                   |                                                                                                                                  TGGCTGACTGGCTGCTGTCGAAGATGGTGGAGCTGGCCTGGTTTCCCCAGGAAAGGCCGTAAATCTTCCCTGCTGAAATGGCGGTCATCGGCCCTAAGGACTGCAAGTTCCAGGTGTCGAGCATTCCACCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCG  >  NC_000913/308480‑308723

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