Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,857,1440CT52.2% 3.8 / 48.2 29L18F (CTC→TTC) ‡mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/5);  minor base G (6/5);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCA  >  NC_000913/2857034‑2857254
                                                                                                                |                                                                                                               
tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGAt                                                                                    >  2:227250/1‑139 (MQ=255)
                             aTTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCA‑‑CTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGAttttt                                                       <  2:598859/139‑1 (MQ=255)
                              tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTa                                                      >  1:737817/1‑139 (MQ=255)
                                  atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  <  1:227250/139‑1 (MQ=255)
                                  atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  <  2:109154/139‑1 (MQ=255)
                                  atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  >  1:347903/1‑139 (MQ=255)
                                  atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  <  2:740179/139‑1 (MQ=255)
                                    atCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  <  1:88652/137‑1 (MQ=255)
                                    atCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACt                                                >  2:73979/1‑139 (MQ=255)
                                    atCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  >  2:88652/1‑137 (MQ=255)
                                      cAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGt                                              >  1:142158/1‑139 (MQ=255)
                                      cAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGt                                              >  1:302027/1‑139 (MQ=255)
                                        gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGttt                                            >  1:669143/1‑139 (MQ=255)
                                        gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGttt                                            >  1:641905/1‑139 (MQ=255)
                                                           aTGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcgcg                         >  2:364485/1‑139 (MQ=255)
                                                               gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTc                     >  1:57393/1‑139 (MQ=255)
                                                                 gCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGc                   <  1:364485/139‑1 (MQ=255)
                                                                 gCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGc                   <  2:669143/139‑1 (MQ=255)
                                                                  cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCa                  >  2:281245/1‑139 (MQ=255)
                                                                  cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                                 <  1:296702/108‑1 (MQ=255)
                                                                  cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                                 >  2:296702/1‑108 (MQ=255)
                                                                  cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCa                  >  2:343521/1‑139 (MQ=255)
                                                                      aGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAActg              <  2:641905/139‑1 (MQ=255)
                                                                      aGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAActg              <  1:699737/139‑1 (MQ=255)
                                                                        aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATga                                       <  2:174476/112‑1 (MQ=255)
                                                                        aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATga                                       >  1:174476/1‑112 (MQ=255)
                                                                                 acgcccatacgcccatGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCgg                                                                    >  1:65013/1‑74 (MQ=255)
                                                                                 acgcccatacgcccatGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCgg                                                                    <  2:65013/74‑1 (MQ=255)
                                                                                  cgcccatacgcccatGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCa  <  2:302027/139‑1 (MQ=255)
                                                                                                                |                                                                                                               
TACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCA  >  NC_000913/2857034‑2857254

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