Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,974,062	0	T	G	54.2%	19.0 / 29.8	38	L684R (CTG→CGG)	sasA_19	Adaptive‑response sensory‑kinase SasA
Reads supporting (aligned to +/- strand):  ref base T (4/7);  new base G (4/9);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-02
Rejected as consensus: Frequency below/above cutoff threshold.

AAGCCTTGATTGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCTGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCTTTTAC  >  ACB122/5973930‑5974187                                                                                                                                      |                                                                                                                               aaGCCTTGATTGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCTg                                                                                                                               <  5:193515/134‑1 (MQ=255)  aGCCTTGATTGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAGGc                                                                                                                             >  1:362404/1‑135 (MQ=255)        gATTGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATCGGCCCGGGcccgcg                                                                                                                        >  5:247797/1‑131 (MQ=255)        gATTGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCCGGCCCGgcgc                                                                                                                      >  7:381680/1‑136 (MQ=255)        gATTGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAGGcccagcg                                                                                                                       >  5:411094/1‑130 (MQ=255)          ttGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCCGGGCCAGCGcc                                                                                                                     >  6:320718/1‑135 (MQ=255)                 tGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GGCCAGGGCATGCGCCACGCCTg                                                                                                              >  4:115727/1‑135 (MQ=255)                  ggCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAGGCCCAGCGCCCGGCGTgg                                                                                                             >  4:401905/1‑135 (MQ=255)                       ccTGCAAACGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCCCCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAGGCCCATCGCCCGCGTTGGCCAgg                                                                                                        >  4:73624/1‑135 (MQ=255)                          gCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCGCGCCATGCGCCCGGCTTGGGCAggggg                                                                                                     >  4:339595/1‑133 (MQ=255)                           cAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCTGGCCCTGCGCCCCGCTTGGCCACGGCGa                                                                                                    >  6:90291/1‑135 (MQ=255)                           cAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAGGCCAATCGGCAGGCTTGTCCAGGGCGaa                                                                                                   >  8:334670/1‑136 (MQ=255)                              ccGGCCCAGACGGGGTGGGGTTTCACTCTTCACCGTGTTATGGCGGCCGCTTAGGAACTCCGGACCGGGCACCACCGCTCGTCGCGGAATGCAGTTATC‑TCCCTGGCCATGCGCCAGGCTTGGCCAGGGGGAACGc                                                                                                <  8:122504/136‑1 (MQ=255)                                 ggCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGAGGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCGGGCCATGCGCCAGGCTGGGCCAGGGCGAACGCAc                                                                                              <  4:346463/135‑1 (MQ=255)                                     aGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAGGCCCTGCGCCCGGCGTGGCCCGGGCGCACGCCCGAgg                                                                                          >  4:411104/1‑135 (MQ=255)                                     aGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAGGCCCAGCGCCACGCGTGGGCAGGGCGCCCGCACGAgg                                                                                          >  4:417951/1‑135 (MQ=255)                                               gggCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAGGCCAAGCGCCACCCGTGGGCCCGGCGCACCCCCCCGGTGTTGTTGGa                                                                                >  1:356424/1‑135 (MQ=255)                                               gggCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAAGCCCAGCCGCACGCGTGGCCAGGGCCCACCCCCGAGGTGGTGTTGGAg                                                                               >  7:120026/1‑136 (MQ=255)                                                                    gATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCTgg                                                                                                                              >  4:69595/1‑68 (MQ=255)                                                                    gATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCTgg                                                                                                                              <  3:69595/68‑1 (MQ=255)                                                                     aTGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCAGGCCCAGCGCCACGCGTGGCCAGGGCGGCCGCACGAGGTGGTGTTAGGGCTGGCCCGCCCCCTA‑CAAccc                                                         >  3:90532/1‑133 (MQ=255)                                                                      tGGCGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCTGGCCATGCGCCAGGCTGGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAAcgcg                                                        <  3:411104/135‑1 (MQ=255)                                                                      tGGCGGCCGCTTATGCACGCACGCCCGTGCACCCGCGCGCCTCGCTGCATGCACGTATC‑GCCCTGGCCATGCGCCCGGCTTGGCCCGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAAcgcg                                                        <  3:417951/135‑1 (MQ=255)                                                                              gCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCTGGCCCTGCGCCACGGTTGGGCAGGGGGCACGCACGCGGTGGTGGTGGCGCTCGCCCCCCCGCGC‑CACCGCGCCACcccc                                                >  2:13361/1‑135 (MQ=255)                                                                              gCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCGGGGCCGGCGCCCGGCGTGGGCCGGGGGCCCGCCCGCGGTGGTGTTCGGGCGCGCCCGCCCGCGA‑CAACGCGccccccccc                                               >  7:45554/1‑135 (MQ=255)                                                                               cTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCGGGGCATTCGCCACGCGTGGGCAGGGCGTGCGCGCGAGGTGTTGTTGGCGCTCGTCACGCAGCGAGCAACGC‑CCAGccccc                                               >  4:317740/1‑134 (MQ=255)                                                                                            cccggggCACCCGCGCGCCGCGCGGCAGGCAGGTATC‑CCCCGGCCCAGGCGCCAGGCTGGGCCAGGGCGAACGCACGAGGTGGTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGc                                  <  3:103585/129‑1 (MQ=255)                                                                                             cccggccACCCGCGCCCGGCGCTGCAGCCACTTACC‑CCCCGGGCCAGGCGCCTGGCGGGGCGTGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCg                                 <  3:274965/129‑1 (MQ=255)                                                                                             cccGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCCGGGCCAGCCCCCCCCGTGGGCCGGGCGCACCCACGCGGGGGTGTTGGAGCTCGCCCCCCCCCGA‑CACCCCGCCCGCCCCCGCCGGGTTGCcgcg                                 >  2:25361/1‑135 (MQ=255)                                                                                             cccGGGCACCCGCGCGCCTCGCTGCATGCACGTATC‑GCCCGGGCCAGGCGCCCGGCTGGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCg                                 <  4:30379/135‑1 (MQ=255)                                                                                                        gcgcgcCTGGCGGCAGGCAGGGAGC‑GCCCGGGCCAGGCGCCAGGCTTGGCCCGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGaa                     <  8:75887/136‑1 (MQ=255)                                                                                                           cgcCTCGCTGCATGCACGTATC‑GGCCTGGCCCAGCGCCCGGCGTGGGCCGGGGGCACCCCCCAGGGGGTGTTGGAGCTCGCCAGCCCCCTA‑CAACGCCCCCCCCCCAGCGGGGTTGGAGCCGCTTCCTGCGAAGCg                  >  8:20536/1‑136 (MQ=255)                                                                                                               tCGCTGCATGCACGTATC‑GCCCTGGCCATGCGCAAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTc               <  2:102844/135‑1 (MQ=255)                                                                                                                cGCTGCATGCACGTCCC‑GCCCGGCCCCGGCGCCTGGGTGGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCt             <  8:69578/136‑1 (MQ=255)                                                                                                                  cggcatgcaCGTACC‑GCCCGGGCCATGCGCCAGGCGGGCCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTcc              <  5:90291/131‑1 (MQ=255)                                                                                                                    gcctgcaCGTAGC‑GCCCGGGCCATGCGCCCGGCTGGGCCCGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAg         <  8:123539/133‑1 (MQ=255)                                                                                                                    gcaggcaCGTAGC‑GCCCGGGCCATGCGCCAGGCTTGGCCCGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAg         <  8:45554/136‑1 (MQ=255)                                                                                                                            gTATC‑GCCCTGGTCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCTTTTAc  <  3:47666/135‑1 (MQ=255)                                                                                                                                      |                                                                                                                               AAGCCTTGATTGGCGCTGGCTTCCTGCAACCGGGCCAGAACGGCTTGGGCTTCACTCATGACCGGGTGATGGAGGCCGCTTATGCACTCACGCCCGTGCACCAGCGCGCCTCGCTGCATGCACGTATC‑GCCCTGGCCATGCGCCAGGCTTGGCCAGGGCGAACGCACGAGGTGTTGTTCGAGCTCGCCAGCCAGCTA‑CAACGCGCCAGCCCCTGCGGGTTTGCAGCGGATGACTGCGAAGCGTTCCTGCAGCTTTTAC  >  ACB122/5973930‑5974187

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