Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,679,3570TG43.2% 47.0 / 14.8 45E474D (GAA→GAC) rcsC_4Sensor histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base T (1/24);  new base G (19/0);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.42e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCACCCAGCGTGCGGCGCAGCAGCTCTTCCAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTCAGCAGGTTATTGAAGTCATGCGCCACGCCAGCCGTTAG  >  ACB122/2679223‑2679488
                                                                                                                                      |                                                                                                                                    
ccACCCAGCGTGCGGCGCAGCAGCTCTTCCAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCtt                                                                                                                                     <  8:40545/136‑1 (MQ=255)
 cACCCAGCGTGCGGCGCAGCAGCTCTTCCAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCt                                                                                                                                      <  3:235943/134‑1 (MQ=255)
           gcggcgCAGCAGCTCTTCCAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATg                                                                                                                          <  6:424791/135‑1 (MQ=255)
                       ctctTCCAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTCGC‑GGATATGGGGGAGGCGCGt                                                                                                              >  4:347637/1‑135 (MQ=255)
                            ccAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCGCCCGGCGG‑GGGGTTTGGTGAGGGGCGTGCggg                                                                                                         >  1:243142/1‑135 (MQ=255)
                            ccAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCgg                                                                                                          <  7:409823/134‑1 (MQ=255)
                            ccAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTCGC‑GGATATGGGTGAGGCGAGTGCGGGc                                                                                                        >  7:369432/1‑136 (MQ=255)
                             cAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGcag                                                                                                          <  4:144820/133‑3 (MQ=255)
                               ggccggccAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTCGC‑GGATATGGGGGAGGCGAGTGCCGGcgc                                                                                                      >  2:149057/1‑133 (MQ=255)
                                   ggccAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTCGCGGGATA‑GGGGGAGGCGCGAGCGGGCCGTGTc                                                                                                  >  5:64802/1‑135 (MQ=255)
                                    gccAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCa                                                                                                <  8:300085/136‑1 (MQ=255)
                                    gccAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGGGAGGCGAGTGCGGGCCTTGTCCa                                                                                                <  8:255838/136‑1 (MQ=255)
                                     ccAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCa                                                                                                <  1:260439/135‑1 (MQ=255)
                                     caggggTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCa                                                                                                <  5:241624/132‑1 (MQ=255)
                                       aGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTCGC‑GGCTATGGGGGAGGCGCGTGCGGGCCGTGTGCAgg                                                                                              >  5:183471/1‑135 (MQ=255)
                                         gggaCGGTTGAGGGTCAGCGGCACCGGCGCCAGACGCGGGCGGCGGGAGAAGGCCAGCAACTGCGAGGGCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTc                                                                                            <  6:258522/132‑1 (MQ=255)
                                                 tGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTGGC‑GGCTATGGGTGGGGCGAGTGGGGGCCTTGTACAGGTCGGCCTTGGc                                                                                   >  8:55675/1‑136 (MQ=255)
                                                 tGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCCCGTAGC‑GGCTATTGGTTAGGCGCGTGCGGGCCGTGTACACGTAGGCCGTGGc                                                                                   >  8:97130/1‑136 (MQ=255)
                                                 agAGGTTCAGCGGCACCGGCTCCACACGCTCGCGGCGGGAGAAGGCCAGCAACTGCGAGCTCAGCTTGCCACCACGCTCCCCCGCCTCGC‑GGATATCGCTGTGGCGATCGCGGGCCGTGTCCACGTCGGCCATgg                                                                                    >  6:90496/2‑135 (MQ=255)
                                                  gAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCGCGCGGGGC‑GGCAATTGGTGGGGCCAGTGGGGGGCTTGTGCAGGTCGGGCCTggg                                                                                   >  3:23440/1‑134 (MQ=255)
                                                      ttCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTGGC‑GGGTGTGGGGGGGGCGCGGGCGGGGCGTGGGCAGGGCGGCCGTgggccgcg                                                                              >  7:47332/1‑130 (MQ=255)
                                                      ttCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTGGC‑GGGTATGGGTGAGGTGAGTGCGGGCCGTGGCCCGGTCGGCCGTGGGCAgct                                                                              >  7:58044/1‑134 (MQ=255)
                                                      ttCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTCGC‑GGCTATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCGTGGCCCGa                                                                               >  4:307716/1‑135 (MQ=255)
                                                      ttCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCCCGTCGC‑GGATATGGGTGTGGCGCTTGCGGGCCTTGTCCAGGTCGGGCTTgggccgc                                                                               >  4:103867/1‑130 (MQ=255)
                                                        cAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGAGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATc                                                                             <  2:16528/135‑1 (MQ=255)
                                                        cAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGTCGC‑GGTTATGGGTGAGGCGCGTGCGGGGCTTGTCCAGGTCGGCCGTGGGCAGATcg                                                                            >  8:334002/1‑136 (MQ=255)
                                                           cGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCACCTGCGAGGTCAGCTTGGCACCACGCTCGCGCGCGTGGG‑GGTTATGGGTGGGGCGCGTGCGGGGCGTGTGCCGGTCGGGCGTGGGCGGCTcgcg                                                                          >  3:211529/1‑135 (MQ=255)
                                                             gCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCGGAGG‑GGATTAGGGTGCGGCCCGTGGGGGCCGAGTCCAGGTCGGCCTGGGCCAGATGgcgca                                                                        >  2:103132/1‑134 (MQ=255)
                                                                         gACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTgg                                                            <  1:9838/135‑1 (MQ=255)
                                                                           cGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTGCGC‑GGTTATTGGTGGGGCGAGTGGGGGCCGTGTGCACGTGGGCCGTGGGCAGATCGCGCTCCAGCAACCGGGcc                                                          >  6:412130/1‑134 (MQ=255)
                                                                            gCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGGCAGCTTGGCACCACGCTCGCCCGCTTCGC‑AGATATGGGGGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCAcc                                                        <  7:295282/136‑1 (MQ=255)
                                                                             cTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCAcc                                                        <  4:420907/135‑1 (MQ=255)
                                                                                 cggcggGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCgg                                                      <  4:202481/133‑1 (MQ=255)
                                                                                 cggcggGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCCCGTCGC‑GGAGATGGGGGCGGCGAGTGCGGGCCGTGGCCAGGGCGGGCGTGGGCCGAGCGCGCTCCCGCCAGCGGGCAcccgggg                                                   >  8:310814/1‑131 (MQ=255)
                                                                                   gcggGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGt                                                 <  8:396503/136‑1 (MQ=255)
                                                                                           aGGCCAGCAACTGCGAGGTCAGCTTGGCAGCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGt                                         <  8:149950/136‑1 (MQ=255)
                                                                                           aGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCg                                          <  4:145808/135‑1 (MQ=255)
                                                                                             gCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCGCGCGTCGC‑GGAGTTGGGGGTGGCGCGTGCGGGCCGTTTCCAGGGCGGCCGTGGGCCGATCGCGCGCCAGAAAGCCGGCCCCCGGGCGTCTCACCCTc                                        >  4:286641/1‑135 (MQ=255)
                                                                                             gCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTca                                       <  7:82281/136‑1 (MQ=255)
                                                                                                      tGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTCAGCAGGTTa                               <  3:309788/135‑1 (MQ=255)
                                                                                                      tGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTCAGCAGGTTa                               <  5:406090/135‑1 (MQ=255)
                                                                                                      tGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTCAGCAGGTTa                               <  1:45435/135‑1 (MQ=255)
                                                                                                      tGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCGTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTCAGCAGGTTa                               <  6:89119/135‑1 (MQ=255)
                                                                                                            gTCAGCTTGGCACCACGCTCGCGCGCGGGGC‑GGGTATTGGTGTGGCGCGTGGGGGCCGTGGGCACGGCGGGCGTGGGCCGCTCGCGCGCCCGCAAGCGGGGCCCGGGGAGGTTCACCCGCAGCCGGGTTTTGatt                         >  6:420990/1‑133 (MQ=255)
                                                                                                                         cacGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTCAGCAGGTTATTGAAGTCATGcgccacgcc           <  7:20926/136‑1 (MQ=255)
                                                                                                                              tCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTCAGCAGGTTATTGAAGTCATGCGCCACGCCAGcc       <  2:359346/135‑1 (MQ=255)
                                                                                                                                   cGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTCAGCAGGTTATTGAAGTCATGCGCCACGCCAGCCGTTAg  <  2:269404/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                    
CCACCCAGCGTGCGGCGCAGCAGCTCTTCCAGGCCGGCCAGGGTACGGTTGAGGTTCAGCGGCACCGGCTCCAGACGCTGGCGGCGGGAGAAGGCCAGCAACTGCGAGGTCAGCTTGGCACCACGCTCGCCCGCTTCGC‑GGATATGGGTGAGGCGAGTGCGGGCCTTGTCCAGGTCGGCCTTGGCCAGATCGCGCTCCAGAAAGCTGGCACCGGTGAGTATCACCGTCAGCAGGTTATTGAAGTCATGCGCCACGCCAGCCGTTAG  >  ACB122/2679223‑2679488

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