Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,210,8810TG53.6% 22.7 / 36.2 31V204G (GTC→GGC) oprJ_2Outer membrane protein OprJ
Reads supporting (aligned to +/- strand):  ref base T (11/2);  new base G (0/15);  total (11/20)
Fisher's exact test for biased strand distribution p-value = 3.63e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-06
Rejected as consensus: Frequency below/above cutoff threshold.

AGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTGGTGC  >  ACB122/3210746‑3211012
                                                                                                                                       |                                                                                                                                   
aGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGt                                                                                                                                     >  8:52745/1‑136 (MQ=255)
                       aTCGCCCTGATCGCCGAGGTGAGCCCGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCCGGCGCTGCCCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCCGCGGCGCGCCCCCCGc                                                                                                               >  1:296554/1‑135 (MQ=255)
                                   gCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGGCCGCCCGGCCCCGgcgc                                                                                                  >  7:252616/1‑136 (MQ=255)
                                        ggTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGCTTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCCCCGCCCCCGCGCggggc                                                                                             >  7:164286/1‑134 (MQ=255)
                                        ggTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCCCCCGCCCCGCCCCCGCGCgggc                                                                                              >  1:27020/1‑134 (MQ=255)
                                               cAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCCGCCCGGCCCCCGCGCGGGCCCCCCAag                                                                                      >  6:267934/1‑136 (MQ=255)
                                                         gggCCTTGGCCGGCGCCGCGCGGCGCCAGGCGCTGACCCGGCCCCCCCTGCCCAGCGGGGTGGATCCCCTGGCGTGGGGCGGCCTGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTgg                                                                             <  2:95293/134‑1 (MQ=255)
                                                           gCCTTGGCCGGCGCGGCGCGGCGCCGGGCGCGCCCCCGCCACCCCCGGCCCCGCCGCTTGGAGTCCGGGGCGTGGGGCGGCCTGCGGCGCGCCGCCGGCAGGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGc                                                                          <  5:99217/136‑1 (MQ=255)
                                                           gCCTTGGACGGGGCCGAGCGGCGCAAGGCGCTCACCCTTCAGACCCTGGCCTGCCGCGAGGATTCCCTGGGGTTGGGCGGCCAGCGGGGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGc                                                                          <  8:25724/136‑1 (MQ=255)
                                                           gCCTTGGACGGCCCGGCGCGGCGCCCGGGGCGCCCCCGGCCCCCCGTGGCCCGCCGCGAGGATGCCCTGGGGTTGGGCGGCCCGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGc                                                                          <  6:177327/136‑1 (MQ=255)
                                                           gCCTGGGACGGCGCGGCGCGGCGCAAGGCGCTGCCCCGCCCGCCCCTGCCCAGCCGCGAGTAGACCCGGGCGTGGGGCGGCCAGCGGCGCGCCGCGGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGc                                                                          <  8:332789/136‑1 (MQ=255)
                                                            ccTGGGCCGGGGCCGGGCGGCGCCTGGCGCTCCCCGGTCCGCCCGTGGCCTGCCGGGTGTATGCCGGGGGGTGGGGGGGCCCGCGGCGCGCGGCGGGCAGGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGc                                                                          <  2:32537/135‑1 (MQ=255)
                                                               tGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCCCGGCCACGGCGCGGGGCCCCCAAGCGGCGCGGGGGCggg                                                                       >  2:4513/1‑135 (MQ=255)
                                                               tGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCCCCGGCCCCGCCACCGCGCGGGCCCCCCCAGCGGCGCGggggggggg                                                                      >  3:197502/1‑131 (MQ=255)
                                                                ggccGGCGCCGCGGGGCGCCTGGGGGGGACCCGTCGGCCCCTGGCCTGCCGCGAGGATGCCCGGGCGTGGGGCGGCGAGGGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTc                                                                     <  7:210949/133‑1 (MQ=255)
                                                                ggccGGCGCCGCGCGGCGCCCGGCGCGCACCCGGCCGCCCCGGCCCCGCCGCTAGGATGCCCGGGGGTGGGGGGCCCAGGGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTc                                                                     <  5:207334/133‑1 (MQ=255)
                                                                         cGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCCCCGGCCCGGCCCCCGGGCCGGGCCCCCAAGAGGCGGGGGGGGTGGTCGGGGAGtgc                                                            >  8:53769/1‑134 (MQ=255)
                                                                         cGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCCCCCGCACCGCCCCGGCGCGGGCCCCCCAAGAGGCGCTGGGGCGGGGCGCGCacgg                                                             >  2:35589/1‑132 (MQ=255)
                                                                             cggcgCCCGGCGCGCCCCCGCCCCCCCCGGGCCCGCCGCGCGGACCCCCGGGCGTGGGCCGCCCCGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGGCGAGCAGTCCCGt                                                         <  2:31251/135‑1 (MQ=255)
                                                                                  ccAGGCGCGTCCCCGCCATACCGGGGCCCGCCGTGTGGTTGCTGTGGCGTGGGCCGGCCCGGGGCGCGCCGCGGGCAGGGCGACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGa                                                    <  2:187616/135‑1 (MQ=255)
                                                                                            aCCCGTCCGCCCCTGGCCTGCCGCGTGGATTCCCGGGCGTTGGGCGGCCCGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGt                                           <  6:269731/134‑1 (MQ=255)
                                                                                                      ccgtGGCCAGCCGCGAGCATGCCCTGGCGTTGGTCGGCCTGCGGCGCGCGGCCGGCAGGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGcc                                <  1:279384/132‑1 (MQ=255)
                                                                                                           gCCAGCCGCGAGGATTCCCGGGCGGGGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAgc                          <  3:28448/136‑1 (MQ=255)
                                                                                                            ccAGCCGGGTGCATGCCCGGCCGTGGGCCGGCCCGCGGCGCCCCGCCGCCACGCCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCGGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAgc                          <  1:275430/135‑1 (MQ=255)
                                                                                                            ccAGCCGCGAGGATTCCCGGGCGTGGGGCGCCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAgc                          <  2:177104/135‑1 (MQ=255)
                                                                                                                   cgAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCCCGGCGCGGGCCCACCCAGAGGCGCGGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCCGGcc                  >  8:98740/1‑136 (MQ=255)
                                                                                                                            ccTGGCGTGGGTCGGCCCGGGGCGCGCCGCGGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAAcg           <  1:35589/134‑1 (MQ=255)
                                                                                                                                gCGTGGGGCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTgg     <  7:53769/136‑1 (MQ=255)
                                                                                                                                gCGTGGGGCGGCCAGCGGCGCGCCGCCGGCACGCCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTgg     <  3:334795/136‑1 (MQ=255)
                                                                                                                                gCGTGGGGCGCCCCGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTgg     <  5:267934/136‑1 (MQ=255)
                                                                                                                                  ggtGGGCGGCCAGCGGCGCGCCGCCGGCCCGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTGGTg   <  8:252616/134‑1 (MQ=255)
                                                                                                                                    tGGTCGGCCAGCGGCGCGCCGCCGGCCCGGCCCCGGCGCTGGCCCCCCCAGACGCGCGGGGGCTGGGGGGGCCGTCCCGGGGCGGGCGCGCAGTCCCCGTGCGCCCAACGCGCCAGGCCGACCACGCGCgggggg  >  2:64862/1‑132 (MQ=255)
                                                                                                                                       |                                                                                                                                   
AGGAAGTCGGGCGTGCCGCACGCATCGCCCTGATCGCCGAGGTGAGCCAGGCCTACCTGGCCTTGGACGGCGCCGAGCGGCGCCAGGCGCTGACCCGTCAGACCCTGGCCAGCCGCGAGGATTCCCTGGCGTTGGTCGGCCAGCGGCGCGCCGCCGGCACGGCCACGGCGCTGGACCACCAAGAGGCGCTGGGGCTGGTCGAGCAGTCCCGTGCCGAGCTCGAAGTCACCGTGCGCCAACAGCGCCAGGCCTACAACGCGCTGGTGC  >  ACB122/3210746‑3211012

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: ‑