Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,857,1430TA56.2% 2.4 / 47.1 32Y17* (TAT→TAAmutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (7/7);  new base A (10/8);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTG  >  NC_000913/2857027‑2857257
                                                                                                                     |                                                                                                                    
aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTGGCTGAAAGCCCAg                                                                                                  <  2:763745/135‑1 (MQ=255)
               aaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGt                                                                               >  2:947253/1‑139 (MQ=255)
                    aaaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTAc                                                                          >  1:950424/1‑139 (MQ=255)
                           tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATgg                                                                   <  1:947253/139‑1 (MQ=255)
                                 cccATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCa                                                                                                >  2:824418/1‑104 (MQ=255)
                                 cccATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCa                                                                                                <  1:824418/104‑1 (MQ=255)
                                     tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTa                                                         >  2:803934/1‑139 (MQ=255)
                                         atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                     <  2:950424/139‑1 (MQ=255)
                                         atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                     <  1:1009356/139‑1 (MQ=255)
                                         atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                     <  1:803934/139‑1 (MQ=255)
                                              aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                     >  1:583113/1‑134 (MQ=255)
                                              aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                     <  2:583113/134‑1 (MQ=255)
                                              aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGtt                                                >  1:610421/1‑139 (MQ=255)
                                                       gACATAACCCCTTGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                     <  1:1041996/125‑1 (MQ=255)
                                                       gACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                     >  2:1041996/1‑125 (MQ=255)
                                                          aTAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATgg                                                                   <  1:296475/108‑1 (MQ=255)
                                                          aTAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATgg                                                                   >  2:296475/1‑108 (MQ=255)
                                                                 cATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGa                                                               <  2:55472/105‑1 (MQ=255)
                                                                 cATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGa                                                               >  1:55472/1‑105 (MQ=255)
                                                                  aTGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcgcg                            >  1:91505/1‑139 (MQ=255)
                                                                      gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTc                        >  2:587983/1‑139 (MQ=255)
                                                                      gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctgctg                                                                                >  1:803644/1‑83 (MQ=255)
                                                                      gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctgctg                                                                                <  2:803644/83‑1 (MQ=255)
                                                                      gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGAGGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTc                        >  1:791592/1‑139 (MQ=255)
                                                                       tGCAATAGAAAATTTCGACGCCCATACGCCCATGATGTAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCg                       >  2:759482/1‑139 (MQ=255)
                                                                       tGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                     <  2:681425/109‑1 (MQ=255)
                                                                       tGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCg                       >  2:955443/1‑139 (MQ=255)
                                                                       tGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                     >  1:681425/1‑109 (MQ=255)
                                                                        gCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGCGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGc                      <  1:703806/139‑1 (MQ=255)
                                                                         cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCa                     >  2:15437/1‑139 (MQ=255)
                                                                               aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAActgct               <  2:610421/139‑1 (MQ=255)
                                                                                            ccatacgcccatGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTg  <  1:587983/139‑1 (MQ=255)
                                                                                                                     |                                                                                                                    
ATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCACTG  >  NC_000913/2857027‑2857257

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