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.9%	47.4 / 26.5	41	T509P (ACC→CCC)	hyfB	hydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (21/2);  new base C (0/18);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 9.40e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGC  >  NC_000913/2603195‑2603487                                                                                                                                                    |                                                                                                                                                  cgcAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATAccc                                                                                                                                                <  1:605744/151‑1 (MQ=255)  gcAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATAcccc                                                                                                                                               <  2:687727/151‑1 (MQ=255)             gCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTg                                                                                                                                    >  2:506065/1‑151 (MQ=255)                gaggagcaacaGGAGGTGCAAAATACGATGTTCGTCGCCATGCTACTGCCCGCGGCATTCTGTTTATTATTTTCGTTTGGGGCTTGTGGGCTGCCCCCGAAGATAATGCATATTGCCCTTGCGTTTCCCAACCCCCCTGCCGCCACTGt                                                                                                                                   <  2:555409/141‑1 (MQ=255)                        acaGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCgg                                                                                                                          >  1:688964/1‑150 (MQ=255)                           ggaggtccaaaaTACGATGTTCGCCCCCATCCTCCTCCCCGCGGCCCTCTGCTTATTAATTGCGTTTAGTGCGGGTGGGCGGGCACGGAAAATAATGAATATTCCCCTTGCTTTTACCAACCCCCCCCCCGCCACTGTCGCCAGCGGa                                                                                                                         <  1:788376/139‑1 (MQ=255)                               gTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCAc                                                                                                                  >  1:207340/1‑151 (MQ=255)                                    aaaTAGGATGTTCGTCGCCTTTCTCCTGCTCGCGGCTCTCTGTTTTTTAATTGCGTTTAGTGCGATTGGGCTGCCCCAGAAGATAATGAATTTCGCCCGTGCTTTTCCAACCCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTg                                                                                                               <  1:365759/149‑1 (MQ=255)                                      aatcGATGACCGTCCCCATGCTACCGCTCGGGGCACTCTGCGTTTTAATTGCCTTTAGTGCCTGTTGGCTGGCACAAAAAATAATGCATTTTGCCCGCGCTTTTCCCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTAcc                                                                                                           <  1:177767/148‑1 (MQ=255)                                        aCGATGTCCGTCGCCATGCTCCTGCCCGCGCCACTCTGCTTATTATTTGCGTTTGGTGCAGGTGGGCGGCCACGAAAAAAAATGCATATTCCCCATGCTTTTACCAACCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCg                                                                                                         <  1:683136/151‑1 (MQ=255)                                         cGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCgg                                                                                                        >  1:711009/1‑151 (MQ=255)                                         cGATGACCGCCCCCTTCCTCCTGCCCGCGCCACTCGGCTTATTATTTGCTTTTGGTGCCAGTTGGCTGGCCCAAAAAATAATGAATATCCCCCTTGTTTTTACCAATCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCgg                                                                                                        <  1:501596/151‑1 (MQ=255)                                           atgatCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGc                                                                                                       >  1:766207/1‑150 (MQ=255)                                            tggtcgccgCCATGCTCCTCCCCGCGGCCCTCTGTGTATTATTGGCGTTTGGTGCCGGTGGGCTGGCCCCGAAGATAATGCATTTTCCCCTTGCGTTCACCAATCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCgg                                                                                                        <  1:7233/145‑1 (MQ=255)                                              atCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGt                                                                                                   >  1:121126/1‑151 (MQ=255)                                              atCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACCGTCGCCAGCGGAATAGCACTTGTACCCGGCACGt                                                                                                   >  1:195638/1‑151 (MQ=255)                                                    gCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCAt                                                                                              >  1:562590/1‑150 (MQ=255)                                                                cgcgcgCCCCTCTGCGTATAATTTGCGTTTAGTGCTTGTGGGCTGCCACAAAAGATAATGCATTTGGCCCTTGGTTTAACCATCCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCa                                                                                    <  1:770937/146‑1 (MQ=255)                                                                  cgcgcccttcggCGTATTAATTGCTTTTAGTGCCTGTTGGCGGGCACCGAAAATAATGCATATTGCCCATGCTTTTACAACCCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTc                                                                                     <  2:354231/134‑1 (MQ=255)                                                                       ctctctGCGTATTATTTGCGCTTAGTGCGAGTGGGCTGCCACCGAAGATAATGCATATTGCCCATGCTTTTACCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCt                                                                              <  2:129620/145‑1 (MQ=255)                                                                           ctGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTg                                                                      >  2:219033/1‑151 (MQ=255)                                                                                  ttAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGt                                                               >  2:462818/1‑151 (MQ=255)                                                                                    aaTTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGtt                                                              >  2:765216/1‑150 (MQ=255)                                                                                      ttGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACAAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCCTTACTGTTGCTGTtact                                                           >  2:662532/1‑151 (MQ=255)                                                                                        gcgtttgttccggtttggttggCACAAAAGATATGTAATTTTCCCCTTGCTTTTCCCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTtacta                                                          <  1:276143/136‑1 (MQ=255)                                                                                         cgttagggcccggTTGGCTGCCACGAAAAATACTGTATATCCCCAAGGCTTTTCCCAACCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCTCTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTtact                                                           <  2:95315/136‑1 (MQ=255)                                                                                         cgcTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTactact                                                        >  1:255959/1‑151 (MQ=255)                                                                                            ttggTGCCGGTGGGCGGGCACAAAAAATAATGCATTTTCCCCATGTTTTTACCAACCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTactact                                                        <  1:36362/145‑1 (MQ=255)                                                                                            ttAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGc                                                     >  1:639055/1‑151 (MQ=255)                                                                                              aGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCg                                                    >  1:28119/1‑150 (MQ=255)                                                                                                      tgggctggcACCGAAGATATTGCATATCGCCCTTGCTTTCCCCAACCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCg                                                    <  1:538423/140‑1 (MQ=255)                                                                                                       tgtctggcACCGAAGATAATCCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTtg                                           >  2:507999/4‑150 (MQ=255)                                                                                                        ggctgccACCAAAAATAATTCTTTTTGCCCTTGCTTTCCCCAATCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCttt                                            <  1:272168/148‑1 (MQ=255)                                                                                                           tggcACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGcc                                      >  2:431567/1‑151 (MQ=255)                                                                                                            gccACCAAAGATAATGAATATTCCCCGTGCTTTTACCACCCCCCCCCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTtgctg                                        <  1:417089/146‑1 (MQ=255)                                                                                                                ccGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGcc                                 >  1:478929/1‑151 (MQ=255)                                                                                                                  aaaGATAATGCATCTTGCCCAGGCTTTCCCCACTCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCtt                               <  1:321189/150‑1 (MQ=255)                                                                                                                   aaaaTAATGCATATTCCCCGTGCTTTCCCCAACCCCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTgg                                   <  2:689042/143‑1 (MQ=255)                                                                                                                        aaTGCATATTGGCCATGCGTTTACCAATACCCCTCCCGCCCCTGACGCCAGCGGAATACCCCTTGTACCCGGCACGATTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCGTTGCTGCCTGGCCTTTACTgg                         >  2:616358/1‑151 (MQ=255)                                                                                                                                ttGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTgtggtg                 >  1:76749/1‑151 (MQ=255)                                                                                                                                               ccAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCCTACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTcgcgccgcgc  >  2:824078/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CGCAGTACACACGCTGAAGAGGCACAGGAAGTGCCAAATACGATGATCGTCGCCATGCTACTGCTCGCGGCACTCTGCGTATTAATTGCGCTTAGTGCCAGTTGGCTGGCACCGAAGATAATGCATATTGCCCATGCGTTTACCAATACCCCTCCCGCCACTGTCGCCAGCGGAATAGCACTTGTACCCGGCACGTTTCATACACAGGTCACCCCCTCATTACTGTTGCTGTTACTACTGGCGATGCCTTTGCTGCCTGGCCTTTACTGGCTGTGGTGTCGTTCGCGCCGCGC  >  NC_000913/2603195‑2603487

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