Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,857,1420AG42.1% 19.7 / 64.6 57Y17C (TAT→TGT) mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base A (15/18);  new base G (14/10);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCA  >  NC_000913/2857024‑2857254
                                                                                                                       |                                                                                                                  
cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGc                                                                                              >  2:1130038/1‑139 (MQ=255)
      tgATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCg                                                                                        <  2:1632921/139‑1 (MQ=255)
                 aaaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctgctg                                                                             <  1:1130038/139‑1 (MQ=255)
                         aaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGt                                                                            <  2:423260/132‑1 (MQ=255)
                         aaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGt                                                                            >  1:423260/1‑132 (MQ=255)
                              tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTa                                                                        <  2:119228/131‑1 (MQ=255)
                              tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTa                                                                        >  1:119228/1‑131 (MQ=255)
                               cacaCCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctg                                                                                >  2:1578966/1‑122 (MQ=255)
                               cacaCCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctg                                                                                <  1:1578966/122‑1 (MQ=255)
                               cacaCCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATggg                                                               >  1:1465218/1‑139 (MQ=255)
                                      cATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGAtttt                                                        >  2:1218228/1‑139 (MQ=255)
                                       aTTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGAttttt                                                       <  1:1218228/139‑1 (MQ=255)
                                        tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCctgct                                                                              <  2:1480855/115‑1 (MQ=255)
                                        tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCctgct                                                                              >  1:1480855/1‑115 (MQ=255)
                                           aaTATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGa                                                   >  1:822586/1‑139 (MQ=255)
                                            atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAGATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  <  1:1589106/139‑1 (MQ=255)
                                            atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  >  1:1471823/1‑139 (MQ=255)
                                            atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  <  2:865088/139‑1 (MQ=255)
                                            atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  <  2:106405/139‑1 (MQ=255)
                                            atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  >  1:865088/1‑139 (MQ=255)
                                            atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  <  2:1465218/139‑1 (MQ=255)
                                            atatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa                                                  <  2:822586/139‑1 (MQ=255)
                                                cAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCct                                                                                 <  1:1496158/104‑1 (MQ=255)
                                                cAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGt                                              >  1:653229/1‑139 (MQ=255)
                                                cAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCct                                                                                 >  2:1496158/1‑104 (MQ=255)
                                                   ggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCctg                                                                                <  2:489060/102‑1 (MQ=255)
                                                   ggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGa                                                                   >  1:1420269/1‑115 (MQ=255)
                                                   ggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctg                                                                                >  1:489060/1‑102 (MQ=255)
                                                   ggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGa                                                                   <  2:1420269/115‑1 (MQ=255)
                                                              tAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCa                                <  2:1471823/139‑1 (MQ=255)
                                                                    cATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcgc                          >  1:1648269/1‑139 (MQ=255)
                                                                     aTGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGc                                                                                              >  1:1264953/1‑70 (MQ=255)
                                                                     aTGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGc                                                                                              <  2:1264953/70‑1 (MQ=255)
                                                                       gAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctg                                                                                >  2:1610170/1‑82 (MQ=255)
                                                                       gAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCctg                                                                                <  1:1610170/82‑1 (MQ=255)
                                                                         gTGCAATAGTAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTc                     >  2:1572652/1‑139 (MQ=255)
                                                                         gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCg                                                                                        <  2:60468/72‑1 (MQ=255)
                                                                         gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCg                                                                                        >  1:60468/1‑72 (MQ=255)
                                                                          tGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCg                    >  1:986518/1‑139 (MQ=255)
                                                                          tGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                                 >  2:1283495/1‑110 (MQ=255)
                                                                          tGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTCATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc                                                 <  1:1283495/110‑1 (MQ=255)
                                                                            cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCctgct                                                                              <  1:697829/79‑1 (MQ=255)
                                                                            cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATggg                                                               <  1:1532501/94‑1 (MQ=255)
                                                                            cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATggg                                                               >  2:1532501/1‑94 (MQ=255)
                                                                            cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCa                  >  2:387713/1‑139 (MQ=255)
                                                                            cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCctgct                                                                              >  2:697829/1‑79 (MQ=255)
                                                                            cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATc                                                                                   >  2:953194/1‑74 (MQ=255)
                                                                            cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATc                                                                                   <  1:953194/74‑1 (MQ=255)
                                                                             aaTAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCaa                 <  2:1292274/139‑1 (MQ=255)
                                                                                  aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcg                           <  2:4079/124‑1 (MQ=255)
                                                                                  aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcg                           >  1:4079/1‑124 (MQ=255)
                                                                                  aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAActgct            >  2:667028/1‑139 (MQ=255)
                                                                                    aaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTgg          <  2:1648269/139‑1 (MQ=255)
                                                                                    aaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTgg          >  1:1142468/1‑139 (MQ=255)
                                                                                         cGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGa                                                            <  2:1263512/84‑1 (MQ=255)
                                                                                         cGACGCCCATACGCCCATGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGa                                                            >  1:1263512/1‑84 (MQ=255)
                                                                                            cgcccatacgcccatGATGCAGCAGTGATAGTC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCa  <  2:653229/139‑1 (MQ=255)
                                                                                                                       |                                                                                                                  
CTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTGGATATTTCA  >  NC_000913/2857024‑2857254

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