Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,580,0950AC22.4% 108.3 / 12.9 49G456G (GGT→GGGhsdStype I restriction enzyme EcoKI specificity protein
Reads supporting (aligned to +/- strand):  ref base A (29/9);  new base C (0/11);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 5.76e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGA  >  NC_000913/4579958‑4580229
                                                                                                                                         |                                                                                                                                      
aaaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAcc                                                                                                                                      <  1:709933/140‑1 (MQ=255)
 aaTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAccc                                                                                                                                     >  1:537596/1‑140 (MQ=255)
  aTAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTAcccc                                                                                                                                    <  1:158353/140‑1 (MQ=255)
           aaTGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGCCCTTTTTACCCCCGCTAGCTg                                                                                                                           >  1:583751/1‑140 (MQ=255)
                gCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGt                                                                                                                      >  1:662948/1‑140 (MQ=255)
                  aaTTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTc                                                                                                                    >  1:145186/1‑140 (MQ=255)
                  aaTTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTc                                                                                                                    >  1:585852/1‑140 (MQ=255)
                       cTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCtt                                                                                                               <  1:688949/140‑1 (MQ=255)
                               aTTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGAttttt                                                                                                       >  1:512423/1‑140 (MQ=255)
                                         atTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAAcg                                                                                             >  2:329011/1‑140 (MQ=255)
                                         atTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAAcg                                                                                             >  2:588336/1‑140 (MQ=255)
                                         atTTATGGGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGGGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAAcg                                                                                             >  2:368286/1‑140 (MQ=255)
                                            tATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcgg                                                                                          <  1:99516/140‑1 (MQ=255)
                                              tgtgGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggcg                                                                                        <  2:145186/140‑1 (MQ=255)
                                              tgtgGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACgcggcg                                                                                        <  2:585852/140‑1 (MQ=255)
                                                   gTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCt                                                                                   >  1:84193/1‑140 (MQ=255)
                                                   gTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCt                                                                                   >  1:741847/1‑140 (MQ=255)
                                                   gTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCt                                                                                   >  2:653793/1‑140 (MQ=255)
                                                             gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg                                                                         >  1:3696/1‑140 (MQ=255)
                                                             gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg                                                                         <  2:203998/140‑1 (MQ=255)
                                                             gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg                                                                         >  1:405197/1‑140 (MQ=255)
                                                             gCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCg                                                                         >  1:376290/1‑140 (MQ=255)
                                                                               cgcACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGttt                                                       >  1:207070/1‑140 (MQ=255)
                                                                                            tGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACt                                          <  2:405197/140‑1 (MQ=255)
                                                                                              cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  2:565735/1‑140 (MQ=255)
                                                                                              cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  2:461313/1‑140 (MQ=255)
                                                                                              cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  1:613698/1‑140 (MQ=255)
                                                                                              cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  1:614600/1‑140 (MQ=255)
                                                                                              cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  1:18910/1‑140 (MQ=255)
                                                                                              cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  1:463353/1‑140 (MQ=255)
                                                                                              cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  1:733262/1‑140 (MQ=255)
                                                                                              cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  2:154975/1‑140 (MQ=255)
                                                                                              cgcCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTAATCAAATCCGGGTTTTCGGCCCGCCACTgg                                        >  2:605974/1‑140 (MQ=255)
                                                                                                ccAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGc                                      >  1:400483/1‑140 (MQ=255)
                                                                                                caaaaaaaTAATGTTCAGGATTTTTTAGGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGTTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCCCCGCTGATCAAATCCGGTTTTTCGGCCCGCCACTGGGc                                      <  2:741847/136‑1 (MQ=255)
                                                                                                  aGAAAATAATGTTCAGGATTTTTTACGTGAGGTTTTTTTCCCCCCGCTAGCTGCGCTTTCAGCTTTGATTTTTCCCAGCAACGCGGCGGCGCTTTTTTCTCCGCTGATAAAATCCGGGTTTTCGGCCCGCCACTGGGCgg                                    <  2:752005/140‑1 (MQ=255)
                                                                                                        taatGTCCAGGTTTTTTTACGGGGGGCTTTTTTCCCCCCGCTAGGTGCGGTTTTTTTTTTGATTTTTTCAGGAACCGCGGGGGCGCTGTTTCCTCCGCTGATCAAACCCGGTTTTTCGGCCCGCCACTGGGCGGTAAGtt                              <  2:376290/140‑1 (MQ=255)
                                                                                                             ttCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTcacca                         >  1:464996/1‑140 (MQ=255)
                                                                                                                aGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACgg                      >  1:234389/1‑140 (MQ=255)
                                                                                                                aGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACgg                      >  2:447784/1‑140 (MQ=255)
                                                                                                                aGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACgg                      >  1:27194/1‑140 (MQ=255)
                                                                                                                 ggtttttttACGTGAGGTTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGa                     <  1:723819/137‑1 (MQ=255)
                                                                                                                 ggATTTTTTACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGa                     <  1:447784/140‑1 (MQ=255)
                                                                                                                      ttttACGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCt                <  1:274149/140‑1 (MQ=255)
                                                                                                                      ttttACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGGGGCGGCGCTGTTTTCTCCGCTGATCAAATCTGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCt                <  1:263008/140‑1 (MQ=255)
                                                                                                                          aCGTGAGGTTTTTTTCCCCCCGCTAGCTGCGGGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTg            <  2:613698/140‑1 (MQ=255)
                                                                                                                          aCGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCCCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTg            <  1:54221/140‑1 (MQ=255)
                                                                                                                           cGTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGc           <  1:565735/140‑1 (MQ=255)
                                                                                                                            gTGAGGCTTTTTTCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc          <  1:374307/140‑1 (MQ=255)
                                                                                                                            gTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc          <  1:653793/140‑1 (MQ=255)
                                                                                                                            gTGAGCCCTTTTTCCCCCCGCTAGCTGCGCGTTCAGTTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGTTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGcc          <  2:583751/140‑1 (MQ=255)
                                                                                                                                    tttttCCCCCCGCTGCCTGCGCGTTCAGCTTTGTTTTTTTCCAGCAACGCGGCGGCGCTGTTTTCCCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggatgga  <  2:27194/140‑1 (MQ=255)
                                                                                                                                    tttttCCCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAggatgga  <  2:234389/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
AAATAGAATAAAATGGGCAATTTCTGGAATGATTTAAATATATTTATGTGGGTTATGATTGGCGTGAAATAATAAAAAGCGCACCGGAAAGGTGCGCCAGAAAATAATGTTCAGGATTTTTTACGTGAGGCTTTTTTACCCCCGCTAGCTGCGCGTTCAGCTTTGATTTTTTCCAGCAACGCGGCGGCGCTGTTTTCTCCGCTGATCAAATCCGGGTTTTCGGCCCGCCACTGGGCGGTAAGTTCACCACGGAACGCTTTTGCCAGGATGGA  >  NC_000913/4579958‑4580229

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