Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,390,993 T→G 30.3% E566D (GAA→GAC)  bvgS ← Virulence sensor protein BvgS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,390,9930TG30.3% 30.8 / 2.3 43E566D (GAA→GAC) bvgSVirulence sensor protein BvgS
Reads supporting (aligned to +/- strand):  ref base T (20/9);  new base G (12/2);  total (32/11)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTG  >  ACB143/2390865‑2391116
                                                                                                                                |                                                                                                                           
gCCTGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTGAGCGCGCGGTCGgcgc                                                                                                                      >  4:149131/1‑136 (MQ=255)
   tGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCgcgct                                                                                                                   >  2:285651/1‑136 (MQ=255)
        ggCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCgcgccggcgc                                                                                                              >  4:194193/1‑136 (MQ=255)
        ggCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTGCAGCGCGCGTTAGGCgcgctggcgc                                                                                                              >  5:256678/1‑136 (MQ=255)
             cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCTCATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCCGCGGGGGGTCCGGGcgcgggcgcggggc                                                                                                         >  8:287413/1‑134 (MQ=255)
             cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCGCATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTAAGCGCGCGTTCGGCGCGcaggcgctggcg                                                                                                         >  2:3669/1‑136 (MQ=255)
             cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGcgctggcgctgg                                                                                                           <  7:278522/134‑1 (MQ=255)
             cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCgctggcgcgggc                                                                                                          >  3:311161/1‑135 (MQ=255)
             cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCCCGGTCGGCGccctggcgcggggc                                                                                                         >  3:71496/1‑134 (MQ=255)
             cAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCCCCTTCGGCCCGcaggcgcgggcg                                                                                                         >  6:112357/1‑136 (MQ=255)
                 gTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATc                                                                                                     <  4:97761/136‑1 (MQ=255)
                       gTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTGAGCGCGCGGTCGGCGCGCGGGCGCGGGGGAATCTgcagg                                                                                               >  1:103791/1‑136 (MQ=255)
                       gTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCCGCGCGCGTTGGGCGCGCAGGCGCTGGCGAATCTgcagg                                                                                               >  7:128368/1‑136 (MQ=255)
                            tGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTgcaggcgca                                                                                           <  4:137187/135‑1 (MQ=255)
                              cAACTTTTCCGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTGCGGCGCGCTGGCGCTGGCGAATCTgcaggcgcaggc                                                                                        <  2:259477/136‑1 (MQ=255)
                              cAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTgcaggcgcaggc                                                                                        <  5:65074/136‑1 (MQ=255)
                                            aCCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTAAGCGCGCGTTCGGCGCGCTGGCGCGGGCGAATATGCAGGCGCAGGCGGGCATTACagat                                                                           >  8:56094/1‑134 (MQ=255)
                                            aCCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTAAGCGCCCGTTGGGCGCCCAGGCGCTGGGCAAAATGCACGCCCACGCCGGCCTTACAGAGc                                                                          >  6:83621/1‑136 (MQ=255)
                                                                       gtggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTGCCGCCCGCGGTCGGCGCGCGTGCGCGGGCGCATCTGCAGGCGCCGGCGGGCAGTCCCGAGCCGGGCCCGCCACACACACACGGCCCTg                                               >  6:245640/1‑136 (MQ=255)
                                                                         ggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGc                                             >  2:198948/1‑136 (MQ=255)
                                                                         ggggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTGAAGCGCGCGTTAGGCGCGCTGGCGCGGGCGAATCTGCAAGCCGAGGGGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCCCCGCCGTGGc                                             >  7:134229/1‑136 (MQ=255)
                                                                           ggCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATATGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCa                                           >  5:127947/1‑136 (MQ=255)
                                                                             cGTACGGTTGAGCTGGGCGCGCATGAATGCCACGGGGTCGTGCTGGGGGGGTTCGGCGGGGTGGCGCTGGCGAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGc                                         <  2:323556/136‑1 (MQ=255)
                                                                             cGTACCGTTGCGCAGGGCGCGCATGAATGCCAAGTGGGCGTTCAGCGCGCGTTCGGCGCGGTGGCGCTGGCTAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGAAACAAAAGCACGGCCGTGGCCAGc                                         <  5:111764/136‑1 (MQ=255)
                                                                             cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCAGGCGCTGGCGAATCAGCAGGCGCCGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAACCACGGCCCTGGCCCGc                                         >  8:72090/1‑136 (MQ=255)
                                                                             cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTAGGCGCGCTGGCGCTGGCGAATATGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGc                                         >  7:58630/1‑136 (MQ=255)
                                                                             cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTGGGCGCGCAGGCGCGGGGGCATCTGCAGGCGCAGGCGGGCCTTCCAGAGCAGGGCCCGCAACAAAAGCACCGCCGTGGcccc                                          >  5:88486/1‑133 (MQ=255)
                                                                             cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCGGGCGAATCTGCACGCCCACGCGGGCATTCCAGAGCACGGCCAGCAAAAAAACCACCGCCGAGGCCac                                          >  5:311926/1‑134 (MQ=255)
                                                                             cGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTAAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATATGCAGGCGCAGGCGGGCATTCCAAAGCAGGGCCAGCAACAAAAGCACCGCCGTGGCCac                                          >  3:129876/1‑134 (MQ=255)
                                                                               tACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCGGGCGCATATGCAGGCGCCGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGccc                                       >  4:135552/1‑134 (MQ=255)
                                                                               tACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTAAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTGCCGGCGCAGGCGGGCAATCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGGTGCCCGcac                                       >  8:205108/1‑136 (MQ=255)
                                                                               tACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTAACCGCGCGTTAGGCCCCCTGGCGCGGGCGAATCTGCAGGCGCAGGCGGGGATTCCAGAGCCGGGCCCGCAAAAAAAACACCGCCGTGGCCAGccc                                       >  2:43001/1‑134 (MQ=255)
                                                                                 cctttGCGCCGGGCGCGCATTATTGCCACGTGGTCGTGCTGCGCGCGTGGGGGGCGGGGGCGCTGGCGAATCTGCCGGCGCAGGCGGGCTTTCCAGCGCAGGCCCAGCAACAAAACCACGGCCGTGGCCAGcacca                                     <  4:71496/133‑1 (MQ=255)
                                                                                   gTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTAGGCGCGCTGGCGCTGGCGCATATTCACGCGCAGGCGGGCATTACAGAGCAGGGCCAGCAACAAAAGCACGGCCCTGGCCAGcgcccc                                    >  4:282358/1‑133 (MQ=255)
                                                                                   gTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGGTCGGCGCGCTGGCGCTGGCGAATATGCCGGCGCCGGCGGGGCTTACAGAGCAGGGGCCGCAACAAAAGCACGGCCCTGGCCCGccccac                                    >  5:15049/1‑135 (MQ=255)
                                                                                   gTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTAACCGCGCGTTGGGCGCCCTGGCGCGGGCCAATCAGCAGGCGCAGGCGGGCCTTACAGAGCACGGGCAGCAACAAAATCACGGCCGTGGCCCCCaccccc                                   >  1:297752/1‑136 (MQ=255)
                                                                                            ggCGCGCATGAATGCCAACTGGTCGTTCACCCCGCGGTAGGCGCGCGGGCGCGGGGGTAAAGTCACGGGCACGCCGGCCTTAACGACCACGGCCAGCAACAAAAACAACGCCGGGGGCCGCACCCCCCGCCGGGcg                          >  4:157795/1‑135 (MQ=255)
                                                                                               gcgcATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACCGCCGTGGCCCGCACCCCCTGCAGGGCCAgg                       >  2:188097/1‑136 (MQ=255)
                                                                                                   aTGAATGCCACCTGGGCGTGCCGCGCGCGGGCGGCGCGCTGGCGCTGGCGAAGCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAggcggc                   <  5:49897/136‑1 (MQ=255)
                                                                                                     gAATCCCAACTGGTCGTTCCGCGCGCGTGCGGCGCGCGGGCGCTGGCGAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCt                  <  3:149131/135‑1 (MQ=255)
                                                                                                         gccacgtGGTCTTTGCGCGGGCGGGCGGCGCTCTGGCGCTGGCTATTCTGCAGGCGCTGGCGGGCATGCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGc             <  4:233943/130‑1 (MQ=255)
                                                                                                         gCCAACTGGTCGTTCAGCCCGCGGTCGGCCCGCATGCGCTGGGGAATATTCACGCGCAGGCGCGCATTACAGATCAAGGACAGCCCCAAAACCACCGCCGTGGCCAGCACCCCCATCAAGGCCAGGCGGCTATAgg             >  6:122076/1‑135 (MQ=255)
                                                                                                                    gTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATATGCAGGCGCAGGCGGGGATTACAGAGCAGGGCCCGCAACAAAAGCACGGCCGTGGCCAGCCCCCCCAGCAGGGCCAGGCGGCTGTAGCCGTTGGGGTTg  >  8:167000/1‑136 (MQ=255)
                                                                                                                                |                                                                                                                           
GCCTGCACGGCTTCAAGGTAGCTGTAGTTGCAACTTTTCAGGCAACCCTCGCGGTCGCGCACATACATGGGGTGGGGCGTACCGTTGAGCAGGGCGCGCATGAATGCCAACTGGTCGTTCAGCGCGCGTTCGGCGCGCTGGCGCTGGCGAATCTGCAGGCGCAGGCGGGCATTCCAGAGCAGGGCCAGCAACAAAAGCACGGCCGTGGCCAGCACCACCTGCAGGGCCAGGCGGCTGTAGCCTTTGGAGTTG  >  ACB143/2390865‑2391116

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