Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,138,764	0	A	T	55.6%	6.2 / 34.8	27	V286E (GTG→GAG)	hybO	hydrogenase 2 small subunit
Reads supporting (aligned to +/- strand):  ref base A (6/6);  new base T (7/8);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACA  >  NZ_CP009273/3138677‑3138848                                                                                        |                                                                                     tGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCTcc                                                                                    <  2:227528/90‑1 (MQ=255)        cGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCTCCGGCcaca                                                                             >  2:664299/1‑90 (MQ=255)        cGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCAcc                                                                                    >  2:492902/1‑83 (MQ=255)        cGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCAcc                                                                                    <  1:492902/83‑1 (MQ=255)          aCGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCcacaca                                                                           <  2:601910/90‑1 (MQ=255)          aCGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCcacaca                                                                           <  2:254003/90‑1 (MQ=255)            gTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACAcc                                                                         >  1:433555/1‑90 (MQ=255)                gCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCTCCGGCCACACACCGCCa                                                                     <  1:729438/90‑1 (MQ=255)                  aaGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAAc                                                                   <  1:227344/90‑1 (MQ=255)                                         aGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGc                                            >  1:227966/1‑90 (MQ=255)                                         aGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGc                                            >  2:264114/1‑90 (MQ=255)                                         aGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGc                                            >  2:672111/1‑90 (MQ=255)                                                   ttcttcGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTa                                  <  1:41117/90‑1 (MQ=255)                                                          ttACAGCCATAGCAAGGGTGACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCt                           <  1:664299/90‑1 (MQ=255)                                                             cAGCCATAGCAAGGGTGACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGc                        <  2:148742/90‑1 (MQ=255)                                                                         gggTGACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGt                                         >  2:672522/1‑61 (MQ=255)                                                                         gggTGACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGt                                         <  1:672522/61‑1 (MQ=255)                                                                         gggTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCg            >  2:597808/1‑90 (MQ=255)                                                                             gACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGa                                              >  1:847891/1‑52 (MQ=255)                                                                             gACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGa                                              <  2:847891/52‑1 (MQ=255)                                                                             gACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAggt        >  2:554573/1‑90 (MQ=255)                                                                             gACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAggt        >  2:370239/1‑90 (MQ=255)                                                                             gACCAATCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAggt        >  1:711734/1‑90 (MQ=255)                                                                             gACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGa           >  1:116612/1‑87 (MQ=255)                                                                             gACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGa           <  2:116612/87‑1 (MQ=255)                                                                                 aaTCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTa    <  2:573919/90‑1 (MQ=255)                                                                                   tCGCCTCCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACa  <  2:789175/90‑1 (MQ=255)                                                                                        |                                                                                     TGATTTTCGACGTTGGCAAGCTGATGGATGCCTTTATGGAAGCCGATACCTTCTTCGTTACAGCCATAGCAAGGGTGACCAATCGCCACCGGCCACACACCGCCAACATCGCAGAATTGCAGCGTTGAGCAGTTGCCGTAAGTTTCTGGCCCTTTACAGCCGAGGTGGTACA  >  NZ_CP009273/3138677‑3138848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑