Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,029,304 G→T 30.1% intergenic (‑278/‑340) yihA ← / → yihI GTP‑binding protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,029,3040G→T30.1% 13.3 ‑2.1 0.038intergenic (‑278/‑340)yihA/yihIGTP‑binding protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/1):  ref base (16/10):  total (27/11)
Fisher's exact test strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAAATGTGCAAATACTGATGGCGG.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGCACCCATAGAGC  >  REL606/4029270‑4029339
                                   |                                   
aaaaTGTGCAAATACTGATGGCGG.TTGAttttttgt                                    >  1:7659368/1‑36 (MQ=255)
aaaaTGTGCAAATACTGATGGCGG.TTGAttgtttgt                                    <  1:1041272/36‑1 (MQ=255)
aaaaTGTGCAAATACTGATGGCGG.TTGAttgtttgt                                    <  1:6290414/36‑1 (MQ=255)
 aaaTGTGCAAATACTGATGGCGG.TTGATtgtttgtt                                   <  1:4361190/36‑1 (MQ=255)
  aaTGTGCAAATACTGATGGCGG.TTTATTtttttttt                                  >  1:3010249/1‑36 (MQ=21)
   aTGTGCAAATACTGATGGCGG.TTGATTGTTTGTTTa                                 <  1:6733239/36‑1 (MQ=255)
    tgtgCAAATACTGATGGCGG.TTGATTGTTTGTTTaa                                <  1:5572480/36‑1 (MQ=255)
     gtgCAAATACTGATGGCGG.TTTATTGTTTTTTTaaa                               >  1:7668720/1‑36 (MQ=37)
     gtgCAAATACTCATGGCGT.TTGATTGTTTTTTGaac                               <  1:5659920/36‑2 (MQ=16)
      tgCAAATACTGATGGCGG.TTTATTTTTTTTTTTAAg                              >  1:4304700/1‑36 (MQ=16)
      tgCAAATACTGATGGCGG.TTTATTGTTTTTTTTAAg                              >  1:1085452/1‑36 (MQ=21)
       gCAAATACTGATGGCGG.TTGATTGTTTGTTTAAAGc                             <  1:7352806/36‑1 (MQ=255)
       gCAAATACTGATGGCGG.TTGATTGTTTGTTTAAAGc                             <  1:3660253/36‑1 (MQ=255)
        cAAATACTGATGGCGG.TTGATTGTTTGTTTAAAGCa                            <  1:184349/36‑1 (MQ=255)
          aaTACTGATGGCGG.TTGATTTTTTTTTTTAATCaaa                          >  1:2694889/1‑36 (MQ=19)
          aaTACTGATGGCGG.TTGATTGTTTTTTTAAAGCaaa                          >  1:2548598/1‑36 (MQ=255)
           aTACTGATGGCGG.TTGGTTTTTTTTTTATAGCaaaa                         >  1:1272745/1‑35 (MQ=16)
           aTACTGATGGCGG.TTGATTGTTTTTTTTAAGCAAAg                         >  1:5782836/1‑36 (MQ=37)
            tACTGATGGCGGTTTGTTTTTTTTTTTAAAAAAAAg                         >  1:3822480/1‑36 (MQ=0)
            tACTGATGGCGG.TTGATTTTTTTTTTTAAGCaaaag                        >  1:6163921/1‑34 (MQ=21)
             aCTGATGGCTG.TTGATTGTTTGTTTAAAGCAAAGGc                       <  1:5803415/36‑1 (MQ=255)
             aCTGATGGCGG.TTTATTGTTTGTTTAAAGCAAAGGc                       >  1:6530605/1‑36 (MQ=255)
             aCTGATGGCGG.TTGATTGTTTTTTTTAAGCAAAggg                       >  1:2455429/1‑35 (MQ=37)
             aCTGATGGCGG.TTGATTGTTTGTTTAAAGCAAAggg                       >  1:910484/1‑35 (MQ=255)
                 aTGGCGG.TTGATTGTTTGTTTAAAGCAAAGGCGTaa                   <  1:2579094/36‑1 (MQ=255)
                  tGGCGG.TTGATTGTTTTTTTAAAGAAAAGGCGTaaa                  >  1:6264189/1‑36 (MQ=37)
                    gCGG.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGt                <  1:7798097/36‑1 (MQ=255)
                    gCGG.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGt                >  1:1306976/1‑36 (MQ=255)
                    gCGG.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGt                >  1:2875771/1‑36 (MQ=255)
                     cGG.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTa               >  1:873333/1‑36 (MQ=255)
                      gg.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAg              >  1:1845979/1‑36 (MQ=255)
                      gg.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAg              >  1:4303655/1‑36 (MQ=255)
                       g.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGc             >  1:1918588/1‑36 (MQ=255)
                       g.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGc             >  1:4880400/1‑36 (MQ=255)
                       g.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGc             <  1:5147772/36‑1 (MQ=255)
                       g.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGc             >  1:1060441/1‑36 (MQ=255)
                       g.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGc             >  1:2697786/1‑36 (MQ=255)
                         ttGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGCa            >  1:7841869/1‑36 (MQ=255)
                         ttGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGCa            >  1:1867527/1‑36 (MQ=255)
                          tGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGCAc           >  1:1488081/1‑36 (MQ=255)
                            aTTGTTTGTTTAAAGCAAAGGCGTAAAGTAGCAccc         >  1:5367846/1‑36 (MQ=255)
                            aTTGTTTGTTTAAAGCAAAGGCGTAAAGTAGCAccc         <  1:832800/36‑1 (MQ=255)
                            aTTGTTTGTTTAAAGCAAAGGCGTAAAGTAGCAccc         >  1:2934757/1‑36 (MQ=255)
                              tgtttgttTAAAGCAAAGGCGTAAAGTAGCACCCAt       <  1:1286377/36‑1 (MQ=255)
                               gtttgtttAAAGCAAAGGCGTAAAGTAGCACCCATa      <  1:6743195/36‑1 (MQ=255)
                               gtttgtttAAAGCAAAGGCGTAAAGTAGCACCCATa      <  1:6305328/36‑1 (MQ=255)
                               gtttgtttAAAGCAAAGGCGTAAAGTAGCACCCATa      <  1:2652453/36‑1 (MQ=255)
                               gtttgtttAAAGCAAAGGCGTAAAGTAGCACCCATa      >  1:4810127/1‑36 (MQ=255)
                                 ttgtttAAAGCAAAGGCGTAAAGTAGCACCCATaga    >  1:1237895/1‑36 (MQ=255)
                                 ttgtttAAAGCAAAGGCGTAAAGTAGCACCCATaga    <  1:3759228/36‑1 (MQ=255)
                                  tgtttAAAGCAAAGGCGTAAAGTAGCACCCATagag   >  1:7462853/1‑36 (MQ=255)
                                  tgtttAAAGCAAAGGCGTAAAGTAGCACCCATagag   >  1:1630776/1‑36 (MQ=255)
                                   gtttAAAGCAAAGGCGTAAAGTAGCACCCATAGAGc  <  1:2000064/36‑1 (MQ=255)
                                   gtttAAAGCAAAGGCGTAAAGTAGCACCCATAGAGc  <  1:5636354/36‑1 (MQ=255)
                                   gtttAAAGCAAAGGCGTAAAGTAGCACCCATAGAGc  >  1:4893897/1‑36 (MQ=255)
                                   gtttAAAGCAAAGGCGTAAAGTAGCACCCATAGAGc  >  1:1364327/1‑36 (MQ=255)
                                   |                                   
AAAATGTGCAAATACTGATGGCGG.TTGATTGTTTGTTTAAAGCAAAGGCGTAAAGTAGCACCCATAGAGC  >  REL606/4029270‑4029339

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG