Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,584,539 A→C 23.4% G137G (GGT→GGG)  rhaS_1 ← HTH‑type transcriptional activator RhaS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,584,5390AC23.4% 73.6 / 3.6 49G137G (GGT→GGG) rhaS_1HTH‑type transcriptional activator RhaS
Reads supporting (aligned to +/- strand):  ref base A (29/7);  new base C (2/11);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 5.86e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCAGGCGGTAGAGAATCTCGCGGCGGAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGC  >  ACB122/1584409‑1584672
                                                                                                                                  |                                                                                                                                     
gcagcaGGCGGTAGAGAATCTCGCGGCGGAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCAc                                                                                                                                   >  6:453110/1‑135 (MQ=255)
 cagcagGCGGTAGAGAATCTCGCGGCGGAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCCCATACAACCCACg                                                                                                                                  >  8:52581/1‑135 (MQ=255)
   gccgggggTAGAGAAGCGCGCGGCGCAACCGCGGCGCCAGCTTTGCGATGGCGCGTGGGGTGCCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGAACCGCGGCAGCTGTCCTTCCCACAAACAACCCACGGc                                                                                                                                <  6:165233/129‑1 (MQ=255)
        cGGTAGAGAATCTCGCGGCGGAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGaa                                                                                                                           >  4:277431/1‑135 (MQ=255)
            agagAATCTCGCGGCGGAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGt                                                                                                                        >  8:316982/1‑134 (MQ=255)
               gaATCTCGCGGCGGAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCAGGAAGGTCGAt                                                                                                                    >  1:25058/1‑135 (MQ=255)
                   ctcGGGGCGGACCGGGGGGGCCTGCATGGCGTGGGCGGGTGGGGGGTCCAGCCGCCGCAGCACCCCCCGCCGGGCACCGCGTACCGCGGAATCTGTCCTTACCACAAACAACCCACGGCCGGAAGGTCGATTggg                                                                                                                <  2:336470/135‑1 (MQ=255)
                   cccgcgCCGAACCACCGGCGCCCGCATTGGGAGGTCGCGTGGGGTGCCCACCCCCCGCCGCAACCGCAGCAGGGCATCGATAACCGCGGCAGCTGCCCTTCCCACAAACAACCCACGGCCGGAAGGTCGATTggg                                                                                                                <  8:202318/133‑1 (MQ=255)
                        ggcggAACAGCGGCGCCTGCATTGGGATGTCGCGAGGGGTGCCCAGCAGCCGCAGAAACCGCAGCAGGGCATCGATAAAGGCGGCATCTGTCCTTTCCAAAACCAAACCACGGCCGGAAGGTCGATTGGGTACCa                                                                                                           <  7:342570/135‑1 (MQ=255)
                        ggcggAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCa                                                                                                           >  7:55975/1‑135 (MQ=255)
                        ggcggAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCa                                                                                                           >  7:322191/1‑135 (MQ=255)
                        ggcggAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCa                                                                                                           <  7:316982/135‑1 (MQ=255)
                        ggcggAACAGCGGCGCCAGCATTGCGATGGCGCGAGGGGTGGCCAGCACCCGCAGCAACCGCAGCAGGGCAGCGAGTAACGCGGCATCTGTCCTTCCCACAAACAAACCACGGCCGGAAGGTCGATTGGGTACCa                                                                                                           <  7:57979/135‑1 (MQ=255)
                                     cgcCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGct                                                                                              >  5:249856/1‑135 (MQ=255)
                                                  aTGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGt                                                                                 >  2:288467/1‑135 (MQ=255)
                                                        cgAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTg                                                                           >  7:392624/1‑135 (MQ=255)
                                                        cgAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCAATGGCCAGCTCTCGGCGATCAGTTGGCTg                                                                           >  4:4332/1‑135 (MQ=255)
                                                              gtggCCAGCCGCCGCAGCAACCCCCGCAGGGAATCGAGTCACGGCGCTTCTGCCCTAACCACAAAAAACCCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTcc                                                                     <  5:350992/135‑1 (MQ=255)
                                                                  ccAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGCTCAGTTGGCTGATCTccgccg                                                                 >  4:282314/1‑135 (MQ=255)
                                                                  ccAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTccgccg                                                                 >  5:231707/1‑135 (MQ=255)
                                                                   cagcagCCGCAGCACCCGCAGCAGGGCATCGAGAAACGCGGCATCTGTCCTTACCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCgg                                                                <  8:289818/135‑1 (MQ=255)
                                                                    agcagcCGCAGCAACCGCAGCAGGGCATCGATAAACGCGGCATCTGGCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCggg                                                               <  6:313855/135‑1 (MQ=255)
                                                                    agcagcCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCggg                                                               >  5:129303/1‑135 (MQ=255)
                                                                    agcagcCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCggg                                                               >  5:424966/1‑135 (MQ=255)
                                                                    agcagcCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCAATGGCCCGCTCTCGTCGATCAGTTGGCTGATCTAAGCCggg                                                               >  2:320712/1‑135 (MQ=255)
                                                                            cagcaACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAACCCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGggtccaggt                                                        >  4:262529/1‑134 (MQ=255)
                                                                              gcaACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGggtccaggtcca                                                     >  6:204737/1‑135 (MQ=255)
                                                                               cacccgccgCGGGGCATCGAGCACCGCGGCATGTGTCCTTTCCACAAACAACCCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCAt                                                    <  3:340262/132‑1 (MQ=255)
                                                                               caACCGCAGCAGGGCAGCGAGAAACGCGGCAGGTGTCCTTTCCACAAACACACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCAt                                                    <  3:152135/135‑1 (MQ=255)
                                                                                aaCCGCAGCCGGGCAGCGATTACCGCGTCATCTGTCATTTCCACATACAACCCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATg                                                   <  5:98187/135‑1 (MQ=255)
                                                                                     cagcagGGCTGCGAGTAACGCGGCATCTGTCCTGTCCACATACACCCCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAg                                              <  8:77221/135‑1 (MQ=255)
                                                                                     cagcagGGCATCGAGTACCGCGGCATCTGTCCTTACCACAAACAACCCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAg                                              <  1:117215/135‑1 (MQ=255)
                                                                                          gggCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCa                                         >  5:309116/1‑135 (MQ=255)
                                                                                                   gTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGCTCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCCGGCTTGGGt                                >  1:223683/1‑135 (MQ=255)
                                                                                                     aaCGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTtc                              >  1:305825/1‑135 (MQ=255)
                                                                                                     aaCGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTACATGCGCAGGCCCAGGCTTGGGTtc                              >  3:90706/1‑135 (MQ=255)
                                                                                                         cgGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCgg                          >  2:453202/1‑135 (MQ=255)
                                                                                                               cTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCACCTCCATGCGCAGGCCCAGGCATGCGGACTCGGGGCtcg                    >  4:254142/1‑133 (MQ=255)
                                                                                                                tggccTTGCCACATACAACCCTCGGCCGGAAGGGCGATGGGGGACCAGCATTGGCCGGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCTGGTCCATGCGCAGGCCCAGGCTTGGGTGCTCGGGGCTGGc                   <  8:440334/132‑1 (MQ=255)
                                                                                                                    cTTCCCACATACACCCCACGGCCGGAAGGGCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATcc               <  3:223705/135‑1 (MQ=255)
                                                                                                                         cacaTACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGcc           >  3:414614/1‑134 (MQ=255)
                                                                                                                         cacaTACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGcc           >  7:371417/1‑134 (MQ=255)
                                                                                                                         cacaTACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCt          >  7:97712/1‑135 (MQ=255)
                                                                                                                         cacaAACAACCCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGcc           <  3:167475/134‑1 (MQ=255)
                                                                                                                          acaTACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTc         >  6:202881/1‑135 (MQ=255)
                                                                                                                          acaTACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTc         >  3:26044/1‑135 (MQ=255)
                                                                                                                          acaTACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGATGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCGTGCATCCAGccac         >  2:238835/1‑133 (MQ=255)
                                                                                                                               cAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGgc    <  1:288467/135‑1 (MQ=255)
                                                                                                                               cAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGgc    >  8:53793/1‑135 (MQ=255)
                                                                                                                                 aaCCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGgcgc  <  1:320712/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
GCAGCAGGCGGTAGAGAATCTCGCGGCGGAACAGCGGCGCCAGCATTGCGATGTCGCGAGGGGTGTCCAGCAGCCGCAGCAACCGCAGCAGGGCATCGAGTAACGCGGCATCTGTCCTTTCCACATACAAACCACGGCCGGAAGGTCGATTGGGTACCAGCATTGGCCCGCTCTCGGCGATCAGTTGGCTGATCTCCGCCGGGTCCAGGTCCATGCGCAGGCCCAGGCTTGGGTTCTCGGGGCTGGCATCCAGCCTCACGGCGC  >  ACB122/1584409‑1584672

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