Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,516,335 A→T 14.2% intergenic (+166/+417) yjhR → / ← yjhS frameshift suppressor contains phospolipase D/nuclease; b4308_2/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,516,3350A→T14.2% 7.3 ‑4.4 0.039intergenic (+166/+417)yjhR/yjhSframeshift suppressor contains phospolipase D/nuclease; b4308_2/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/0):  ref base (12/21):  total (18/21)
Fisher's exact test strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAAT  >  REL606/4516300‑4516370
                                   |                                   
gCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTa                                     <  1:6761036/36‑1 (MQ=255)
 cATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAg                                    <  1:6104447/36‑1 (MQ=255)
 cATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAg                                    <  1:3785844/36‑1 (MQ=255)
  aTAATGCGCCTGGGATTAATGGCGTGTGCTTTTAgg                                   <  1:6785474/36‑1 (MQ=255)
   tAATGCGCCTGGGGTTAATGGCGGGGGGTTTTGGGa                                  >  1:640053/1‑36 (MQ=16)
   tAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGa                                  <  1:3257740/36‑1 (MQ=255)
   tAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGa                                  <  1:8291050/36‑1 (MQ=255)
    aaTGCGCCTGGGATTAATGGCGGGTGGTTTTTggga                                 >  1:3380735/1‑34 (MQ=25)
     aTGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAAt                                <  1:5646642/36‑1 (MQ=255)
       gcgcCTGGGATTAATGGCGGGTGCTTTTAGGAATga                              <  1:4033675/36‑1 (MQ=255)
       gcgcCTGGGATTAATGGCGGGGGCTTTTTGGAGTga                              >  1:1693529/1‑36 (MQ=25)
        cgcCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                             <  1:5324734/36‑1 (MQ=255)
        cgcCTGGGATTAATGGCGGGTGCTTTTAGGAATgaa                             <  1:4082173/36‑1 (MQ=255)
        cgcCTGGGATTAATGGCGGGGGGTTTTTGGGATgta                             >  1:5761154/1‑34 (MQ=16)
            tGGGATTAATGGCGGGGGGTTTTTTGGATtaaagaa                         >  1:2355912/1‑36 (MQ=9)
             gggATTAATGGCGGGTGCTTTTAGGAATgaaagaaa                        <  1:571439/36‑1 (MQ=255)
              ggATTAATGGCGGGTGCTTTTAGGAATGAAAGaaaa                       <  1:4373006/36‑1 (MQ=255)
               gATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAAt                      <  1:6638744/36‑1 (MQ=255)
                aTTAATGGCGGGTGGTTTTTTGGATTAAAGaaaaag                     >  1:7491346/1‑34 (MQ=12)
                 ttAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgc                    <  1:5365818/36‑1 (MQ=255)
                 ttAATGGCGGGTGCTTTTAGGAATGAAAGAAAATgc                    <  1:2122208/36‑1 (MQ=255)
                     tGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagc                >  1:4003164/1‑36 (MQ=255)
                     tGGCGGGTGCTTTTAGGAATGAAAGAAAATgcaagc                <  1:5711909/36‑1 (MQ=255)
                      ggCGGGTGCTTTTAGGAATGAAAGAAAATgcaagca               >  1:8342974/1‑36 (MQ=255)
                      ggCGGGTGCTTTTAGGAATGAAAGAAAATgcaagca               <  1:7465955/36‑1 (MQ=255)
                        cGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCaat             >  1:4187128/1‑36 (MQ=255)
                         gggTGCTTTTAGGAATGAAAGAAAATGCAAGCaata            >  1:7343854/1‑36 (MQ=255)
                         gggTGCTTTTAGGAATGAAAGAAAATGCAAGCaata            <  1:7614798/36‑1 (MQ=255)
                          ggTGCTTTTAGGAATGAAAGAAAATGCAAGCaataa           >  1:8086601/1‑36 (MQ=255)
                          ggTGATTTTAGGAATGAAAGAAAATGCAAGCaataa           <  1:1955637/36‑1 (MQ=255)
                           gTGCTTTTAGGAATGAAAGAAAATGCAAGCaataat          >  1:1601437/1‑36 (MQ=255)
                            tGCTTTTAGGAGTGAAAGAAAATGCAAGCAataata         >  1:1459309/1‑36 (MQ=255)
                            tGCTTTTAGGAATGAAAGAAAATGCAAGCAataata         >  1:879112/1‑36 (MQ=255)
                             gCTTTTAGGAATGAAAGAAAATGCAAGCAataatat        >  1:2453872/1‑36 (MQ=255)
                              cTTTTAGGAATGAAAGAAAATGCCAGCAataatata       <  1:6970931/36‑1 (MQ=255)
                              cTTTTAGGAATGAAAGAAAATGCAAGCAataatata       <  1:559840/36‑1 (MQ=255)
                              cTTTTAGGAATGAAAGAAAATGCAAGCAataatata       >  1:5211939/1‑36 (MQ=255)
                              cTTTTAGGAATGAAAGAAAATGCAAGCAataatata       >  1:3602493/1‑36 (MQ=255)
                                tttAGGAATGAAAGAAAATGCAAGCAataatataat     >  1:3299922/1‑36 (MQ=255)
                                 ttAGGAATGAAAGAAAATGCAAGCAAtaatataata    <  1:7121318/36‑1 (MQ=255)
                                   aGGAATGAAAGAAAATGCAAGCAATAATATaataat  <  1:6215798/36‑1 (MQ=255)
                                   |                                   
GCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAAT  >  REL606/4516300‑4516370

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG