Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091364,6040CT37.5% 7.5 / 31.0 34E393K (GAA→AAA) polBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base C (11/9);  new base T (5/7);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCA  >  NC_000913/64368‑64840
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
agtaacaATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCTGACGGGATCAATCAGAAAGGTGCGGATGATCGTCGGGTACAGGCTTTTATAGTCCAGCACCAGCACTT‑AATCATAAAGCCCTTGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGgcg                                                                                                                                                                                                                                <  2:236861/247‑1 (MQ=255)
aGTCACAATCTCCGGCAGGCACTGGTTTTCTCGCGAGCACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCAACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGATTTTATAGTCCAGCACTAGCCATG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGTGCTTCGCCGAGTATAGgcg                                                                                                                                                                                                                                >  2:223211/1‑251 (MQ=255)
     caATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTTGCCGAGATTAGGCGCGACa                                                                                                                                                                                                                           <  2:216718/251‑1 (MQ=255)
          tcCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAAcc                                                                                                                                                                                                                      >  1:268715/1‑251 (MQ=255)
                  gcaACGTTTTGCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTTGCCGAGATTAGGCGCGACATAACCAGCGCGAt                                                                                                                                                                                                              <  2:233687/251‑1 (MQ=255)
                       gTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTTGCCGAGATTAGGCGCGACATAACCAGCGCGATGCAtt                                                                                                                                                                                                         <  2:223448/251‑1 (MQ=255)
                             ctagcaagcaCCAGGCATCGAGAACACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGTATCGATCAGAACGGCGCGGATGCTCGCCGGGTACAGGATTTTATAGGCCAGCACCAGCACTG‑AATCATCAAGCCATGGCTGTGAATCCATCAAGTAGCCGACAGGGCTGGCGTTCGGCCGCACTTCGCAGAGATTAGGCGCGACATAACCAGCGCGATGGAttacaggg                                                                                                                                                                                                   >  2:63420/10‑245 (MQ=255)
                                    gaACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCATCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCTTCACGTAGCCGCCAGGGCGGGCGTGCGGCGTCATTTCGCCGAGATTAGGTGCGACATAACCAGCGCGATGCATTCGCGGAAAATaga                                                                                                                                                                                            <  1:154852/251‑1 (MQ=255)
                                    gaACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCAATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTTGCCGAGATTAGGCGCGACCTAACCAGCGCGCTGCATTCGCgggatcgga                                                                                                                                                                                             >  1:28856/1‑243 (MQ=255)
                                    gaACCAGGCATCGAGAAAACCCTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGGCCAGCACCAGCACTG‑AATCATAAAGCCCGGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACGTTGCCGATATTAGGCGCGACATAACCAGCGCGATGCATTcgcg                                                                                                                                                                                                     <  2:28856/242‑1 (MQ=255)
                                       ccAGGCATCGAGAAAAAATTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCATTCCACCAGACCGACGGGCTCAATCAGAAAGGGGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAATCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCCGCGAGATGCAGTCGCGTAAACTag                                                                                                                                                                                             >  2:4829/1‑247 (MQ=255)
                                                     aaCCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCCGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGCGATGACCAAATgccgcc                                                                                                                                                                           >  2:177691/1‑251 (MQ=255)
                                                                              aGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGCCATAACCAgcgc                                                                                                                                                                                                                 >  2:241781/1‑188 (MQ=255)
                                                                              aGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAgcgc                                                                                                                                                                                                                 <  1:241781/188‑1 (MQ=255)
                                                                                                   caGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACGG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGCCCAAAGGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGccc                                                                                                                             <  1:120877/251‑1 (MQ=255)
                                                                                                       ccGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGaa                                                                                                                                                                                                  >  1:73624/1‑178 (MQ=255)
                                                                                                       ccGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTAGAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGaa                                                                                                                                                                                                  <  2:73624/178‑1 (MQ=255)
                                                                                                        cGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTTGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCg                                                                                                                        <  1:38645/251‑1 (MQ=255)
                                                                                                         gACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCgg                                                                                                                       >  2:152089/1‑250 (MQ=255)
                                                                                                             ggATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTTAGGCGGCACTTTGCCGCGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCAACGTTGCCCGCTCgggtc                                                                                                                    >  2:225280/1‑246 (MQ=255)
                                                                                                             ggATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTa                                                                                                                    >  1:153029/1‑250 (MQ=255)
                                                                                                               atcaatcaAAAAGGTGGGGTTGATCGACGGGTACAGGTGTTTAAAGCCCAGCCCAAGCCCTG‑AATCTGAACGCCCTGGCCGTGAACCCATCACGTAGCCGCCAGGGCTGGCGGGAGGCGGCACTGCGCCGAGATAAGGCGCGACATACCCAGCGCGATGCACTCGCTGAAAATAGAGATGCCCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTGCACCGTTGCCCGTTCGAGTaaaa                                                                                                                 <  2:60664/251‑1 (MQ=255)
                                                                                                                        aaaGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGa                                                                                                        <  1:152089/251‑1 (MQ=255)
                                                                                                                        aaaGGGGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCCCCAGC‑CTGGAATCATAAAGCCCGGGCGGTGATCCCACCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTGCGCCGAGATTAGGCGCGACATACCCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGa                                                                                                        <  2:153029/251‑1 (MQ=255)
                                                                                                                             ggCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGGGGCACTTTGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAAGAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCa                                                                                                   <  1:26704/250‑1 (MQ=255)
                                                                                                                                       ccttttccgatcttttttATAGTCCAGCACCAGCCCGG‑AACCAAAATGCCTGGGCGGTGAACCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCTCTTGGCCGAGATTAGGCGCGCAATACCCAGCGCGATGCATCCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGa                                                                                         <  2:170354/237‑1 (MQ=255)
                                                                                                                                        gACGTGTAAAGGCTTGTATACTCCCGCACCCGCACTA‑AATCAAAAAGCCCTGGCCGTGAATCCAGCACGAAGCCGCCAGGGCTTGCGTGGGGCGGCTCTTGGCCGAGATTAGGCGCGACATAACCAGCGCGATGCTTTCGCGGAAAATAGAGATGACCAAGTGCCGCCCCCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGaa                                                                                        <  2:88067/251‑1 (MQ=255)
                                                                                                                                                    cTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTTGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGa                                                                                         >  1:170354/1‑238 (MQ=255)
                                                                                                                                                         ataGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTTGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCt                                                                       >  1:229197/1‑251 (MQ=255)
                                                                                                                                                            gTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGc                                                                    >  1:244723/1‑251 (MQ=255)
                                                                                                                                                                           tG‑AATCATAAAGCACTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGTAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTtat                                                     >  2:282953/1‑251 (MQ=255)
                                                                                                                                                                             gaagaaTAAAGCCCTGGCCGTGAATCCCTCACGTGGCCGCCAGGGCTGGCGGGCGGCGGCTCTTTGCCGAGATTAGGCGCGCCATAACCAGCGCGATGCATTCGGGGAAAATGGAGATGCCCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTtat                                                     <  2:43456/245‑1 (MQ=255)
                                                                                                                                                                                        cccTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTTGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAgcg                                         >  1:143062/1‑251 (MQ=255)
                                                                                                                                                                                                                               gTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTACGTCCa  <  1:91634/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                            
AGTCACAATCTCCGGCAGGCAATGTTTTTCTCGCGAGAACCAGGCATCGAGAAAACCTTCGGTACTGTGCTCTGGATCAGGCTGCGCCATGCCTTCCACCAGCCCGACGGGATCAATCAGAAAGGTGCGGATGATCGACGGGTACAGGCTTTTATAGTCCAGCACCAGCACTG‑AATCATAAAGCCCTGGCCGTGAATCCATCACGTAGCCGCCAGGGCTGGCGTGCGGCGGCACTTCGCCGAGATTAGGCGCGACATAACCAGCGCGATGCATTCGCGGAAAATAGAGATGACCAAATGCCGCCACCGAACCGCCGTGTCGGTCCACCGGCAGGCCGTTCACCGTTGCCCGTTCGAGTAAAAATGGCATGATTTCAGTTTTGTGGAAGATCTGCGTCACCAGCTCGCAATCTTTCAGGTTATAAGTTGCCAGCGCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCA  >  NC_000913/64368‑64840

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