Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	305,502	C→T	R439H (CGC→CAC)	ecpD ←	fimbrial adhesin EcpD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	305,502	0	C	T	97.8%	144.4 / ‑3.4	45	coding (1316/1644 nt)	ecpD	fimbrial adhesin EcpD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/23);  minor base . (0/1);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

CATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGG  >  NC_000913/305365‑305635                                                                                                                                          |                                                                                                                                      cATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAgtgc                                                                                                                                     >  1:969455/1‑140 (MQ=255)  aTCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGtgcc                                                                                                                                    <  1:725152/140‑1 (MQ=255)  aTCGGTCATATCCCGCCAGCTAACATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGtgcc                                                                                                                                    <  2:614270/140‑1 (MQ=255)       tCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCa                                                                                                                               <  1:149856/140‑1 (MQ=255)         atatCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGtgccaatcacc                                                                                                                             >  1:429407/1‑133 (MQ=255)            tCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCAc                                                                                                                          >  2:760612/1‑140 (MQ=255)            tCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCAc                                                                                                                          <  1:1213666/140‑1 (MQ=255)                gCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTgc                                                                                                                      >  1:850130/1‑140 (MQ=255)                 ccAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTgcg                                                                                                                     <  2:969455/140‑1 (MQ=255)                 ccAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTgcg                                                                                                                     <  1:1145192/140‑1 (MQ=255)                 ccAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTgcg                                                                                                                     <  1:1188567/140‑1 (MQ=255)                 ccAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTgcg                                                                                                                     <  1:764488/140‑1 (MQ=255)                 ccAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTgcg                                                                                                                     >  1:617918/1‑140 (MQ=255)                          atcGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCgg                                                                                                            <  1:1061900/140‑1 (MQ=255)                          atcGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCgg                                                                                                            <  2:671493/140‑1 (MQ=255)                          atcGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCgg                                                                                                            <  1:220179/140‑1 (MQ=255)                               aCCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCAcc                                                                                                       <  1:522549/140‑1 (MQ=255)                                cccGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCgg                                                                                                            <  2:537143/134‑1 (MQ=255)                                cccGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCgg                                                                                                            >  1:537143/1‑134 (MQ=255)                                cccGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCt                                                                                                      >  1:401129/1‑140 (MQ=255)                                                       cTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCt                                                                               >  1:343710/1‑140 (MQ=255)                                                        tCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTg                                                                              >  1:314101/1‑140 (MQ=255)                                                            ttGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGtttt                                                                          >  1:866079/1‑140 (MQ=255)                                                               cGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAG‑‑‑‑TGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCtg                                                                   <  2:429407/140‑1 (MQ=255)                                                                gggTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCg                                                                      <  2:1111971/140‑1 (MQ=255)                                                                   taataaAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCtg                                                                   <  2:376079/140‑1 (MQ=255)                                                                   taataaAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCtg                                                                   <  2:351096/140‑1 (MQ=255)                                                                     ataaAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCtggt                                                                 <  2:343710/140‑1 (MQ=255)                                                                         aGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGtt                                                             <  2:866079/140‑1 (MQ=255)                                                                            aaaGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACg                                                          >  1:97568/1‑140 (MQ=255)                                                                                tgtCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTa                                                            <  2:934154/134‑1 (MQ=255)                                                                                tgtCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTa                                                            >  1:934154/1‑134 (MQ=255)                                                                                tgtCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGatat                                                      >  1:262817/1‑140 (MQ=255)                                                                                tgtCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGatat                                                      >  1:1173126/1‑140 (MQ=255)                                                                                tgtCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGatat                                                      >  1:825394/1‑140 (MQ=255)                                                                                    gCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGcc                                                  >  1:295486/1‑140 (MQ=255)                                                                                                       ccTTACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTcg                               <  2:295486/140‑1 (MQ=255)                                                                                                         ttACCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTcgcc                             >  1:1198014/1‑140 (MQ=255)                                                                                                           aCCGTCATCGGAGCTGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTcgccgc                           <  2:825394/140‑1 (MQ=255)                                                                                                                         tGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATAtccc             >  1:1043755/1‑140 (MQ=255)                                                                                                                         tGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATAtccc             >  1:1204057/1‑140 (MQ=255)                                                                                                                         tGAAGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATAtccc             >  1:899353/1‑140 (MQ=255)                                                                                                                           aaGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCt                          >  1:956329/1‑125 (MQ=255)                                                                                                                           aaGACACAATAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCt                          <  2:956329/125‑1 (MQ=255)                                                                                                                                   attggAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGGCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGGGTGCCGACTCGCGGCTACCAACATATCCCTCCCGCGACg   <  2:800080/137‑1 (MQ=255)                                                                                                                                   aTAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCCCCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACg   <  2:262817/140‑1 (MQ=255)                                                                                                                                   aTAGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACg   <  2:1127669/140‑1 (MQ=255)                                                                                                                                     aGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACgg  <  2:777148/139‑1 (MQ=255)                                                                                                                                     aGGAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACgg  >  1:777148/1‑139 (MQ=255)                                                                                                                                       gAGTGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACg   <  2:893954/136‑1 (MQ=255)                                                                                                                                       gAGTGCCCGCCAATCACCTGCGCGGGGCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGGGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACg   >  1:893954/1‑136 (MQ=255)                                                                                                                                          |                                                                                                                                      CATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGG  >  NC_000913/305365‑305635

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