Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,857,1450TA64.0% 21.3 / 26.8 25L18H (CTC→CAC) ‡mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base A (10/6);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTG  >  NC_000913/2857034‑2857242
                                                                                                                  |                                                                                                 
tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGAt                                                                        >  2:63108/1‑139 (MQ=255)
    aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCctg                                                                    >  1:552712/1‑139 (MQ=255)
    aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctg                                                                    <  1:618309/139‑1 (MQ=255)
    aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctg                                                                    >  1:34628/1‑139 (MQ=255)
                          cccATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctg                                                                    >  2:1009809/1‑117 (MQ=255)
                          cccATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctg                                                                    <  1:1009809/117‑1 (MQ=255)
                              tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTa                                          >  1:932620/1‑139 (MQ=255)
                                 aaTATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGa                                       <  2:34628/139‑1 (MQ=255)
                                   tatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                     <  2:552712/139‑1 (MQ=255)
                                   tatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                     >  1:530733/1‑139 (MQ=255)
                                        gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGttt                                                              <  1:472749/109‑1 (MQ=255)
                                        gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGttt                                                              >  2:472749/1‑109 (MQ=255)
                                        gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                     <  2:458521/134‑1 (MQ=255)
                                        gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                     >  1:458521/1‑134 (MQ=255)
                                        gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                      <  1:1314436/133‑1 (MQ=255)
                                        gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTAACGGATGGGTGATTTTTATGaa                                      >  2:1314436/1‑133 (MQ=255)
                                           aaCCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTAt                             >  2:509540/1‑139 (MQ=255)
                                              cGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATgac                          <  2:34041/139‑1 (MQ=255)
                                                           aTGAGTGTAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                      <  2:279047/114‑1 (MQ=255)
                                                           aTGAGTGTAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                      >  1:279047/1‑114 (MQ=255)
                                                                  cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                      <  1:231328/107‑1 (MQ=255)
                                                                  cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                      >  2:231328/1‑107 (MQ=255)
                                                                      aGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAActg  <  2:530733/139‑1 (MQ=255)
                                                                        aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATg                                        <  1:272842/99‑1 (MQ=255)
                                                                        aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATg                                        >  2:272842/1‑99 (MQ=255)
                                                                                                                  |                                                                                                 
TACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTG  >  NC_000913/2857034‑2857242

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