Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,585,043 T→G 32.3% intergenic (‑94/‑109) rhaS_1 ← / → CJ019_01492 HTH‑type transcriptional activator RhaS/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,585,0430TG32.3% 25.0 / 3.1 32intergenic (‑94/‑109)rhaS_1/CJ019_01492HTH‑type transcriptional activator RhaS/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (5/14);  new base G (8/2);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCC  >  ACB122/1584912‑1585176
                                                                                                                                   |                                                                                                                                     
aGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTg                                                                                                                                      <  6:391892/133‑1 (MQ=255)
     cAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCGGCAGGa                                                                                                                                 >  1:271626/1‑133 (MQ=255)
          tcggttcggtGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCAGAAACAGCCCCCCAGCAGGCACAGGc                                                                                                                           >  5:394967/1‑134 (MQ=255)
               tagctGGACGTGGCGGTGGTCATGATGGCTATCCTATGCTGACCGACACATCATAGATTTGGGCAAAGACGCGCTCTAGTCGAAAGATGCATTTGGAATGCAATATACAGCCCCCCAGCAGGAACAGGCAATCGa                                                                                                                     >  4:155169/5‑135 (MQ=255)
                      aCGGGGCGGTGGTCATGATGGCTTCCCTCTGCTGACCGACGCAGAATAGGTTTGGGCAACGACGGGCGCTAGTCTAAAGCGGCATGTGGATTGCCTGATCCGGCCGCGCTGCAGGATCAGGCAATCGACATGCAg                                                                                                              <  6:394967/135‑1 (MQ=255)
                       cGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCAGATACTTCCCCCCCGCACGCTCAGGCAATACACATGCAGt                                                                                                             >  7:81021/1‑135 (MQ=255)
                             ggtggtCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCAGATACTGCCGCGCTGCAGGCTCAGGCAATCGACAGGCAGTAATgg                                                                                                        >  6:319458/1‑134 (MQ=255)
                                       atgGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATAGACAAGCAGTAAATGCCTAAcgcgcc                                                                                             >  4:307584/1‑135 (MQ=255)
                                                 tctGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCGGGTACTGCCCCGCGGCCGGGTCCGGCCATCGCCATGAAGTAATGGCCCAACCCGCCGCGCGGGTGAg                                                                                  >  8:217846/1‑136 (MQ=255)
                                                 tctGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCGGATCCTGCCGCGCGGCCCGCTAAGGGAATCGCCAGTCACTAAAGGGCGAACCCCCCGCGCgggggg                                                                                   >  3:374774/1‑132 (MQ=255)
                                                           gaggcagcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGCTCCGGCCGCGCGGGCGGCTCAGGCCATCGACCTGGAGTTATGGCCTAACGCGCCCCGCGGGTGGGCGCCGAAcc                                                                         >  3:198926/4‑135 (MQ=255)
                                                           gACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGCAACAGCCCCGCGGCACGTTGCGGCAATCGACAAGCACTAATGGCCGTACGCGCCGCGCGGGTGCGCGCCGAAcca                                                                        >  8:296410/1‑135 (MQ=255)
                                                                agcaTAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCAGATCAGGCCGCGCGGGAGGCTCAGGCAAACGACATGCAGTATTTGCCTAACGCGCCGCGCGGGGGGGCGCCGAACCGTGGc                                                                    >  6:317711/1‑135 (MQ=255)
                                                                 gcaTAGGTTTGGGCAAAGACGGGCGCTAGTCTAAAGCTGCATGTGGATTGCCTGATCCGGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCa                                                                   <  6:167158/135‑1 (MQ=255)
                                                                          tGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATAGCCTGATACTGCCGCGCGGCAGGTTCAGGCCATCGCCATGGAGTAATGGCCGAACGCGCCGCGCGGGTGCGCGCCGAACCTTGGCAGCCTGGcgt                                                          >  4:222906/1‑133 (MQ=255)
                                                                               aaaGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATACTGCCGCGCTGCAGGATAATGCAATCGACATTCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCAAACCTTGGCAGCCTGGCTTTtccgt                                                    >  8:287903/1‑136 (MQ=255)
                                                                                aaGACGGGCGCTAGTCGAAAGCTGCATGTGGATGGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtc                                                   <  7:121428/136‑1 (MQ=255)
                                                                                aaGACGGGCGCAAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtc                                                   <  7:316257/136‑1 (MQ=255)
                                                                                 agccGGGCGCTAGGCGAAAGCTGCATGTGGATGGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtc                                                   <  1:34253/132‑1 (MQ=255)
                                                                                 aGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTtccgtc                                                   <  2:271626/135‑1 (MQ=255)
                                                                                    cGGGCGCTAGTCGAAAGCTGCATGTGGATTTCCGGGTACTGCCGCGCGGCCGGAGCAGGGCATAGCCCTGGAGGTAAGGCCTAAACCCCCGCGCGGGGGAGCGCCGAAACATTGCCCCCTGGCGTTTACGGCCGc                                                >  6:332584/1‑135 (MQ=255)
                                                                                    cGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCGGATGCTGCCGCCCTGCCGGCTACCGCAATCGCCATTCACTAATGGCCGAACGCCCCGCCCCGGTGTGCGCCGAACCGTTGCAGCCGGGCGTTtccgtgccc                                                >  4:334010/1‑131 (MQ=255)
                                                                                              tCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCAACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCa                                        <  3:18452/133‑1 (MQ=255)
                                                                                               cGAAAGCTGAATGTGGATTGCCTGATCCTGCCGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAgg                                     <  5:319458/135‑1 (MQ=255)
                                                                                                   cgctgctTGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGAAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGgaga                                  <  4:374774/133‑1 (MQ=255)
                                                                                                        aaTGTGGATTGCCTGATCGTGCCGGTATTCATTTTCAGCAAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGc                            <  6:337446/134‑1 (MQ=255)
                                                                                                           gtgGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAAt                         <  4:346043/135‑1 (MQ=255)
                                                                                                                      gttCCTGCCGCTGTGATGTATAAGTCAATCGACATGCAGTAATTGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaac              <  3:316327/133‑1 (MQ=255)
                                                                                                                      gagccTGCCGCGCGGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAaac              <  3:222906/132‑1 (MQ=255)
                                                                                                                              cgcgCTGCTGGATCAGGAAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTc      <  5:332584/135‑1 (MQ=255)
                                                                                                                                cgcTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcac    >  2:316745/1‑135 (MQ=255)
                                                                                                                                  cTGCAGGATCCGTCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGAAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTcaccc  <  5:68282/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
AGTTCCAGGGTCGGTTCGGTGGACGTGGCGGTGGTCATGATGGCTTTCCTCTGCTGACCGACGCAGCATAGGTTTGGGCAAAGACGGGCGCTAGTCGAAAGCTGCATGTGGATTGCCTGATCCTGCCGCGCTGCAGGATCAGGCAATCGACATGCAGTAATGGCCTAACGCGCCGCGCGGGTGAGCGCCTAACCTTGGCAGCCTGGCTTTTCCGTCCGCCTGCTTCAAGGAGATTGTGCAATGACCCTTAAACAACCGGTCACCC  >  ACB122/1584912‑1585176

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: