Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 456,440 C→A 8.4% A4S (GCT→TCT)  ramA_1 ← (R)‑stereoselective amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122456,4400CA 8.4% 111.9 / 2.3 44A4S (GCT→TCT) ramA_1(R)‑stereoselective amidase
Reads supporting (aligned to +/- strand):  ref base C (13/27);  new base A (2/2);  total (15/29)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

TTGTAGCCGGACAGGAACATCTCCGGGCACACCAGCAACTGTGCGCCGCGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGTGGTTCATGATACCCCACACCGCATTC  >  ACB122/456307‑456572
                                                                                                                                     |                                                                                                                                    
ttgtAGCCGGACAGGAACATCTCCGGGCACACCAGCAACTGTGCGCCGCGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCCCATCCAGTGGCTTGGGTGCGCCCGGGTACAGcgcg                                                                                                                                     >  5:126879/1‑135 (MQ=255)
    aGCCGGACAGTAACATCTCCGGGCACACCAGCAACTGTGCGCCGCGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATg                                                                                                                                  <  6:36861/134‑1 (MQ=255)
    aGCCGGACAGGAACATCTCCGGGCACACCAGCAACTGTGCGCCGCGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATgc                                                                                                                                 <  8:364382/135‑1 (MQ=255)
    aGCCGGACAGGAACATCTCCGGGCACACCAGCAACTGTGCGCCGCGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATg                                                                                                                                  <  6:41779/134‑1 (MQ=255)
         gACAGGAACATCTCCGGGCACACCAGCAACTGGGCGCCGCGGGCGGCGGCCAACTGCGCCTGGGGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATg                                                                                                                            <  3:206546/135‑1 (MQ=255)
                     tcCGGGCACACCAGCTACTGGGCGCCGCGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGTTGCGCATGGTGCCTCCTTg                                                                                                                 <  7:376267/134‑1 (MQ=255)
                       cGGGCACACCAGCAACTGTGCGCCGCGGTCGGCTGCCAACATCGCCTGGTGGCGCCGCCGTTGCAGGTTGCCCGGCCCATCCAGTGGCTTGGGTGCGCCCTGGTTCAGAGCGATGCGCATGGGGCCTCCGTGtct                                                                                                              >  1:10953/1‑135 (MQ=255)
                        gggCACACCAGCAACTGTGCGCCGCGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCGGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGtc                                                                                                               <  4:128452/133‑1 (MQ=255)
                                 aGCAACTGTGCGCCGCGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCa                                                                                                    >  7:67721/1‑135 (MQ=255)
                                           cgccgcGGTCGGCTGCCAACTGCGCCTGGTGGCGCATCCCTTGCCGGTTGCCCGGGACATCCAGTGGCTTGGGTGCGCCCTGGTAAAGAGAGATGCGCATGGTTCCTACGTGTCTCAAAACGCCACGGCAATCgg                                                                                          >  4:118420/1‑135 (MQ=255)
                                           cgccgcGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGGCCGGTGCCGGGCCCCTCCAGTGGCTTGGGGGCGCCCGGGTACCGAGCGCTGCGGATGGTGCCTCCGTGTCTACATCCGCCAGGGCGGTCgg                                                                                          >  5:370965/1‑135 (MQ=255)
                                           cgccgcGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGCGCGCTGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCgg                                                                                          >  4:118694/1‑135 (MQ=255)
                                           cgccgcGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCgg                                                                                          >  7:256558/1‑135 (MQ=255)
                                           cgccgcGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCGTGTCTCAATCCGCCAGGGCGATCgg                                                                                          >  3:119997/1‑135 (MQ=255)
                                           cgccgcGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGAGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCgg                                                                                          >  8:265832/1‑135 (MQ=255)
                                                 ggTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTTCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGCGCGATTCGCATGGTGCCGCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGAt                                                                                    >  1:15058/1‑135 (MQ=255)
                                                 ggTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTTCAGCGCGAGGCGCATGGTGCCTGCGTGTCGCAATCCGCCAGGGCGCTCGGGCCGAt                                                                                    >  1:8173/1‑135 (MQ=255)
                                                               tGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGGTGAACACAt                                                                      <  3:121068/135‑1 (MQ=255)
                                                                  gcCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATc                                                                     <  2:17807/133‑1 (MQ=255)
                                                                  gcCGGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATc                                                                     <  4:201296/133‑1 (MQ=255)
                                                                  gcCGGGGGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATc                                                                     <  2:296083/133‑1 (MQ=255)
                                                                              aGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGt                                                       <  5:7767/135‑1 (MQ=255)
                                                                                   ttGGAGGTTGCCGGGCACATCCAGTGGCTGGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGGCTCAATCCGCCAGGGCGAGCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCg                                                  <  2:337689/135‑1 (MQ=255)
                                                                                           tGCGGGGCACATCGAGTGGCTTGGGGGCGCCCTGGTACAGAGCGAGGCGCAGGGTGCCTCCTTGTCTCAATCCGCCCGGGCGATCGGGCCGTTTTCGTTGATCACATCGCCCGGGCCGGGGTTGTCGgggtgggt                                          <  1:342410/135‑1 (MQ=255)
                                                                                           tGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGgggtgggt                                          <  6:271255/135‑1 (MQ=255)
                                                                                           tGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGgggggggt                                          <  3:175348/135‑1 (MQ=255)
                                                                                             ccGGGCACATCCTGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGgtgt                                        <  8:210079/135‑1 (MQ=255)
                                                                                              cGGGCACATCCAGTGGCTTGGTTGCGCCCTTGTACAGTGAGATGCGCTTGGTGCCTCCTTGTCTGAATCCGCCAGGGCTATCGGTCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGgtgt                                        <  3:118420/134‑1 (MQ=255)
                                                                                              cGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGc                                                                                                >  2:328246/1‑78 (MQ=255)
                                                                                              cGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGc                                                                                                <  1:328246/78‑1 (MQ=255)
                                                                                              cGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGAGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGgtgt                                        <  7:265832/134‑1 (MQ=255)
                                                                                                  cacaTCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGcca                                   <  3:356137/135‑1 (MQ=255)
                                                                                                     ttCCAGTGGCTGGGGGGCGCCCGGGTACAGAGCGAGGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGGGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGccacca                                <  5:240974/134‑1 (MQ=255)
                                                                                                     aTCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGccacca                                <  2:32350/135‑1 (MQ=255)
                                                                                                         aGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGGGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGt                            <  8:190208/135‑1 (MQ=255)
                                                                                                         aGTGGCTGGGGTGCGCCCGGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGt                            <  2:266256/135‑1 (MQ=255)
                                                                                                              cTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGTGGTTc                       >  6:62409/1‑135 (MQ=255)
                                                                                                                   gTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGTGGTTCATGAt                  >  4:205265/1‑135 (MQ=255)
                                                                                                                        cccTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGTGGTTCATGCTAcccc             >  5:367611/1‑135 (MQ=255)
                                                                                                                               aCAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTTGCCACCAAAGTGGTTCATGCTACCCCACACCGc      >  6:121946/1‑135 (MQ=255)
                                                                                                                                   cgCGATGGGCAGGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTGGGTTGCACACATGGCCCGGGCCGGGGTTGTCGGGGTGGGGGTGGCCACCAAAGTGGTTCATGATACCCCACACCGCATTc  <  4:141072/134‑1 (MQ=255)
                                                                                                                                   agCGATTCGCATGGTGCCTCCTGGTCTCAATCGGCCAGGGCGCTCGGGCCGATTTCGTTGAACACATCGCCCTGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGTGGTTCATGATACCCCACACCGCATTc  <  1:312194/135‑1 (MQ=255)
                                                                                                                                   agCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGTGGTTCATGATACCCCACACCGCATTc  <  6:92978/135‑1 (MQ=255)
                                                                                                                                   agCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGTGGTTCATGATACCCCACACCGCATTc  <  3:331950/135‑1 (MQ=255)
                                                                                                                                   agCGATGCGCATGGTGCCGCCTTGTCTCAATCCGCCCGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGTGGTTCATGATACCCCACACCGCATTc  <  8:325540/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
TTGTAGCCGGACAGGAACATCTCCGGGCACACCAGCAACTGTGCGCCGCGGTCGGCTGCCAACTGCGCCTGGTGGCGCAGCCGTTGCAGGTTGCCGGGCACATCCAGTGGCTTGGGTGCGCCCTGGTACAGAGCGATGCGCATGGTGCCTCCTTGTCTCAATCCGCCAGGGCGATCGGGCCGATTTCGTTGAACACATCGCCCGGGCCGGGGTTGTCGGGGTGGGTGTGGCCACCAAAGTGGTTCATGATACCCCACACCGCATTC  >  ACB122/456307‑456572

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