Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913566,3320TG42.0% 15.6 / 81.7 50pseudogene (45/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base G (12/9);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGT  >  NC_000913/566199‑566455
                                                                                                                                     |                                                                                                                           
tttCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTc                                                                                                                        >  1:404379/1‑139 (MQ=255)
 ttCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCg                                                                                                                       >  2:14119/1‑139 (MQ=255)
        cGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccac                                                                                                                >  1:513314/1‑139 (MQ=255)
         gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAaccacc                                                                                                               >  1:507988/1‑139 (MQ=255)
         gCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAaccacc                                                                                                               >  2:268448/1‑139 (MQ=255)
             aaTCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGc                                                                                                           >  2:764497/1‑139 (MQ=255)
                        aaTGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCAt                                                                                                >  2:555924/1‑139 (MQ=255)
                         aTGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATg                                                                                               >  1:784450/1‑139 (MQ=255)
                         aTGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATg                                                                                               <  2:214054/139‑1 (MQ=255)
                          tGCAGGCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGa                                                                                              >  1:850312/1‑139 (MQ=255)
                               gCCTTCACGCTTCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATcc                                                                                         >  2:608688/1‑139 (MQ=255)
                                   tCACGCTTCATGCGCGCGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCCGAATTTCATGAAATCCCggg                                                                                     <  2:685194/139‑1 (MQ=255)
                                         ttCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaa                                                                               <  1:876507/139‑2 (MQ=255)
                                         ttCATACGCGGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaa                                                                               <  2:780212/139‑2 (MQ=255)
                                              acgcgGGTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                                                                          <  1:14119/138‑7 (MQ=255)
                                                   ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAGGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                                                                     >  1:799601/1‑128 (MQ=255)
                                                   ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                                                                     <  1:555924/139‑12 (MQ=255)
                                                   ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                                                                     <  2:849687/139‑12 (MQ=255)
                                                   ggTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacgggcattt                                                                     >  2:914544/1‑128 (MQ=255)
                                                         aGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATTCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGGGGTACGAGTATTGCCGGAAGc                                                               <  2:507988/139‑1 (MQ=255)
                                                           ttGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGt                                                             <  2:823021/139‑1 (MQ=255)
                                                              gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCAGCTGGGCCATGTAAGCCGATTCTATTGCCGCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGCATTGCCGGAAGCGTGGc                                                          <  1:578735/139‑1 (MQ=255)
                                                              gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                          <  2:994323/139‑1 (MQ=255)
                                                              gCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGc                                                          >  1:994323/1‑139 (MQ=255)
                                                                  aGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGt                                                      >  2:572887/1‑139 (MQ=255)
                                                                             tCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                           <  1:927325/139‑1 (MQ=255)
                                                                             tCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                           <  1:91668/139‑1 (MQ=255)
                                                                                         aCCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGt                                       <  1:985284/131‑1 (MQ=255)
                                                                                         aCCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGt                                       >  2:985284/1‑131 (MQ=255)
                                                                                                         gCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa                                                                                    <  2:905320/70‑1 (MQ=255)
                                                                                                         gCACCTGGGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGa                                                                                    >  1:905320/1‑70 (MQ=255)
                                                                                                             cTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTgg                                                           <  1:912118/91‑1 (MQ=255)
                                                                                                             cTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTgg                                                           >  2:912118/1‑91 (MQ=255)
                                                                                                               gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCgg                                                <  1:189618/100‑1 (MQ=255)
                                                                                                               gggCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCgg                                                >  2:189618/1‑100 (MQ=255)
                                                                                                               gggCCATGTAAGCCGA‑TTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                            >  2:658675/1‑103 (MQ=255)
                                                                                                               gggCCATGTAAGCCGA‑TTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGc                                            <  1:658675/103‑1 (MQ=255)
                                                                                                                    aTGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                                                                          <  2:612413/69‑7 (MQ=255)
                                                                                                                    aTGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGTaaacggg                                                                          >  1:612413/1‑63 (MQ=255)
                                                                                                                    aTGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtgg       <  1:956428/136‑1 (MQ=255)
                                                                                                                    aTGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtgg       >  2:956428/1‑136 (MQ=255)
                                                                                                                     tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCAt                                                                                                <  1:530314/46‑1 (MQ=255)
                                                                                                                     tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCAt                                                                                                >  2:530314/1‑46 (MQ=255)
                                                                                                                     tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                  <  1:446829/124‑1 (MQ=255)
                                                                                                                     tGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAAcg                  >  2:446829/1‑124 (MQ=255)
                                                                                                                      gTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACtggttggt  >  1:397328/1‑139 (MQ=255)
                                                                                                                           ccGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                           <  2:721504/93‑1 (MQ=255)
                                                                                                                           ccGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCa                                           >  1:721504/1‑93 (MQ=255)
                                                                                                                                 ttATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGt                                                                                 <  2:953632/49‑1 (MQ=255)
                                                                                                                                 ttATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCGGGAGGt                                                                                 >  1:953632/1‑49 (MQ=255)
                                                                                                                                     |                                                                                                                           
TTTCATTCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATGCGCGGGTCATAGTTGGCAAAGTACCAGGCATCTTTTCGCGTCACCCACATGCTGTACTGCACCTGGGCCATGTAAGCCGATTTTATTGCCTCGAAACCACCGAGCCGGAATTTCATGAAATCCCGGGAGGTACGAGTATTGCCGGAAGCGTGGCCTGTATCCGGATGCAGAGTCTTATCCGTGGAAATCGAACGCGCATTACTGGTTGGT  >  NC_000913/566199‑566455

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