Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	322,719	0	T	C	72.0%	66.7 / 40.1	50	G80G (GGA→GGG)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (7/7);  new base C (18/18);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCG  >  NZ_CP009273/322583‑322841                                                                                                                                         |                                                                                                                           ggTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGccc                                                                                                                          <  2:155783/139‑1 (MQ=255)  gtagtaGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCg                                                                                                                         <  1:587754/139‑1 (MQ=255)             ggcagTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCa                                                                                                                       <  2:395849/130‑1 (MQ=255)             ggcagTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCa                                                                                                                       >  1:395849/1‑130 (MQ=255)             ggcagTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATgcg                                                                                                              >  2:398610/1‑139 (MQ=255)              gcagTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATgcgg                                                                                                             >  1:47045/1‑139 (MQ=255)                  tCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGt                                                                                                         <  1:326282/139‑1 (MQ=255)                  tCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGt                                                                                                         <  1:334720/139‑1 (MQ=255)                              cAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGAACTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATaaa                                                                                               <  2:70628/137‑1 (MQ=255)                                   ctcCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGCTCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGtt                                                                                          >  2:228624/1‑137 (MQ=255)                                      cAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGt                                                                                     <  1:310907/139‑1 (MQ=255)                                         aCCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTcc                                                                                  >  1:766429/1‑139 (MQ=255)                                           cGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGt                                                                                >  1:501404/1‑139 (MQ=255)                                             ggTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGttt                                                                              >  1:756665/1‑139 (MQ=255)                                                     cgcCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGcc                                                                      <  1:476954/139‑1 (MQ=255)                                                            ttATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGt                                                               <  2:756665/139‑1 (MQ=255)                                                            ttATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTg                                                                >  1:563956/1‑138 (MQ=255)                                                            ttATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGt                                                               <  2:83165/139‑1 (MQ=255)                                                            ttATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGt                                                               <  2:317550/139‑1 (MQ=255)                                                            ttAGCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGt                                                               <  2:502760/139‑1 (MQ=255)                                                                 gagaTCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGt                                                          <  1:631754/139‑1 (MQ=255)                                                                 gagaTCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGt                                                          <  2:501404/139‑1 (MQ=255)                                                                       cAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTAcc                                                     <  1:313613/138‑1 (MQ=255)                                                                       cAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTAcc                                                     <  1:187430/138‑1 (MQ=255)                                                                        aGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAAc                                                                                              >  2:406650/1‑96 (MQ=255)                                                                        aGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAAc                                                                                              <  1:406650/96‑1 (MQ=255)                                                                          cgcATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTTTTCATAAACGGTTCCGGGTCCGGTTCATAGGCCCCGTTGTAGCGTTTAACCCGt                                                 >  2:631754/1‑139 (MQ=255)                                                                          cgcATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGt                                                          <  1:200108/130‑1 (MQ=255)                                                                          cgcATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGt                                                          >  2:200108/1‑130 (MQ=255)                                                                                            gCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCa                               <  2:563956/139‑1 (MQ=255)                                                                                                aTGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCCGCgg                           >  2:32030/1‑139 (MQ=255)                                                                                                aTGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCgg                           >  1:200166/1‑139 (MQ=255)                                                                                                 tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATg                                                                                                                <  1:414494/53‑1 (MQ=255)                                                                                                 tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATg                                                                                                                >  2:414494/1‑53 (MQ=255)                                                                                                 tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCggg                          >  1:432956/1‑139 (MQ=255)                                                                                                 tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGGAGCGGGAATGCCAGGGCAGCggg                          >  2:328566/1‑139 (MQ=255)                                                                                                 tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTa                                                       >  1:96329/1‑110 (MQ=255)                                                                                                 tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTa                                                       <  2:96329/110‑1 (MQ=255)                                                                                                 tGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCggg                          >  1:596488/1‑139 (MQ=255)                                                                                                  gCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGt                                                               <  2:391454/101‑1 (MQ=255)                                                                                                  gCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGt                                                               >  1:391454/1‑101 (MQ=255)                                                                                                  gCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCggg                          >  1:713746/1‑138 (MQ=255)                                                                                                  gCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGc                         >  2:684348/1‑139 (MQ=255)                                                                                                    cGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCAt                       <  2:596488/139‑1 (MQ=255)                                                                                                            aGCGACGATCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCCTCTGGGTGc               >  2:699602/1‑139 (MQ=255)                                                                                                              cgacgaTCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCgg             <  2:766429/139‑1 (MQ=255)                                                                                                                   aTCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGt         >  1:717159/1‑138 (MQ=255)                                                                                                                   aTCCACCCAGACCTTTACCGCGCCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGt         >  2:534501/1‑138 (MQ=255)                                                                                                                    tCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCa       <  2:432956/139‑1 (MQ=255)                                                                                                                             aCCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCg  >  2:403182/1‑135 (MQ=255)                                                                                                                             aCCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCg  <  1:403182/135‑1 (MQ=255)                                                                                                                                         |                                                                                                                           GGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTTCCGTTTCATAGGCCCAGTTGTAGCGTTTACCCTGTAGCGGGAATGCCAGGGCAGCGGGCATCTGGGTGCGGAAGTCAAAGCG  >  NZ_CP009273/322583‑322841

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