Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24402,391,0260TG58.8% 3.4 / 11.9 19Q561P (CAG→CCG) bvgSVirulence sensor protein BvgS
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (10/0);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 5.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTT  >  KT2440/2390905‑2391132
                                                                                                                         |                                                                                                          
gTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCGGGGGGATATgccgg                                                                                              >  8:179133/1‑136 (MQ=255)
gTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTGGGCGCGCGGGCGCGGGGGCATCTgacgg                                                                                              >  7:176217/1‑132 (MQ=255)
gTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCCCGCGGTCGGCGCGCTGGCGCAGGCGAAACAgcagg                                                                                              >  1:16347/1‑136 (MQ=255)
       cAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTgcaggcgcaggc                                                                                       <  7:71817/136‑1 (MQ=255)
       cAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGGGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTgcaggcgcaggc                                                                                       <  1:46863/136‑1 (MQ=255)
            tttCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCCGTCGGCGCGCAGGCGCGGGCGAATCTGCAGGCGCCGGCGGGcc                                                                                  >  4:28155/1‑135 (MQ=255)
                                     agatacagGGGGGGGGGCGTCCCTTTGAGCCGGGCGCGCATGATCGCCATGGGGCCTGACCGGGCGGGTGCGGCGCGGTGGCGCTGGCGAATCTGCAGGCGCTGGCGGGCATTCCAGCGCAGGGCCAGCAACaaac                                                         <  4:158794/134‑2 (MQ=255)
                                          catGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTAGGCGCGCGGGCGCAGGCGCATATGGAGGCGCACGCCGGCCATACACAGCAGGGCCAGCCACAAAACCAcc                                                    >  8:203554/1‑135 (MQ=255)
                                                gtggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCTGGCGAATCTGCAGGCGCCGGCGGGCCTTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTg                                              >  4:211296/1‑136 (MQ=255)
                                                gtggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTGGAGCGCGCGGGCGGCGCGCGGGCGCGGGCGCCAATGCAGGCGCCCGCCGGGAGTCCAGACCACGGCCACCACCAAAAACACCGcccgg                                              >  2:105219/1‑133 (MQ=255)
                                                gtggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTGCAGAGCGCCTTCGGCGCGCACGCGCGGGGGAAACTGCCGGCGCCCGCGGGCATTACAGAGCAGGGCCAGGACCAAAATCACCGCCATg                                              >  6:47655/1‑136 (MQ=255)
                                                  ggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTGAGCGCGCGTTCGGCGCGCGGGCGCGGGCGCATCTGCCGGCGCCCGCGGGGCTTACAGAGCAGGGCCCGCCACAAAAGCACCGCCGGGGc                                            >  6:236081/1‑136 (MQ=255)
                                                        tACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCGGGCGCATCTGCAGGCGCAGGCGGGCAGTCCCGAGCAGGGCCAGCAACCAAAGCACCGCCGTGGCCCGccc                                      >  2:152223/1‑134 (MQ=255)
                                                                aGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCGGGCGGATCTGCAGGCGCAGGCGGGGATTACAGCGCAGGGGCCGCACCAAAAGCACGGCCGTGGCCCGCCCCCCCAGcc                              >  2:61799/1‑135 (MQ=255)
                                                                        gcgcATGAATGCCAACTGGTCGTGCAGCGCGCGGGGGGCGCGCGGGCGCGGGGGACTCGGCAGGCGCCCGCGGGCCTTGCAGCGCACGGGCCGCCACAAAAGCCCCGCCGTGGGCCGCGCCCCCCGCCGGGCCCgg                      >  2:27317/1‑136 (MQ=255)
                                                                               aaggCCACCTGGGGGTTCAGCGGGGGTGCGGCGCGCGGGCGCTGGCGAAGCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGt               <  7:30334/133‑1 (MQ=255)
                                                                                           gcGTTCAGCGGGGGTGCGGCGCGCGGGCGCTGGCGAAGCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGt   <  7:179133/135‑1 (MQ=255)
                                                                                            ggTTCCGCGCGGGGGCGGCGCGGTGGCGCTGGCGAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGtt  <  2:112007/135‑1 (MQ=255)
                                                                                            cGTTCAGCCCGCGTTCGGCGCGCTGGCGCGGGCGAATCTGCAGGCGCCGGCGGGCATTCCCGAGCCGGGCCCGCAACAAAACCACGGCCGTGGCCAGCCCCCCCAGCCGGGCCCGGCGGCGGTTGCCGTTGGAGtg  >  8:141772/1‑135 (MQ=255)
                                                                                                                         |                                                                                                          
GTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTT  >  KT2440/2390905‑2391132

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