Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,6050TG54.8% 8.2 / 26.4 31K97Q (AAA→CAA) borDprophage lipoprotein BorD
Reads supporting (aligned to +/- strand):  ref base T (6/8);  new base G (8/9);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCC  >  NC_000913/578380‑578848
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
ggtttGCACATAAAACCAGTCTGTCGTCCAGGGCGGAAAGCCCCCCCCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGAATTATTTCACCCTCAGAGAGAGGCGGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATGAATGAGCAGTGGAGATAGAGTTCCCCATATCGATGGGCAACTAATGCAATTATTGTGAGCAATACACACGCGCTTCCAGc                                                                                                                                                                                                                            <  1:438423/248‑1 (MQ=25)
               ncAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCt                                                                                                                                                                                                             >  1:266656/2‑251 (MQ=255)
                                                             gATAATGCAGCCTCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAACAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAACGTAATAAAACCGAGCAATCCATTTACGAATGTTTGc                                                                                                                                                                      >  1:253125/1‑244 (MQ=255)
                                                                aaTGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACTCGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGtttn                                                                                                                                                              <  1:269890/250‑2 (MQ=255)
                                                                    naGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGttt                                                                                                                                                         >  1:265313/2‑251 (MQ=255)
                                                                                      cAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTCTATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTgcgc                                                                                                                                       >  1:194271/1‑251 (MQ=255)
                                                                                                                  gcTCAGAGAGAGGCTGCGCACTATGCAAAAACAACTGGAAGGCACCCAGCAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGCTGGGCAACTCATGCAATTATTGCGAGCAATACACACGCGCTTNCAGNNNNNtann                                                                                                                                                                                                                    >  2:264548/2‑144 (MQ=255)
                                                                                                                         gagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACTGAGCAATCCACTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTttcttc                                                                                                    <  1:61663/251‑1 (MQ=255)
                                                                                                                            agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGTGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGATTATAAATGCCTAAAGTAATAAAACCGAGCAAGCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAg                                                                                                                        <  1:291991/228‑1 (MQ=255)
                                                                                                                                 gATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTgg                                                                                                                                                                   <  1:157862/180‑1 (MQ=255)
                                                                                                                                 gATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGAATATTAATGAGCAGTGCAGATAGCGCTGACCATATCGATGGGCCCCTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGt                                                                                                                                                                          >  1:180902/1‑173 (MQ=255)
                                                                                                                                           cAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAAGAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGc                                                                                                 <  1:199158/236‑1 (MQ=255)
                                                                                                                                                      ggaaggaaCCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTATTAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAAttttn                                                                                                                         <  1:281137/201‑2 (MQ=255)
                                                                                                                                                      ggaaggaaCCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGANCAANNNNnttnnnnnnnnnn                                                                                                                                                                        >  2:281137/1‑144 (MQ=255)
                                                                                                                                                      ggaaggaaCCCAGAAGTATATTAATGAGCAGTGCAGAGAGATTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAatgatn                                                                       <  1:133068/251‑2 (MQ=255)
                                                                                                                                                           gaaCCCCGCAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAGCTCATGCAATTATTTTGAGCAGTACACACGCGCTTCCAg                                                                                                                                                                                                                             >  1:57073/1‑96 (MQ=255)
                                                                                                                                                            aaCCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGAt                                                                 >  1:246012/1‑251 (MQ=255)
                                                                                                                                                                                   aGTGCAGATAGATCTTACTGTATCGATGGGCAACTCATGCAAGTATTTGGAGCAATACACACGCGCGGCCCGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGTATGTTTGCGGGGTGTCTTTGTTAACAACGGTTTCTAGGCCGCCACAAATTGTAGCTTCACCGACAGTGTTCTGCTGCCCAAATCCAGAAACGAAGAAATGATGGGTGAGGGTTTCCTTTGGTGCTActgctg                                          <  1:321496/251‑1 (MQ=255)
                                                                                                                                                                                   aGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCAGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTActgctg                                          >  1:33578/1‑251 (MQ=255)
                                                                                                                                                                                    gTGCAGATAGAGGTTCCCATATCGAGGGGCAACTCATGCAATTATTGTGATCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATc                                                                                                                                                                                       <  1:64951/109‑1 (MQ=255)
                                                                                                                                                                                        aGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTTCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACg                                                                                  <  1:256156/206‑1 (MQ=255)
                                                                                                                                                                                          aTAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCGTTGGTGCTActgctgccggttt                                   <  1:121044/251‑8 (MQ=255)
                                                                                                                                                                                                  gACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGGCACAAATTTTAGCTGCATCGACAGTTttcttc                                                                                                    <  1:377164/178‑1 (MQ=255)
                                                                                                                                                                                                  gACCATATCGATGGGCAACTCATGCAAGTATGTTGAGCAATACACAcgcg                                                                                                                                                                                                                                    <  1:357152/50‑1 (MQ=255)
                                                                                                                                                                                                    ccATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATc                                                                                                                 >  1:425669/1‑163 (MQ=255)
                                                                                                                                                                                                               ggCAACTCATGCAATTATTTTGAGCAATACACAAGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTttcttc                                                                                                    >  1:127272/1‑165 (MQ=255)
                                                                                                                                                                                                                 naaCTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTc                                                                                          >  1:131624/2‑173 (MQ=255)
                                                                                                                                                                                                                      cATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTCACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGc        >  1:205718/1‑250 (MQ=255)
                                                                                                                                                                                                                      cATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAAGAAAACCGAGCAATCCATTTACGAATGGTGGCTGGGTGTCTGTTTTAACAACATTTTCTGCGCAGCCACAAAGTTTGGCTGCATCGACAGTTTTCTTCTGTCCAATTCCAGAAACGAAGAAAGGATGGGTGATTGTGTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGc        <  1:405053/250‑1 (MQ=255)
                                                                                                                                                                                                                          cAATTATTGTGAGCAATACACCCGCGCTGCCAGCGGAGTATAAATGCCTAAATTAATAAAACCGAGCAAGCCATTTACGAATGTTTGCGGGGTTTCTGTTTTAGCAACATTTTCTGCGCCGCCACAAATGTTGGCTTCATCGACAGTTTTCTTCTGc                                                                                                 <  1:262730/157‑1 (MQ=255)
                                                                                                                                                                                                                           aaTTATTGTGAGCAATACGCACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGGGATGGTTTCCTGTGGTGCTACTGCTGCCGGTTTGTTTTGAACAGTAAACGTCTGTTGAGCACATcc  <  1:202204/251‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
CGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAACAACATTTTCTGCGCCGCCACAAATTTTAGCTGCATCGACAGTTTTCTTCTGCCCAATTCCAGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGTCTGTTTGTTTTGAACAGTAAATGTCTGTTGAGCACATCC  >  NC_000913/578380‑578848

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