Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,004,497 G→A 9.1% W20* (TGG→TGA)  rhaR_6 → HTH‑type transcriptional activator RhaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,004,4970GA 9.1% 135.3 / 3.7 44W20* (TGG→TGA) rhaR_6HTH‑type transcriptional activator RhaR
Reads supporting (aligned to +/- strand):  ref base G (5/35);  new base A (2/2);  total (7/37)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTGTGCTATTCCATAAAAAACCTGCACTTTTACGAAAAGCGTTCGATTAACGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAGCTGCTGTATGTACAGCGTGGCCAG  >  ACB122/4004368‑4004617
                                                                                                                                 |                                                                                                                        
aaGTGTGCTATTCCATAAAAAACCTGCACTTTTACGAAAAGCGTTCGATTAACGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCg                                                                                                                      <  4:83017/134‑1 (MQ=255)
    gtgCTATTCCATAAAAAACCTGCACTTTTACGAAAAGCGTTCGATTAACGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCAcc                                                                                                                 <  3:83245/135‑1 (MQ=255)
          ttCCATAAAAAACCTGCACTTTTACGAAAAGCGTTCGATTAACGGACAGGCCAACCCCAGATGAAATCCACCATCCCCGGCGTCCCGTTGTTAACGCGGGATGGCGAGAATCAGGCCTGGCCCGGCCCCGACctg                                                                                                           >  5:167758/1‑135 (MQ=255)
               tAAAAAACCTGCACTTTTACGAAAAGCGTTCGATTAACGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCa                                                                                                      <  5:160081/135‑1 (MQ=255)
                       cTGCACTTTTACGAAAAGCGTTCGATTCACGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGt                                                                                              <  6:209702/135‑1 (MQ=255)
                          cACTTTTACGAACAGCGTTCGATTAACGGACAGGCCAACCCCAGATGCAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGa                                                                                           <  8:183917/135‑1 (MQ=255)
                               ttACGAAACGCGTTCGATTAACGGACAGGCCACCCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCTATCCCg                                                                                      <  7:147789/135‑1 (MQ=255)
                                      aaGCGTTCGATTAACGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCa                                                                               <  1:344142/135‑1 (MQ=255)
                                               aTTAACGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTg                                                                       <  5:362/134‑1 (MQ=255)
                                               aTTAACGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGc                                                                      <  3:256504/135‑1 (MQ=255)
                                                    cGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCa                                                                  <  5:57933/134‑1 (MQ=255)
                                                    cGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCAc                                                                 <  5:94924/135‑1 (MQ=255)
                                                        cAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGga                                                            <  7:58834/136‑1 (MQ=255)
                                                         aGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGga                                                            <  4:53232/135‑1 (MQ=255)
                                                               aCCCCAGATGAAATCCACCCTCGCCGGCGTCGCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCTCTGGGAGATCAAg                                                     <  7:208270/136‑1 (MQ=255)
                                                                ccccAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTCTGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAg                                                     <  3:275574/135‑1 (MQ=255)
                                                                         aaaTCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGc                                           <  8:210429/136‑1 (MQ=255)
                                                                               aCCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGcc                                      <  4:207051/135‑1 (MQ=255)
                                                                               aCCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGcc                                      <  5:248375/135‑1 (MQ=255)
                                                                               aCCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGcc                                      <  1:176761/135‑1 (MQ=255)
                                                                                       cGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTc                             >  7:277542/1‑136 (MQ=255)
                                                                                       cGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTc                             <  7:368705/136‑1 (MQ=255)
                                                                                        ggCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTc                             <  5:65319/135‑1 (MQ=255)
                                                                                         gCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGGACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCa                           <  8:76163/136‑1 (MQ=255)
                                                                                          cGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCa                           <  4:95448/135‑1 (MQ=255)
                                                                                          cGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCa                           <  6:137179/135‑1 (MQ=255)
                                                                                          cGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCa                           <  6:167758/135‑1 (MQ=255)
                                                                                           gTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGcc                                      <  8:127403/123‑1 (MQ=255)
                                                                                           gTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGcc                                      >  7:127403/1‑123 (MQ=255)
                                                                                            tCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAgct                        <  7:178768/136‑1 (MQ=255)
                                                                                            tCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAgc                         <  4:201615/135‑1 (MQ=255)
                                                                                               cGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAgctgc                      <  6:344361/135‑1 (MQ=255)
                                                                                               cGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAgctgc                      <  3:127431/135‑1 (MQ=255)
                                                                                                        aGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAGCTGCTGTATGTACa            <  7:14191/136‑1 (MQ=255)
                                                                                                        aGCTGGATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAGCTGCTGTATGTAc             <  2:7317/135‑1 (MQ=255)
                                                                                                             tATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAGCTGCTGTATGTACAGCGt        <  8:59383/135‑1 (MQ=255)
                                                                                                              aTGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACT‑GGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAGCTGCTGTATGTACAGCGt        <  3:63228/133‑1 (MQ=255)
                                                                                                                 gCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAGCTGCTGTATGTACAGCGTgg      >  3:187526/1‑133 (MQ=255)
                                                                                                                 gCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACGGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAGCTGCTGTATGTACAGCGTgg      >  6:249173/1‑133 (MQ=255)
                                                                                                                 gCGAGAATCAGGCCTGACCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAgc                         >  2:3502/1‑114 (MQ=255)
                                                                                                                 gCGAGAATCAGGCCTGACCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAgc                         >  2:10753/1‑114 (MQ=255)
                                                                                                                 gCGAGAATCAGGCCTGACCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAgc                         <  1:3502/114‑1 (MQ=255)
                                                                                                                 gCGAGAATCAGGCCTGACCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAgc                         <  1:10753/114‑1 (MQ=255)
                                                                                                                   gagaATCCGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGGCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGCTCCAGCTGCTGTATGTACAGCGTGGCCAg  <  5:56661/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                        
AAGTGTGCTATTCCATAAAAAACCTGCACTTTTACGAAAAGCGTTCGATTAACGGACAGGCCAACCCCAGATGAAATCCACCCTCCCCGGCGTCCCGTTGTTCAAGCTGTATGGCGAGAATCAGGCCTGGCCCGGCACCGACCTGCTGCACTGCGAGTCGATCCCGGCACGCAGCCGCCTGCATCACTGGGAGATCAAGCCGCACCAGCATGCCGAGCTGTTCCAGCTGCTGTATGTACAGCGTGGCCAG  >  ACB122/4004368‑4004617

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