Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,208,669 T→C 27.0% N296D (AAC→GAC)  yeiN ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,208,6690T→C27.0% 2.5 ‑2.4 ‑0.017886yeiNhypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (5/7):  total (10/7)
Fisher's exact test strand distribution p-value = 8.85e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.59e-01

GTATTCTTTGGCAATTTCACTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATATTGGATTTCAGGCTGT  >  REL606/2208634‑2208704
                                   |                                   
gTATTCTTTGGCAATTTCACTCGCCAGAATGGCgtt                                     >  1:177812/1‑36 (MQ=255)
gTATTCTTTGGCAATTTCACTCGCCAGAATGGCgtt                                     >  1:5477703/1‑36 (MQ=255)
   ttCTTTGGCAATTTCACTCGCCAGAATGGCgttgtt                                  >  1:103330/1‑36 (MQ=255)
   ttCTTTGGCAATTTCACTCGCCAGAATGGCgttgtt                                  <  1:357603/36‑1 (MQ=255)
   ttCTTTGGCAATTTCACTCGCCAGAATGGCgttgtt                                  <  1:228178/36‑1 (MQ=255)
        tGGCAATTTCACTCGCCAGAATGGCGTTGTTGAaca                             >  1:944691/1‑36 (MQ=255)
         ggCAATTTCACTCGCCAGAATGGCGTTGTTTAacac                            >  1:2771435/1‑36 (MQ=255)
            aaTTTCACTCGCCAGAATGGCGTTGTTGAACACCAg                         >  1:672657/1‑36 (MQ=255)
              tttCACTCGCCAGAATGGCGTTGTTGAACACCAGCt                       <  1:3190777/36‑1 (MQ=255)
               ttCACTCGCCAGAATGGCGTTGTTGAACACCAGCTg                      >  1:2996512/1‑36 (MQ=255)
               ttCACTCGCCAGAATGGCGTTGTTGAACACCAGCTg                      <  1:598567/36‑1 (MQ=255)
                tCACTCGCCAGAACAGCGTCGTTGAACACCAGCTgg                     >  1:1921548/1‑36 (MQ=25)
                  aCTCGCCAGAATAGCGTAGTTGAACACAAGCTgggt                   >  1:1078004/1‑34 (MQ=21)
                   cTCGCCAGAATGGCGCCGTTGAACACCAGCTGGata                  >  1:2812260/1‑36 (MQ=37)
                    tCGCCAGAATGGCGTTGTTGAACACCAGCTGGatat                 <  1:4407326/36‑1 (MQ=255)
                      gCCAGAATGGCGTTGTTGAACACCAGCTGGATATTg               <  1:2671849/36‑1 (MQ=255)
                      gCCAGAATGGCGCCGTTGAACACCAGCTGGATATTg               >  1:3796012/1‑36 (MQ=37)
                        cAGAATGGCGTTGTTGAACACCAGCTGGATATTGGa             <  1:1920365/36‑1 (MQ=255)
                        cAGAATGGCGTTGGTGAACACCAGCTGGATATTGGa             >  1:5843722/1‑36 (MQ=255)
                          gAATGGCGTTGTTGAACACCAGCTGGATATTGGAtt           <  1:4476225/36‑1 (MQ=255)
                            aTGGCGTCGCTGAACACCAGCCGGATATTGGATTTc         >  1:687044/1‑36 (MQ=25)
                               gCGTTGTTGAACACCAGCTGGATATTGGATTTCAgg      <  1:1789623/36‑1 (MQ=255)
                                 gttgttGAACACCAGCTGGATATTGGATTTCAGGCt    >  1:742033/1‑36 (MQ=255)
                                  ttgttgAACACCAGCTGGATATTGGATTTCAGGCTg   >  1:203059/1‑36 (MQ=255)
                                   tgttgAACACCAGCTGGATATTGGATTTCAGGCTGt  >  1:3821929/1‑36 (MQ=255)
                                   |                                   
GTATTCTTTGGCAATTTCACTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATATTGGATTTCAGGCTGT  >  REL606/2208634‑2208704

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG