Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,587,323	C→T	A138T (GCC→ACC)	ydeQ ←	putative fimbrial adhesin protein YdeQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,587,323	0	C	T	97.1%	107.9 / NA	34	A138T (GCC→ACC)	ydeQ	putative fimbrial adhesin protein YdeQ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/12);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACC  >  NC_000913/1587085‑1587567                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     agCTTTTGTTCGATTGAGCAATATACGCCAAGGGGTATTTATGCGCTTCCAGGAAAGTCAGGGAGATCCACGGTGAAATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCAACGGTCGCCCAAACTgg                                                                                                                                                                                                                                          >  2:419372/1‑251 (MQ=255)  gCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTggg                                                                                                                                                                                                                                         >  2:124273/1‑251 (MQ=255)  gCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTggg                                                                                                                                                                                                                                         >  1:387608/1‑251 (MQ=255)               tGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGtt                                                                                                                                                                                                                            >  2:78701/1‑251 (MQ=255)                    aaTATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAAc                                                                                                                                                                                                                       >  1:310608/1‑251 (MQ=255)                                            cgcTTCCAGGCAAGTCCAGGACATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCTCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGATTTAACGGTAATGGCATCGGCGTCTTATc                                                                                                                                                                                                >  1:401174/1‑250 (MQ=255)                                            cgcTTCCAGGAAAGTCCGGGAGATCCACGGTGGCATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTGGCAATGTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGGTTTAACGGTAATGGCATCGGTGTCTTATCa                                                                                                                                                                                               <  1:522896/251‑1 (MQ=255)                                                      aaaGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAAAATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATcc                                                                                                                                                                                      >  1:530003/1‑250 (MQ=255)                                                            cGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGt                                                                                                                                                                               >  1:909/1‑251 (MQ=255)                                                                            aCATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCaa                                                                                                                                                               <  1:170732/251‑1 (MQ=255)                                                                                                             ggCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAg                                                                                                                              <  1:124273/251‑1 (MQ=255)                                                                                                              gCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGc                                                                                                                             <  1:78701/251‑1 (MQ=255)                                                                                                                   acaacaTTGTTATTTGTAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATcc                                                                                                                        >  2:41458/1‑251 (MQ=255)                                                                                                                       caTTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTta                                                                                                                    <  1:463012/251‑1 (MQ=255)                                                                                                                                          atTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCt                                                                                                  <  2:202795/250‑1 (MQ=255)                                                                                                                                           ttccaattaacattcagcgGATTTCCGCTTCCTCATGTTGCACTTTTATCCATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACACCTACACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCACTACGTTGGTATTTGTCGTCACCGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGACTGc                                                                                                >  2:387831/17‑251 (MQ=255)                                                                                                                                             ccACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATTAATACCGGCAATGACTGCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATTGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTgtc                                                                                                                                                                    <  2:331247/181‑1 (MQ=255)                                                                                                                                             ccACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGAAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTgtc                                                                                                                                                                    >  1:331247/1‑181 (MQ=255)                                                                                                                                                    aaaTTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGa                                                                                       <  1:337829/251‑1 (MQ=255)                                                                                                                                                                           ttaTGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCACGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGCAGGGGTCGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGTAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGAt                                                                >  2:329075/3‑251 (MQ=255)                                                                                                                                                                             aTGTGGCAATTTGATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCAGCGGTGTCTTATCACCAATATCCAATTCGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTAGATGATCa                                                              <  2:373018/251‑1 (MQ=255)                                                                                                                                                                                    aaTTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGt                                                       >  1:157010/1‑251 (MQ=255)                                                                                                                                                                                    aaTTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGt                                                       <  2:316292/251‑1 (MQ=255)                                                                                                                                                                                             cATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAgccgc                                              >  1:78727/1‑251 (MQ=255)                                                                                                                                                                                                atGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCACCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGt                                           >  1:358342/1‑251 (MQ=255)                                                                                                                                                                                                   aaTACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTtat                                                                     >  1:278117/1‑222 (MQ=255)                                                                                                                                                                                                   aaTACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTAGATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCATGTTtat                                                                     <  2:278117/222‑1 (MQ=255)                                                                                                                                                                                                       cGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCAtt                                    >  1:437525/1‑251 (MQ=255)                                                                                                                                                                                                               taCTTCCCCCGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAAc                            >  2:500183/2‑251 (MQ=255)                                                                                                                                                                                                               gACTTCCCCCGCTTTTATGACAACTCCCCCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAAc                            >  2:297601/1‑251 (MQ=255)                                                                                                                                                                                                                 cgtCCCCCGCTTTTATGACAACTCCCCCGGCCGCCCAAACTGGGGTAATATAGAGTTTTACCGGGAATGGCAACGGTGTTTTATCCCCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGAGCCTT‑ATATTACTGTAGGGAGCCTGCGAAAGCTGGGCCTTGTACCAGGTTTATATGAGCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACtt                          <  2:402002/248‑1 (MQ=255)                                                                                                                                                                                                                       ccGCTTTTATGACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGc                    >  2:435966/1‑251 (MQ=255)                                                                                                                                                                                                                                 gACAACTCCACCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCCCCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGCCTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCGGCCGCCGTAATCATTCCAGCAACTTATATGCTGAGACAAGt          <  2:909/251‑1 (MQ=255)                                                                                                                                                                                                                                          aCCGGTCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAAcc  >  2:17967/1‑250 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     AGCTTTTGTTCGCTTGAGCAATATACGCCAAGGGGTATTTCTGCGCTTCCAGGAAAGTCCGGGAGATCCACGGTGACATTTCGACTATCAACGGTGCAGCCGCCTGTTGGCATAACAACATTGTTATTTGAAATAATATTCCACGTAAAATTACGCGGATTTCCGCTTCCTAATGTTGCAATTTTATACATATGAATACGGGCAATGACTTCCCCCGCTTTTATGACAACTCCACCGGCCGCCCCAACTGGGGTAATATAGAGTTTTAACGGTAATGGCATCGGTGTCTTATCACCAATATCCAATACGTTGGTATTTGTCGTCAACGGGAAGGGGTAGGTCACATTATTCCAGTACAGCGATCCTTTATATGACTGTAGGGAGCCTGCGAAAGCTGAGCCTTGTACCAGGTTTATATGATCAGTGTCGTACCAGCCGCCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACC  >  NC_000913/1587085‑1587567

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