Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479696,5260TA23.8% 40.0 / 11.2 23intergenic (+4/‑125)h16_A0653/h16_A0654conserved hypothetical protein/Hemolysin‑coregulated protein
Reads supporting (aligned to +/- strand):  ref base T (7/9);  new base A (2/3);  total (10/13)
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.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGAGGAACTGGAAGGCAACCCCGGCTACTACAGCGCCAAGTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTGGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGCTCTTGTCTGCAGTACCTGTCTTAACGTCAAGGAGTTTTGCAATGG  >  AM260479/696393‑696654
                                                                                                                                     |                                                                                                                                
cgAGGAACTGGAAGGCAACCCCGGCTACTACAGCGCCAAGTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTGGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGt                                                                                                                                  >  2:354844/1‑134 (MQ=255)
 gAGGAACTGGAAGGCAACCCCGGCTACTACAGCGCCAAGTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTGGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGt                                                                                                                                  <  1:310070/133‑1 (MQ=255)
            aGGCAACCCCGGCTACTACAGCGCCAAGTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTGGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGACGCTCAGGCGtt                                                                                                                      <  2:730540/134‑1 (MQ=255)
                    ccGGCTACTACAGCGCCAAGTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTGGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTa                                                                                                              <  1:354844/134‑1 (MQ=255)
                                     aaGTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTTGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGt                                                                                                      >  1:242568/1‑125 (MQ=255)
                                     aaGTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTTGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGt                                                                                                      <  2:242568/125‑1 (MQ=255)
                                       gTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTGGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGACGCTCAGGCg                                                                                                                        >  1:245083/1‑105 (MQ=255)
                                       gTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTGGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGACGCTCAGGCg                                                                                                                        <  2:245083/105‑1 (MQ=255)
                                                                                     cTGCGGCTTGTGGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGctc                                             >  1:273771/1‑134 (MQ=255)
                                                                                               tgGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGCTCTTGTCTGCa                                    <  1:678167/133‑1 (MQ=255)
                                                                                               tgGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGCTCTTGTCTGCa                                    <  2:273771/133‑1 (MQ=255)
                                                                                                 gcgcgcCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCt                                                    >  1:864696/1‑115 (MQ=255)
                                                                                                 gcgcgcCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCt                                                    <  2:864696/115‑1 (MQ=255)
                                                                                                                   aaGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTgg                                                                   >  2:410807/1‑82 (MQ=255)
                                                                                                                   aaGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTgg                                                                   <  1:410807/82‑1 (MQ=255)
                                                                                                                   aaGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGCTCTTGTCTGCAGTACCTGTCTTAACGTCAAg                >  2:378869/1‑133 (MQ=255)
                                                                                                                   aaGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGCTCTTGTCTGCAGTACCTGTCTTAACGTCAAg                >  2:826028/1‑133 (MQ=255)
                                                                                                                       aCGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGCTCTTGTCTGCAGTACCTGTCTTAACGTCAAGGAGtt           <  1:701605/134‑1 (MQ=255)
                                                                                                                        cgccgcCTGAGCGCCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGct                                              >  2:1085027/1‑98 (MQ=255)
                                                                                                                        cgccgcCTGAGCGCCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGct                                              <  1:1085027/98‑1 (MQ=255)
                                                                                                                               tGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATcgc                                                                                    >  2:1161186/1‑53 (MQ=255)
                                                                                                                               tGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATcgc                                                                                    <  1:1161186/53‑1 (MQ=255)
                                                                                                                               tGAGCGACGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATcgc                                                                                    <  1:461027/53‑1 (MQ=255)
                                                                                                                               tGAGCGACGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATcgc                                                                                    >  2:461027/1‑53 (MQ=255)
                                                                                                                                 aGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGCTCTTGTCTGCAGTACCTGTCTTAACGTCAAGGAGTTTTGCAATgg  <  1:826028/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                
CGAGGAACTGGAAGGCAACCCCGGCTACTACAGCGCCAAGTTCTTCCTGCGTCCGCACTTCCAGCTCGAGGGCCTGACCGTTCCGCTGCGGCTTGTGGCGCGCCTGCCGTCGGTCAAGGACGCCGCCTGAGCGTCGCTCAGGCGTTGAAACGTACCGGTAGTTCCACGCCGAGGTATCGCGCCCCTGACCTGATTGGGGCGCGGCCCTGGCTGTGGCTCTTGTCTGCAGTACCTGTCTTAACGTCAAGGAGTTTTGCAATGG  >  AM260479/696393‑696654

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