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.

CATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCT  >  NC_007779/357046‑357242
                                                                                                    |                                                                                                 
cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                   >  1:706154/1‑101 (MQ=255)
   ttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcggg                                                                                                >  1:1160370/1‑98 (MQ=255)
                     aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAGCTTTTAGccc                                                                             >  1:482284/1‑101 (MQ=39)
                     aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                             <  3:343615/101‑1 (MQ=255)
                             cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGAAAAAACTTTAAAccccccctgtt                                                                     >  3:623483/1‑93 (MQ=255)
                             cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAGATTTAGCCCGTCCtgtt                                                                     >  1:285294/1‑101 (MQ=255)
                               tCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAAGAAAACCTTTTATACCGTCTtgttgt                                                                   >  1:1052314/1‑101 (MQ=255)
                                      ggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                            <  3:392359/101‑1 (MQ=255)
                                       gcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACATTAAAACCCGCATGTTTTGTGGGaaa                                                           >  3:1138747/1‑100 (MQ=255)
                                         tGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACCTTAAACCCGTCCTGTTGTGTGGGGATCg                                                         >  1:314018/1‑101 (MQ=255)
                                          gTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAATAGCTTTTAACCCCTCCTCTTGTGTGGCGATCgc                                                        >  3:1591495/1‑101 (MQ=255)
                                           tCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAGCTTTAAACCCGTCATGTTGTGTGGGAAACgct                                                       >  4:36124/1‑101 (MQ=255)
                                            cTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAGATTTAAGCCCGGCATGTTGTGTGGCGATCgctg                                                      >  1:1464916/1‑101 (MQ=255)
                                              gCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                    <  1:839651/101‑1 (MQ=255)
                                               ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                   <  1:1484855/101‑1 (MQ=255)
                                                    tATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                              <  3:1542138/101‑1 (MQ=255)
                                                     aTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                             <  3:1215976/101‑1 (MQ=255)
                                                              ccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                    <  1:84305/101‑1 (MQ=255)
                                                                 caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAGCTTTAAAACCCTCATGTTGTGTGGGGATTGCTGCTGGTGACAATGTGGGAAggg                                 >  1:1211270/1‑101 (MQ=255)
                                                                  aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGAAAAAAAGTTAAGACCCCCAAGTTGTGGGGCCAGCTCCGCCTGTAACAATGTGGAAAgggg                                >  3:1490676/1‑100 (MQ=255)
                                                                    caaAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAATAGCATTTAGCCCCCCACTTTTTGTGGGACATGCTGCTGGTGACAATGTGGGAAGGGGAg                              >  1:1547328/1‑101 (MQ=255)
                                                                     aaaaCCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAGCTTTTAACCCCGCCTGTTGTGTGGCAAACGCTGCTTGTGACAAAGGGGGAAGGGGAga                             >  1:1259779/1‑100 (MQ=255)
                                                                      aaaCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                            <  3:713237/101‑1 (MQ=255)
                                                                                tctctcCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                  <  1:5883/101‑1 (MQ=255)
                                                                                      cAGTAGTGGCGGGGCGAGAAGAAGAACTTTAAAGCCGTCTCGTTTTGTGGCAATCCCTGCTGGTAACAATGTGGGAAGGGGAGGGCCGGGGATCAGcctcc            >  1:781222/1‑97 (MQ=255)
                                                                                      cAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa            <  1:1131615/101‑1 (MQ=255)
                                                                                       agtagtGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg           <  1:216690/101‑1 (MQ=255)
                                                                                         tagtGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa         <  1:839635/101‑1 (MQ=255)
                                                                                          agtGGCGGGGCAAGAGAAAAAAATTTAAACCCCCCATGTTTTGTGGGAATCGCTGCTGGTAACAATGTGGGAAGGGGAGGGGGGGGGAACAGCTCCCGaaa        >  1:114679/1‑100 (MQ=255)
                                                                                            tGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg      <  3:1280643/101‑1 (MQ=255)
                                                                                               cGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc   >  3:181591/1‑101 (MQ=255)
                                                                                                ggggCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt  <  3:38923/101‑1 (MQ=255)
                                                                                                    |                                                                                                 
CATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCT  >  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: