Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779357,1460CG43.3% 38.4 / 13.2 30A257P (GCC→CCC) cynRDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (3/14);  new base G (13/0);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 4.79e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCT  >  NC_007779/357046‑357242
                                                                                                    |                                                                                                 
cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                   >  1:706154/1‑101 (MQ=255)
   ttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcggg                                                                                                >  1:1160370/1‑98 (MQ=255)
                     aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAGCTTTTAGccc                                                                             >  1:482284/1‑101 (MQ=39)
                     aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGccc                                                                             <  3:343615/101‑1 (MQ=255)
                             cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGA‑GAAAAAACTTTAAAccccccctgtt                                                                     >  3:623483/1‑93 (MQ=255)
                             cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAGATTT‑AGCCCGTCCtgtt                                                                     >  1:285294/1‑101 (MQ=255)
                               tCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAA‑AGAAAACCTTTTATACCGTCTtgttgt                                                                   >  1:1052314/1‑101 (MQ=255)
                                      ggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                            <  3:392359/101‑1 (MQ=255)
                                       gcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAACATTAAAACCCGCATGTTTTGTGGGaaa                                                           >  3:1138747/1‑100 (MQ=255)
                                         tGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAACCTTAAACCCGTCCTGTTGTGTGGGGATCg                                                         >  1:314018/1‑101 (MQ=255)
                                          gTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAATAGCTTTTAACCCCTCCTCTTGTGTGGCGATCgc                                                        >  3:1591495/1‑101 (MQ=255)
                                           tCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAGCTTTAAACCCGTCATGTTGTGTGGGAAACgct                                                       >  4:36124/1‑101 (MQ=255)
                                            cTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAGATTTAAGCCCGGCATGTTGTGTGGCGATCgctg                                                      >  1:1464916/1‑101 (MQ=255)
                                              gCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                    <  1:839651/101‑1 (MQ=255)
                                               ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                   <  1:1484855/101‑1 (MQ=255)
                                                    tATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                              <  3:1542138/101‑1 (MQ=255)
                                                     aTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                             <  3:1215976/101‑1 (MQ=255)
                                                              ccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                    <  1:84305/101‑1 (MQ=255)
                                                                 caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAGCTTTAAAACCCTCATGTTGTGTGGGGATTGCTGCTGGTGACAATGTGGGAAggg                                 >  1:1211270/1‑101 (MQ=255)
                                                                  aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGA‑GAAAAAAAGTTAAGACCCCCAAGTTGTGGGGCCAGCTCCGCCTGTAACAATGTGGAAAgggg                                >  3:1490676/1‑100 (MQ=255)
                                                                    caaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAATAGCATTTAGCCCCCCACTTTTTGTGGGACATGCTGCTGGTGACAATGTGGGAAGGGGAg                              >  1:1547328/1‑101 (MQ=255)
                                                                     aaaaCCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAGCTTTTAACCCCGCCTGTTGTGTGGCAAACGCTGCTTGTGACAAAGGGGGAAGGGGAga                             >  1:1259779/1‑100 (MQ=255)
                                                                      aaaCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                            <  3:713237/101‑1 (MQ=255)
                                                                                tctctcCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                  <  1:5883/101‑1 (MQ=255)
                                                                                      cAGTAGTGGCGGGGCGAGA‑AGAAGAACTTTAAAGCCGTCTCGTTTTGTGGCAATCCCTGCTGGTAACAATGTGGGAAGGGGAGGGCCGGGGATCAGcctcc            >  1:781222/1‑97 (MQ=255)
                                                                                      cAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa            <  1:1131615/101‑1 (MQ=255)
                                                                                       agtagtGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg           <  1:216690/101‑1 (MQ=255)
                                                                                         tagtGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa         <  1:839635/101‑1 (MQ=255)
                                                                                          agtGGCGGGGCAAGA‑GAAAAAAATTTAAACCCCCCATGTTTTGTGGGAATCGCTGCTGGTAACAATGTGGGAAGGGGAGGGGGGGGGAACAGCTCCCGaaa        >  1:114679/1‑100 (MQ=255)
                                                                                            tGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg      <  3:1280643/101‑1 (MQ=255)
                                                                                               cGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc   >  3:181591/1‑101 (MQ=255)
                                                                                                ggggCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt  <  3:38923/101‑1 (MQ=255)
                                                                                                    |                                                                                                 
CATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCT  >  NC_007779/357046‑357242

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