Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,141,009	0	G	A	31.0%	52.8 / 12.6	29	E28K (GAA→AAA)	yijC	DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (14/6);  new base A (0/9);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 7.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCT  >  REL606/4140934‑4141079                                                                            |                                                                       aTTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAg                                                                        <  1:1652671/76‑1 (MQ=255)  ttCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGa                                                                       <  2:1356307/76‑1 (MQ=255)    ctctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaa                                                                     >  2:939798/1‑76 (MQ=255)     tctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaa                                                                    <  1:1077431/76‑1 (MQ=255)      ctGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaa                                                                   >  1:703066/1‑76 (MQ=255)       tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGaaaaaa                                                                  <  1:173893/76‑1 (MQ=255)           ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg                                                              >  1:557648/1‑76 (MQ=255)           ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg                                                              >  1:481693/1‑76 (MQ=255)           ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg                                                              <  2:106277/76‑1 (MQ=255)           ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccg                                                              >  1:1470547/1‑76 (MQ=255)            taGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgc                                                             <  1:1505803/76‑1 (MQ=255)             aGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACccgcc                                                            <  1:1632283/76‑1 (MQ=255)                        tACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtgtgg                                                 >  1:1266735/1‑73 (MQ=255)                        tACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtg                                                 >  1:349480/1‑76 (MQ=255)                         aCTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtgg                                                >  2:103075/1‑76 (MQ=255)                         aCTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggtgg                                                >  2:559883/1‑76 (MQ=255)                             ttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGc                                            >  2:1621764/1‑76 (MQ=255)                              tGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGcc                                           >  2:657887/1‑76 (MQ=255)                                   aGGATTCAGACAGCGTGAGGGGCGAAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAtt                                      <  2:270196/76‑1 (MQ=255)                                    ggATTAAGACATGGTGAGGGGCGAAAGAGCACAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCAttt                                     <  1:280573/76‑1 (MQ=255)                                       tacAGACAGGGTGAGGGGCGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc                                  <  1:1023601/74‑1 (MQ=255)                                         cAGCCATCGTGATGGGGGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCa                                <  1:1437826/76‑1 (MQ=255)                                          aGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCaa                               >  2:1060892/1‑76 (MQ=255)                                            aCATCGTGAGGGGGGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt                             <  1:1396916/76‑1 (MQ=255)                                              atggTGAGGGGGGAAAGAGCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTaa                           <  1:356567/73‑1 (MQ=255)                                              aTCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTaa                           >  1:1426745/1‑76 (MQ=255)                                                    aggggCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTg                     <  2:1578440/74‑1 (MQ=255)                                                      gggCGTAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGaa                   <  1:1456727/76‑1 (MQ=255)                                                      gggCGAAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGaa                   <  1:1209898/76‑1 (MQ=255)                                                          gaaaGAGCCAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCa               <  2:1620100/74‑1 (MQ=255)                                                            aaGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGc             >  1:43042/1‑76 (MQ=255)                                                                   caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGCGCTGAACGCAGCTTCGCCa      >  2:360595/1‑76 (MQ=255)                                                                    aaaaaaaaaaaaaCCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAg     <  1:906730/73‑1 (MQ=255)                                                                     aaaaaaaaaaaaCCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGc    <  1:410594/74‑1 (MQ=255)                                                                       aaaaaaaaaaCCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCt  <  1:846928/76‑1 (MQ=255)                                                                       aaaaaaaaaaACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCt  <  1:31046/76‑1 (MQ=255)                                                                            |                                                                       ATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCT  >  REL606/4140934‑4141079

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