Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913207,4820TG54.2% 7.5 / 40.0 49V786G (GTA→GGA) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (13/9);  new base G (26/0);  total (39/10)
Fisher's exact test for biased strand distribution p-value = 2.97e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATC  >  NC_000913/207407‑207558
                                                                           |                                                                            
ctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGt                                                                              <  2:4189258/76‑1 (MQ=255)
 tGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTa                                                                             <  2:4820892/76‑1 (MQ=255)
 tGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTa                                                                             <  1:5607797/76‑1 (MQ=255)
  gCGGCCTATGCTTTGGTGTCATATCAACCGTTATGGCTGAAAGCGCGCTCTCCTGCGGAGTTCATGGCGGCGGTaa                                                                            <  2:3848169/76‑1 (MQ=255)
   cGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGaaa                                                                           >  1:2657676/1‑75 (MQ=255)
    ggCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAAta                                                                          <  2:3399247/76‑2 (MQ=255)
    ggCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATg                                                                          <  1:2313338/76‑1 (MQ=255)
       cTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAAAGaac                                                                       >  1:1529960/1‑74 (MQ=255)
          tGCTTGGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAccgcc                                                                    <  2:2166435/76‑1 (MQ=255)
           gCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAccgccg                                                                   <  2:2087855/76‑1 (MQ=255)
            cTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAAAGAccccccc                                                                  >  2:5151605/1‑74 (MQ=255)
            cTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAAAGAACGCCGa                                                                  >  1:5036703/1‑76 (MQ=255)
              ttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGGATGACCCCCGcta                                                                >  2:2761252/1‑76 (MQ=255)
                      atatCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAAAGACCGCCCCTATTGgcaac                                                        >  2:2859117/1‑76 (MQ=255)
                         aaaaaCGTTGGGGCGAAAACCCCACTACCGTGGGGAGTTGAGGGGGGGGGAAACGACCGCCGATATGGACAACAcc                                                     <  1:536518/74‑1 (MQ=255)
                          cAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGGATGGCCGCCGGTATGGGCAACACCg                                                    >  2:3733179/1‑76 (MQ=255)
                           aaaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGGACGCCCGTATGGACAACAcccg                                                   >  1:5815024/1‑74 (MQ=255)
                           aaaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAAAGACCCCCGCGATGGGCCACAACgc                                                   >  2:3921893/1‑75 (MQ=255)
                            aaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAAAGACCGCCGCTAAGGACAACACCcag                                                  >  1:4942209/1‑76 (MQ=255)
                             aCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGACCGCCCATATTGGCAACACCgaga                                                 >  1:4536312/1‑76 (MQ=255)
                               gTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGACCGCCGGTATGGGCAACACCGGGAAg                                               >  1:2965759/1‑76 (MQ=255)
                                ttATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGACAGACCCCCCGTATGGACACAACCGCGACgg                                              >  2:5361863/1‑76 (MQ=255)
                                 tATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGAACGCCGATATGGGCAACACCCAGAAggg                                             >  1:3866524/1‑75 (MQ=255)
                                  ggggCAGAACGCCCACTCTCCGGGGGAGTTTAGGGCGGGGATAATGCCCGCCGATATGGACAACACCGAGAAggtg                                            <  1:3702101/74‑1 (MQ=255)
                                  aTGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGGATGAACGCCGGTATTGGCAACAACGaggagggg                                            >  1:372922/1‑74 (MQ=255)
                                   tGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGACCGCCGATATTGACAACAACGAGGAggggg                                           >  1:1034954/1‑76 (MQ=255)
                                     gCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGACCGCCGGTATGGGCAACACCGAGAAggggggg                                         >  1:1552458/1‑74 (MQ=255)
                                     gCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAAAGACCGCCGATATGGACAACACCCAGAAggggggg                                         >  1:3779707/1‑74 (MQ=255)
                                       tGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGAACGCCGGTATGGACAACACCGAGAAGgggggggg                                       >  1:5017981/1‑76 (MQ=255)
                                        gAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGACCGCCGATATGGACAACACCGAGGAGggggggggg                                      >  1:3098679/1‑75 (MQ=255)
                                        gAAAGCGCACTATCCTGCGGAGTTTATGGCGGAGGGAATGACCGCCGATAAGGACAACAACGAGAAGggggggggg                                      >  1:3148312/1‑75 (MQ=255)
                                          aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGACCGCCGGTATGGGCCACACCGAGGAGGGGGGGGGGCt                                    >  2:781696/1‑76 (MQ=255)
                                          aaGCGCACTATCCTGCGGAGTTTATGGCGGCGGGAATGACCGCCGATATGGACAAACACGAGAAGggggggggggt                                    >  1:1390271/1‑73 (MQ=255)
                                           aGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGGTATGGGCAACACCGCGGAGGGGGGGGGGCtg                                   >  2:4183574/1‑76 (MQ=255)
                                            gcgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCCCCCGTATGGACAACCCCCAGAAGGGGGGGGGGCtgg                                  >  2:3947887/1‑76 (MQ=255)
                                             cgcACTATCCTGCGGAGTTTATGGCGGCGGGAATGACCCCCGATATGGACAACCCCCCGAAGGGGGGGGGGCtggg                                 >  2:5510460/1‑75 (MQ=255)
                                              gcACTATCCTGCGGAGTTTATGGCGGCGGGAAAGACCGCCGGTATGGGCCACACCGGGGAGGGGGGGGGGCtgggg                                >  2:4766759/1‑74 (MQ=255)
                                               cACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGCGAAGGGGGGGGGGCtggggg                               >  1:4980571/1‑76 (MQ=255)
                                               cACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGGAGGGGGGGGGGCtggtgg                               >  1:1539591/1‑76 (MQ=255)
                                                    tCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCCAGAAGGGGGTGGGGCTGGGGGGTGAg                          >  2:3959340/1‑76 (MQ=255)
                                                      cTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGGAGGGGGGGGGGCTGGGGggtgggtg                        >  1:5381152/1‑76 (MQ=255)
                                                          ggAGTTTATGGCGGCGGTAATGACCGCCGGTATGGGCAACACCCGGAAGGGGGTGGGGGTGGGGGGTGAGGGCTgg                    >  1:4996168/1‑76 (MQ=255)
                                                                 aTGGCGGCGGGAATGACCGCCGATATGGACAACACCGAGAAGGTGGGGGGTCTGGTGGGTGAGTGCTGGCGGGTgg             >  1:4078930/1‑76 (MQ=255)
                                                                  tGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGGGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg            >  1:5198910/1‑76 (MQ=255)
                                                                      gggggTAACGCCCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTg        <  1:526081/73‑1 (MQ=255)
                                                                      ggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGGGTGCTGGCGGGTGGGGGTg        >  1:379217/1‑76 (MQ=255)
                                                                      ggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTg        >  1:2271156/1‑76 (MQ=255)
                                                                      ggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGGGGGGCTGGGGGATGAGTGCTGGCGGGGGGGGGTg        >  1:546257/1‑76 (MQ=255)
                                                                       gcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGa       >  1:2437483/1‑76 (MQ=255)
                                                                          gTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGaaaa    >  1:5797082/1‑76 (MQ=255)
                                                                           tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAACAt   >  2:325887/1‑76 (MQ=255)
                                                                           tAATGACCGCCGATATGAAAAAAACCGAGAAGGGGGGGGGCCGGGGGAAGGAGGGC‑GGCGGATGGGGGTGgaaaac  >  1:99454/1‑74 (MQ=255)
                                                                           |                                                                            
CTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATC  >  NC_000913/207407‑207558

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 4 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑