Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,4640GT73.9% ‑6.1 / 12.3 23V136V (GTG→GTTnohBbacteriophage DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (9/8);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCG  >  NC_000913/580374‑580548
                                                                                          |                                                                                    
tggaGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAgcg                                                                                 <  1:589294/93‑1 (MQ=1)
 gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGCGAATTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAgcgg                                                                                >  2:9387174/4‑96 (MQ=0)
 gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAgcgg                                                                                <  1:1430878/93‑1 (MQ=0)
 gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAgcgg                                                                                >  2:7623164/4‑96 (MQ=0)
 gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAgcgg                                                                                >  1:4975064/4‑96 (MQ=0)
 gaggtggtgGAAACCGCATTCTGCACGTTCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAgcgg                                                                                <  1:10446144/93‑1 (MQ=0)
          gAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGTTCGCCCTGTCGGTGCAGCGGCTTTTcctg                                                                      <  1:7609413/96‑3 (MQ=1)
                      agtccTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAAc                                                           <  1:8910366/92‑1 (MQ=1)
                        cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCg                                                         <  1:7727629/92‑1 (MQ=0)
                        cacgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCg                                                         >  1:8899563/5‑96 (MQ=0)
                         acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGa                                                        <  1:7903363/93‑1 (MQ=0)
                         acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGa                                                        <  2:7873803/93‑1 (MQ=0)
                         acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGa                                                        >  2:5235007/4‑96 (MQ=0)
                         acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGa                                                        <  2:4970517/93‑1 (MQ=0)
                         acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGa                                                        <  2:199353/93‑1 (MQ=0)
                         acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGa                                                        >  1:9658544/4‑96 (MQ=0)
                         acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGa                                                        >  1:37677/4‑96 (MQ=0)
                         acgttCGTGCTGTCGCGGATGGCCGGAGAAGTAGCCAGCATTCTTGATGGAGTTCCTCTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGa                                                        <  1:34853/93‑1 (MQ=0)
                                         ggctCGTAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGa                                        <  1:6176300/93‑1 (MQ=0)
                                                agttGAAATTGCCACTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGa                                 <  1:1176939/93‑1 (MQ=0)
                                                                            atgccccTGTCGGTGCAGAGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGcc     >  1:8389948/4‑96 (MQ=0)
                                                                               cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCtgcg  >  2:14306642/4‑96 (MQ=0)
                                                                               cctcTGTCGGTTCAGCGGCGCTTCCCGGAGCTGGAAAACCGACATATTGATTTCCTCAAGAAGGACATCATTAAAGCCATGAACAAAGCAGCtgcg  >  1:5898286/4‑96 (MQ=0)
                                                                                          |                                                                                    
TGAAGTGGTGGAAACCGCATTCTGTACTTTCGTGCTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGATCCCCCTGTCGGTGCAGCGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCG  >  NC_000913/580374‑580548

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