Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	753,105	0	A	T	23.2%	143.5 / 11.1	56	intergenic (‑310/+80)	ybgD/gltA	putative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base A (33/10);  new base T (0/13);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 6.05e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGA  >  NC_000913/752970‑753239                                                                                                                                        |                                                                                                                                       agaCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttt                                                                                                                                    >  1:716168/1‑140 (MQ=255) agaCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttt                                                                                                                                    <  1:81775/140‑1 (MQ=255)   aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                  >  1:271511/1‑140 (MQ=255)   aCGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAAtttttt                                                                                                                                  >  1:383349/1‑140 (MQ=255)    cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                                 >  1:192193/1‑140 (MQ=255)    cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTtg                                                                                                                                 >  1:49983/1‑140 (MQ=255)        aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGtt                                                                                                                             >  1:534128/1‑140 (MQ=255)            tCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAgg                                                                                                                         <  1:95712/140‑1 (MQ=255)            tCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAgg                                                                                                                         <  2:68123/140‑1 (MQ=255)                   gTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGtt                                                                                                                             >  1:318407/1‑129 (MQ=255)                   gTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGtt                                                                                                                             <  2:318407/129‑1 (MQ=255)                      ataAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaaaa                                                                                                               >  1:809907/1‑140 (MQ=255)                            gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                         >  2:161468/1‑140 (MQ=255)                             cGACTACCATGAATTTTTTAATTGAAAGTATTGGGTTGCTGATATTTTGAGCTTTTCTTTTCTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTATTTTTACAATTTTTTGTTTTTAGGAATATAAAAATCAAccc                                                                                                        <  2:346963/140‑1 (MQ=255)                             cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                        >  1:481667/1‑140 (MQ=255)                             cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                        >  2:419784/1‑140 (MQ=255)                             cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                        >  1:540732/1‑140 (MQ=255)                                      tGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGa                                                                                               >  1:346963/1‑140 (MQ=255)                                       gAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGaa                                                                                              >  1:231593/1‑140 (MQ=255)                                                   taaaaGTATTGGTTTGCTAATATTTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTaaa                                                                                  <  2:383349/138‑1 (MQ=255)                                                   tGAAAGTATTGGGTTGCTGATAATTTGAGCTTTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                                  <  1:504734/140‑1 (MQ=255)                                                          tttGGGTTGTTAATAATTTGAGTTTTTTTTTTTTTTTAAAATATTTATATAGGTTTGTTAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTa                                                                           <  2:534128/139‑1 (MQ=255)                                                                ttGCTGATAATTTGAGTTTTTTTTTTTTTTTTAAATATCTATATAGGTCTGTTAATGTATTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                     <  2:49983/140‑1 (MQ=255)                                                                ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATCTTTTCAAAtt                                                                     >  2:654003/1‑140 (MQ=255)                                                                   cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                  >  1:459667/1‑140 (MQ=255)                                                                   cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                  >  1:207766/1‑140 (MQ=255)                                                                     aaTAATTTGAGTTGTTCTATTTTTTTTAAATATCTATATAGGTCTGTTAATGTTTTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCaa                                                                <  1:161468/139‑1 (MQ=255)                                                                        aaTTTGATCTTTTTTTTTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCa                                                             <  2:809907/140‑1 (MQ=255)                                                                                   ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGc                                                  >  1:818780/1‑140 (MQ=255)                                                                                        ttCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                             >  1:722905/1‑140 (MQ=255)                                                                                        ttCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                             >  1:604554/1‑140 (MQ=255)                                                                                           ttttaAAATATCTATATAGGTCTGTTAATGGATTTTTTTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATc                                          <  1:788789/140‑1 (MQ=255)                                                                                             tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGTGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGactgtct                                        >  1:491045/1‑134 (MQ=255)                                                                                             tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                        >  1:149209/1‑140 (MQ=255)                                                                                             tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                        >  1:770485/1‑140 (MQ=255)                                                                                             tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                        >  1:908010/1‑140 (MQ=255)                                                                                                aattatatataTAGGTTTGTAAATGTATTTTTTTTTAAAAATTTTTTGTTTTTGGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCttt                                     <  2:111431/137‑1 (MQ=255)                                                                                                     tCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAg                                <  1:503935/140‑1 (MQ=255)                                                                                                           taGGTTTTTTAATGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                          <  2:604554/140‑1 (MQ=255)                                                                                                           taGGTCTTTTAATGTTTTTTTTTTTAAAATTTTTTTGTTTTTAGGAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                          <  2:231593/140‑1 (MQ=255)                                                                                                              gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                       >  1:734982/1‑140 (MQ=255)                                                                                                               tCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                      >  1:233334/1‑140 (MQ=255)                                                                                                                ttGTTAATGTTTTTTTTTTTTACATTTTTTTGTTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                     <  1:506199/139‑1 (MQ=255)                                                                                                                 tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                    >  1:67300/1‑140 (MQ=255)                                                                                                                 tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                    >  1:563044/1‑140 (MQ=255)                                                                                                                  tttACTGGTTTTTTTTTTTAAATTTTTTTTTGTTTAGATATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  1:694721/139‑1 (MQ=255)                                                                                                                  tttAATGGTTTTTTTTTTTACATTTTTTTTTGTTTAGGAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  2:149209/139‑1 (MQ=255)                                                                                                                  tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  1:419784/139‑1 (MQ=255)                                                                                                                  tttAATGGTTTTTATTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  1:654003/139‑1 (MQ=255)                                                                                                                  tttAATGGTTTTTATTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                   <  1:761220/139‑1 (MQ=255)                                                                                                                   ttgaTGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGGTATCGCTTTTAAAATCGCGTTTTTCata                  >  2:524328/4‑140 (MQ=255)                                                                                                                   ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATTAAACATTAACGCTTGATTTTGCTTTTAAAGTCGCGGTTTTtata                  >  2:694721/1‑140 (MQ=255)                                                                                                                   ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAAGCTTGATATCCCTTTTTAAGACGCGTTTTTCtta                  >  2:219449/1‑140 (MQ=255)                                                                                                                     attGGTTTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAGGTCGCGTTTTTCATATc                <  2:459667/138‑1 (MQ=255)                                                                                                                     aatgttttttattttTACAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                <  2:67300/134‑1 (MQ=255)                                                                                                                     aaTGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGAAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                <  2:192193/140‑1 (MQ=255)                                                                                                                                   ttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATAcctgtct  >  1:835016/1‑134 (MQ=255)                                                                                                                                   ttACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGa  >  1:675067/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       AGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGA  >  NC_000913/752970‑753239

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