Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,603,342	0	A	C	43.5%	51.9 / 30.4	46	T509P (ACC→CCC)	hyfB	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (24/2);  new base C (0/20);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 4.12e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCG  >  NC_000913/2603192‑2603491                                                                                                                                                          |                                                                                                                                                       ccgcGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATa                                                                                                                                                        <  2:495736/151‑1 (MQ=255)      cAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCt                                                                                                                                                  >  1:178541/1‑151 (MQ=255)        gTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTcc                                                                                                                                                <  1:474007/151‑1 (MQ=255)                      ggggCAAAGGAAGTCCAAAATACGATGATCGTCGCCATGCTCCTGCCCCCGCCCCTTTGTTTATTATTT‑GCGCTTAGTGCCAGTGGGCGGCCACGGAAAATAATGAATTTTCCCCTT‑GTGTTTACC‑ACCC‑CCCCTCCCGCCACTGTCGcc                                                                                                                                   <  2:567056/148‑1 (MQ=255)                      gagGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGcc                                                                                                                                   >  2:1980/1‑150 (MQ=255)                       agGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGcc                                                                                                                                   >  2:798299/1‑149 (MQ=255)                          cacaGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCgg                                                                                                                              >  2:615716/1‑151 (MQ=255)                                       aaaTACGTTGCTCGTCCCCTTGCTGCTGCCCGCGGCACTTTGTTTATTATTT‑GCCCTTAGTGCGAGTGGGCTGCCACGAAAAATAATCCATATCCCCCTT‑GC‑TTTACCAAACC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACt                                                                                                                     <  2:616430/147‑1 (MQ=255)                                       aaaTACGATGATCGCCGCCATGCTCCTCCCCGCGGCCCTCTGCGTATTATTT‑GCGTTTAGTGCGGTTGGGCTGCCCCAAAAAATAATGCATATGCCCCAT‑GCTTTTCCC‑AACC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACt                                                                                                                     <  1:281714/147‑1 (MQ=255)                                       aaaTACGATGACCGCCGCCATGCTCCTGCCCGCGCCCCTCTGTTTATTAATT‑GCGTTTGGTGCCGGTGGGCTGGCACGGAAAAAAATGAATTTTCCCCTT‑GCTTTTACC‑ATTC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACtt                                                                                                                    <  1:385159/148‑1 (MQ=255)                                                gatCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACg                                                                                                        >  1:701834/1‑151 (MQ=255)                                                gatCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACg                                                                                                        >  1:443541/1‑151 (MQ=255)                                                      cgCCATGCTCCTGCCCGCGCCTCTCTGTTTATTATTT‑GCTTTTGGTGCCGGTGGGCGGGCACAAAAAATAATGCATTTCCCCCTT‑GTGTTTCCC‑AACC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGtt                                                                                                      <  1:27103/147‑1 (MQ=255)                                                        ccATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATac                                                                                                >  2:294226/1‑151 (MQ=255)                                                              tACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGt                                                                                          >  2:331325/1‑151 (MQ=255)                                                                 tGCGCGCGCCACTTTGCGTTTTATTT‑GCGCTGAGTGCTGGTGGGGTGGCCCAGAAGATATTCTTTTTTGCACAG‑GCTTTTACC‑ACTC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATAcaca                                                                                             <  2:347595/145‑1 (MQ=255)                                                                      gcgGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCt                                                                                  >  2:337881/1‑151 (MQ=255)                                                                         gcctttcgttttttAATTGGCGTTTGGTGCTGGTGGGCGGCCACAGAAAAAAATACATATCCCCCTG‑GCGTTTACA‑ACCC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCa                                                                                <  2:733400/139‑1 (MQ=255)                                                                         cctttttgtttattatTT‑GCTCTTGGTCCCTGTGGGCGGCCACCAAAAATAAT‑CATTATGCCCCTGGCTTTTCCC‑AACC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCAt                                                                               <  2:369717/141‑1 (MQ=255)                                                                           aCTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCAt                                                                               >  2:918348/1‑149 (MQ=255)                                                                             tctGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCATTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACt                                                                           >  1:739229/1‑151 (MQ=255)                                                                               tGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGt                                                                         >  1:614406/1‑151 (MQ=255)                                                                                   ttttATTT‑GCGTTTGGGGCCGGTTGGCGGGCACCAAAAATAATTAATTTTCCCCTT‑GCTTTTACC‑AATC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTg                                                                       <  1:864464/147‑1 (MQ=255)                                                                                   tATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCt                                                                     >  2:909530/1‑151 (MQ=255)                                                                                   tATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATCCATATTGCACAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGACTAGCACTTGTACCCGGCACGTTTCATACACAGTTCACCCCCTCATTACTGt                                                                         >  2:782078/1‑147 (MQ=255)                                                                                    aTTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTg                                                                    >  2:504084/1‑151 (MQ=255)                                                                                        attggcttttggTGCAAGTTGGCGGGCACAAAAGATAAGGAATTTTGCCATT‑GCTTTTCCA‑AGCC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTta                                                                 <  1:494778/140‑1 (MQ=255)                                                                                                ttAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGc                                                         >  2:127104/1‑151 (MQ=255)                                                                                                tgggggcgggTTGGCGGGCACCAAAGATAATAAATATTCCCCTT‑GTTTTTCCC‑AACC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTta                                                                 <  1:446917/134‑1 (MQ=255)                                                                                                   gtgcgggggggcgggcACCAAAAAAAATGCATATTCCCCGT‑GCTTTTACC‑AATC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCg                                                        <  1:1007632/140‑1 (MQ=255)                                                                                                   gTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGAt                                                      >  2:79758/1‑151 (MQ=255)                                                                                                            gccggccacaaaaaaTAATGCATATTCCCCTT‑GCTTTTCCC‑AACC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTtgct                                             <  2:952288/137‑1 (MQ=255)                                                                                                                  cccgaaaGATATTGAATATTGCCCTT‑TCTTTTCCC‑AATC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTg                                        <  1:300400/145‑1 (MQ=255)                                                                                                                      aaaaaaaaTGCATGTCCCCCTT‑GTGTTTCCC‑ACCC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGc                                      <  2:759920/142‑1 (MQ=255)                                                                                                                        aaaTAATCCATATCCCCCTG‑GCTTTTACC‑AACC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTa                                 <  2:755333/149‑1 (MQ=255)                                                                                                                        aGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCttt                                  >  2:518118/1‑150 (MQ=255)                                                                                                                        aGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCttt                                  >  1:495726/1‑150 (MQ=255)                                                                                                                                  tttcccccAT‑GCTTTTCCC‑ACCC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCtgtg                        <  1:660390/145‑1 (MQ=255)                                                                                                                                    ttGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCtgtg                        >  1:374965/1‑148 (MQ=255)                                                                                                                                     tGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGgtgt                    >  1:203234/1‑151 (MQ=255)                                                                                                                                       cccTT‑GCTTTTACC‑AACC‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTc                   <  1:861882/150‑1 (MQ=255)                                                                                                                                       cccAT‑TCTTTTACC‑ACCC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCg                  <  1:615725/151‑1 (MQ=255)                                                                                                                                          ccgggcgttACC‑AACACCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTgtggtg                     <  2:635812/140‑1 (MQ=255)                                                                                                                                          aT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTc               >  1:446198/1‑151 (MQ=255)                                                                                                                                               gTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTcgcgc           >  2:187555/1‑151 (MQ=255)                                                                                                                                                ttcaCC‑AACC‑CCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTcg              <  1:375504/144‑1 (MQ=255)                                                                                                                                                         tA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCg  >  1:23159/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       CCGCGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATT‑GCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCAT‑GCGTTTACC‑AATA‑CCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGCAGCG  >  NC_000913/2603192‑2603491

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