Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,603,342	0	A	C	52.4%	13.6 / 66.7	42	T509P (ACC→CCC)	hyfB	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (16/4);  new base C (0/22);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 2.91e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-06

CGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGT  >  NC_000913/2603193‑2603473                                                                                                                                                        |                                                                                                                                    cgcgCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATAc                                                                                                                                    <  2:148955/151‑1 (MQ=255)   cgcAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATAccc                                                                                                                                  <  1:332938/151‑1 (MQ=255)       gTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTc                                                                                                                               >  2:334066/1‑150 (MQ=255)       gTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTc                                                                                                                               >  2:73855/1‑150 (MQ=255)       gTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGCTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTc                                                                                                                               >  2:108466/1‑150 (MQ=255)             acGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCAc                                                                                                                        <  2:319174/151‑1 (MQ=255)                cTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCa                                                                                                                         >  1:425691/1‑147 (MQ=255)                         cacaGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCgg                                                                                                            >  2:40823/1‑151 (MQ=255)                                      aaaTACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTa                                                                                               >  1:135875/1‑151 (MQ=255)                                       aaTACGATGTCCGTCGCCTTGCTCCCGCCCGCGCCACTCTGCTTATTAATTGCGTTTGGTGCCTTTGGG‑CGGGCCCCAAAAATAATGCATATCGCCCT‑TGCTTTTCCCAACCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTAc                                                                                              <  1:334057/151‑1 (MQ=255)                                       aaTACGATGAGCGCCGCCATGCTACTGCGCGCGCCCCTCTGCGTATTATTTGCGTTTGGTGCAGGTGGG‑CTGGCACAAAAAAAAATGAATATTGCCCA‑TGCGTTTCCCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTAc                                                                                              <  2:269639/151‑1 (MQ=255)                                             atgcccgtcgCCATGCCACTCCGCGCCCCCCTTCCCTTATTAATGGCCCTTGGTGCC‑GTTGGTCTGCCCCAGAAAAAAATGCATATTCCCGT‑TGCTTTAACCAACCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGc                                                                                         <  2:443072/145‑1 (MQ=255)                                                    tcgCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCa                                                                                 >  1:359278/1‑151 (MQ=255)                                                       ccATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTAccc                                                                                            >  2:224637/1‑137 (MQ=255)                                                         ttcctccggcgcgcgCCCCTCTGTTTATTATTTGCTCTTGGTGCCTGTTGG‑CGGCCACAAAAAATATTGCATTTTGCCCT‑TGCGTTTACCACCCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATacac                                                                            <  2:426554/140‑1 (MQ=255)                                                         cgcaaaCTGCTCGCGGCCCGGTGAATTTTACTCGCGTTTAGCGCCTGTGGG‑TAGGCCCAGAAGAAAATGCTTTTCCCCGG‑TGTCTTCAAAAACACCCCGGCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATacac                                                                            <  2:627/146‑1 (MQ=255)                                                         atcctccTGCCCGCGGCCCTCTGCTTATTATTGGCGTTTAGTGCAAGTGGG‑CGGGCCCCAAAGAAAATTCATATTGCCCA‑TGCTTTTCCAAATCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATacac                                                                            <  2:138625/145‑1 (MQ=255)                                                         aTGCTACTGCCCGCGGCCCTCTGCGTTTTAATTGCGCTTAGTGCCAGTGGG‑CGGCCACAAAAAATAATGCATTTTCCCCA‑TGCTTTTACCAACCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATacac                                                                            <  2:425702/151‑1 (MQ=255)                                                          tGCTACTGCTCGCGGCCCTCTGCGTATTATTGGCGCTGAGGGCCTGTGGG‑GTGGCACCGAAGATAATGCATATTCCCCG‑TGCTTTTCCCAATCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATacac                                                                            <  2:224165/150‑1 (MQ=255)                                                            ctccTGCTCGGGGCACTCTGCTTATTATTTGCTCTTGGTGCGAGTTGG‑CTGGCACAGAAGATAATGAATTTTGCCCA‑TGCTTTTCCCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAgg                                                                         <  2:337828/148‑1 (MQ=255)                                                            cTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAgg                                                                         >  2:381558/1‑151 (MQ=255)                                                               cTGCCCGCGGCACTCTGCTTATTATTTGCGTTTAGTGCCAGTTGG‑CTGCCCCGAAAAATAATGCATATTCCCCA‑GGCTTTTACAAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCa                                                                      <  2:245451/151‑1 (MQ=255)                                                                    cgcgCCACTCTGCGTATTATCTGCTCTTGGTCCCGGTTGC‑CCGGCCCCGAAAACAATGCTTATCCCCCA‑TGCTTTTACACCCCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCAccccc                                                                 <  1:297158/151‑1 (MQ=255)                                                                         cACTGTGTGTATTATTTGCCCTTGGTCCCGGTTGG‑CGGCCCCAAAAAAAAATCTATATCCCCCA‑TGCTTTTCCCACCCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCAtt                                                            <  1:86374/151‑1 (MQ=255)                                                                          cctctGCGTATTAATTGCGTTTAGTGCTAGTTGG‑CTGCCCCCGAAGATAATGCATATTCCCCA‑TGCTTTTCCCAATCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTa                                                           <  2:13867/150‑1 (MQ=255)                                                                           cttggcttATTAATTGCGTTTGGGGCCGTTGGG‑CGGGCCCCAAAAATAATGCATATCCCCCA‑TGCTTTTCCCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTAc                                                          <  1:73855/144‑1 (MQ=255)                                                                           ctctGGTTATTAATTCCTCTTGGTGCTGGTTGG‑CTGCCACAAAAAATAATGATTATTGCCCT‑TGCTTTTACCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTAc                                                          <  1:451594/151‑1 (MQ=255)                                                                               gCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGtt                                                      >  2:291403/1‑151 (MQ=255)                                                                               gCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑GGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGtt                                                      >  2:292825/1‑151 (MQ=255)                                                                                cttattatTGGCGTTTAGTGCGTGTGGG‑CGGGCACAGAAAATAATGCATTTTCCCCT‑TGTTTTTCCAAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTg                                                     <  2:409159/149‑1 (MQ=255)                                                                                 ttattatTGGCCCTTGGGGCCGGTGGG‑CGGGCCCCAAAAATAATGCATATTGCCTT‑TGCTTTTCCCACCCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTg                                                     <  1:42347/149‑1 (MQ=255)                                                                                  tATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTg                                                     >  2:80580/1‑149 (MQ=255)                                                                                       attccttttggtcccAGTGGC‑CCGGCACAAAATATATTCCATATTCCCCG‑TGCTTTCACCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTtac                                              <  1:183865/138‑1 (MQ=255)                                                                                          ccccTTGGTCCCAGTGGG‑CTGCCACCAAAGATATTGAATATCCCCCT‑TGCGTTCCCCCACCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTtactac                                           <  1:312879/148‑1 (MQ=255)                                                                                               tggTGCCGGTTGG‑CGGGCACCAAAGATAATGCATTTGCCCCA‑TGCTTTTCCAAATCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCg                                      <  2:451362/149‑1 (MQ=255)                                                                                                            gcggccccggaaaaaaaTGAATA‑TGCCCATTGCGTTCACAAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTtgct                           <  1:278963/135‑1 (MQ=255)                                                                                                                       aGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTa               >  2:27216/1‑151 (MQ=255)                                                                                                                       aGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTa               >  1:364420/1‑151 (MQ=255)                                                                                                                        aatatTGCTTATTCCCCA‑TGCTTTTCCCATCCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTAc              <  2:137417/150‑1 (MQ=255)                                                                                                                        aaTAATGCATATTCCCCA‑TGCTTTTCCCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTAc              <  2:123524/150‑1 (MQ=255)                                                                                                                              gCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTgg           >  1:321521/1‑148 (MQ=255)                                                                                                                                     gCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGgtgt  >  2:196474/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                    CGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGG‑CTGGCACCGAAGATAATGCATATTGCCCA‑TGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGT  >  NC_000913/2603193‑2603473

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