Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,134,614 G→A 8.9% A431A (GCG→GCA)  comEC → ComE operon protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,134,6140GA 8.9% 135.7 / 4.0 45A431A (GCG→GCA) comECComE operon protein 3
Reads supporting (aligned to +/- strand):  ref base G (14/27);  new base A (2/2);  total (16/29)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

CGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAGGGGTAATGCT  >  ACB143/2134481‑2134748
                                                                                                                                     |                                                                                                                                      
cGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCt                                                                                                                                      >  2:47175/1‑136 (MQ=255)
      cGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGgc                                                                                                                                >  2:220836/1‑136 (MQ=255)
            tCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCgcagca                                                                                                                          <  1:241456/136‑1 (MQ=255)
            gcAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCgcagca                                                                                                                          <  5:224317/135‑1 (MQ=255)
                             ccGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTcc                                                                                                         >  7:266754/1‑136 (MQ=255)
                             ccGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTcc                                                                                                         >  7:199874/1‑136 (MQ=255)
                             ccGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTcc                                                                                                         <  3:74014/136‑1 (MQ=255)
                                     gTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGcc                                                                                                 <  7:242125/136‑1 (MQ=255)
                                           gggCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCg                                                                                           >  1:138036/1‑136 (MQ=255)
                                                cGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTtgtg                                                                                      <  1:22378/136‑1 (MQ=255)
                                                    gATGTTGCCGCTGGCGGGTGTGGGCGAGGCGTTGCTCTGGCTGGCAGGGGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGcc                                                                                  <  2:251994/136‑1 (MQ=255)
                                                        tggCCGCTGGCGGGTGGTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCGGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgc                                                                              <  2:202229/134‑1 (MQ=255)
                                                           ccGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgctg                                                                            <  1:6673/135‑1 (MQ=255)
                                                           ccGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgctg                                                                            <  8:199874/135‑1 (MQ=255)
                                                            cGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGTCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgctgct                                                                          <  3:209926/136‑1 (MQ=255)
                                                            cGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgctgct                                                                          <  2:138036/136‑1 (MQ=255)
                                                            cGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgctgct                                                                          <  5:261117/136‑1 (MQ=255)
                                                            cGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgctgct                                                                          <  3:22345/136‑1 (MQ=255)
                                                            cGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgctgct                                                                          <  1:212139/136‑1 (MQ=255)
                                                            cGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgctgct                                                                          <  1:35945/136‑1 (MQ=255)
                                                            cGCTGGCGGGTGGTGGCGAGGCGGTGCTCTGGCTGGCAGGTGGCCTGCTGGATGGGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGgctgct                                                                          <  8:242885/136‑1 (MQ=255)
                                                              cTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTgg                                                                        >  1:274673/1‑136 (MQ=255)
                                                                                gCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGt                                                      <  2:161208/136‑1 (MQ=255)
                                                                                        ctggctggCAGGTGGCCTGCTGGATGTGCTGTTCCGTGTGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGACCTGCCGTTGTGGGCCTGGCTGCTGGTAATTCTGGGTGTGCGGTTGGTTTtg                                              >  4:152965/1‑136 (MQ=255)
                                                                                        ctggctggCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCAGCCGGCCCTGCCGGTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTtg                                              >  4:45762/1‑136 (MQ=255)
                                                                                         tggctggcAGGGGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGtttt                                               <  3:229760/134‑1 (MQ=255)
                                                                                           gctggcAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTtgctg                                           <  6:47284/136‑1 (MQ=255)
                                                                                           gctggcAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTtgctg                                           <  4:218070/136‑1 (MQ=255)
                                                                                           gctggcAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTtgctg                                           <  4:108489/136‑1 (MQ=255)
                                                                                           gctggcAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTtgct                                            <  3:60072/135‑1 (MQ=255)
                                                                                           gctggcAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTtg                                              <  1:21419/133‑1 (MQ=255)
                                                                                                  ggTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCgtg                                    <  8:283636/136‑1 (MQ=255)
                                                                                                   gTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCagcagc                                                                                                                        <  1:106053/51‑1 (MQ=255)
                                                                                                   gTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCagcagc                                                                                                                        >  2:106053/1‑51 (MQ=255)
                                                                                                   gTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCgtgg                                   >  4:122698/1‑136 (MQ=255)
                                                                                                            tgGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCt                          <  2:69903/136‑1 (MQ=255)
                                                                                                                     tggtCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGtt                 <  7:16315/133‑1 (MQ=255)
                                                                                                                     tGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGtt                 >  4:178067/1‑136 (MQ=255)
                                                                                                                         ccGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTggg               >  2:55196/1‑134 (MQ=255)
                                                                                                                            tGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAggg          >  7:150206/1‑136 (MQ=255)
                                                                                                                               gCTGGCACTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCt                          <  7:274621/117‑1 (MQ=255)
                                                                                                                               gCTGGCACTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCt                          <  7:278374/117‑1 (MQ=255)
                                                                                                                               gCTGGCACTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCt                          >  8:274621/1‑117 (MQ=255)
                                                                                                                               gCTGGCACTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCt                          >  8:278374/1‑117 (MQ=255)
                                                                                                                                 tGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCGGGGGGGGGTGGGAGGGGGAAt     >  6:173742/1‑136 (MQ=255)
                                                                                                                                    cgcTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAGGGGTAAtgct  >  8:131968/1‑136 (MQ=255)
                                                                                                                                    cgcTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAGGGGTAAtgct  >  2:173225/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                      
CGGTGCCGTGGGTCAGCCTGGCGGTCTTGCCGTTGGCGTTGCTGGGCACGCTGATGTTGCCGCTGGCGGGTGTTGGCGAGGCGTTGCTCTGGCTGGCAGGTGGCCTGCTGGATGTGCTGTTCCGTGGGCTGGCGCTGGTGGCGCAGCAGCGCCCGGCATGGGTCCCGCCGGCCCTGCCGTTGTGGGCCTGGCTGCTGGTATGTCTGGGTGTGCTGTTGGTTTTGCTGCCCCGTGGTGTGCCGCTGCGCGGGTTGGGAGGGGTAATGCT  >  ACB143/2134481‑2134748

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