Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,083,398 +G 5.3% coding (768/1941 nt) yegI ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,083,3981.→G 5.3% 4.0 ‑1.9 ‑0.238coding (768/1941 nt)yegIhypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (13/23):  total (15/23)
Fisher's exact test strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.64e-01

CCCGGTAACATCGACAGCGGAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGCTGATTATCTGACGCGTA  >  REL606/2083363‑2083433
                                    |                                   
cccGGTAACATCGACAGCGGAATCGATCGTggcggc                                      <  1:4893614/36‑1 (MQ=255)
 ccGGTAACATCGACAGCGGAATCGATCGTGgcggc.g                                    <  1:2269959/36‑1 (MQ=255)
   ggTAACATCGACAGCGGAATCGATCGTGGcggg.ggg                                  >  1:6968247/1‑35 (MQ=255)
   ggTAACATCGACAGCGGAATCGATCGTGGCGGC.GGt                                  <  1:6135567/36‑1 (MQ=255)
    gTAACATCGACAGCGGAATCGATCGTGGcggg.gggt                                 >  1:5055334/1‑34 (MQ=38)
     tAACATCGCCAGCGGAATCGATTGTGGCGGC.GGttt                                <  1:7492568/36‑1 (MQ=37)
     tAACATCGACAGCGGAATCGATCGTGGCGGC.GGttt                                <  1:6857347/36‑1 (MQ=255)
       aCATCGACAGCGGAATCGATCGTGGCGGC.GGgtttt                              >  1:3714334/1‑35 (MQ=255)
       aCATCGACAGCGGAATCGATCGTGGCGGC.GGTTTTa                              >  1:6446874/1‑36 (MQ=255)
       aCATCGACAGCGGAACCGATCGTGGCGGC.GGTTTTa                              <  1:3360845/36‑1 (MQ=255)
        cATCGACAGCGGAATCGATCGTGGCGGC.GGTTTTaa                             <  1:4674592/36‑1 (MQ=255)
        cATCAACAGCGGAATCGATCGTGGCGGC.GGTTTTaa                             <  1:6752354/36‑1 (MQ=255)
         aTCGACAGCGGAATCGATCGTGGCGGC.GGTTTTaaa                            <  1:7628326/36‑1 (MQ=255)
         aTCGACAGCGGAATCGATCGTGGCGGC.GGTTTTaaa                            <  1:2120205/36‑1 (MQ=255)
         aTCGACAGCGGAATCGATAGAGGCGGC.GGgttttga                            >  1:5293492/1‑33 (MQ=19)
          ccGACAGCGGAATCGATCGTGGCGGC.GGTTTTAAAc                           <  1:4474798/35‑1 (MQ=255)
           cGACAGCGGAATCGATCGTGGCGGC.GGTTTTAAAcc                          <  1:4238231/36‑1 (MQ=255)
           cGACAGCGGAATCGATCGTGGCGGC.GGTTTTAAAcc                          <  1:7753261/36‑1 (MQ=255)
            gACAGCGGAATCGATCGTGGCGGCGGGTTTTaaaaa                          >  1:1388140/1‑34 (MQ=37)
             aCAGCGGAATCGATCGTGGCGGC.GGTTTTAAACCAc                        <  1:6137186/36‑1 (MQ=255)
             aCAGCGGAATCGATCGTGGCGGC.GGTTTTAAACCAc                        <  1:3576149/36‑1 (MQ=255)
              cAGCGGAATCGATCGTGGCGGCGGGTTTTAAAcctg                        >  1:1266575/1‑34 (MQ=37)
              cAGCGGAATCGATCGTGGCGGC.GGGTTTTAACCCCg                       >  1:6757559/1‑36 (MQ=25)
               aGCGGAATCGATCGTGGCGGC.GGTTTTAAACCACGt                      >  1:7406533/1‑36 (MQ=255)
               aGCGGAATCGATCGTCGCGGC.AGTTTTTAACCACGt                      >  1:10969/1‑36 (MQ=25)
                gCGGAATCGATCGTGGCGGC.GGGTTTTAAACACggg                     >  1:5644907/1‑34 (MQ=25)
                 cGGAATCGATCGTGGCGGC.GGTTTTAAACCACGTCg                    <  1:5010969/36‑1 (MQ=255)
                 cGGAATCGATCGTGGCGGC.GGTTTTAAACCACGTCg                    <  1:3935135/36‑1 (MQ=255)
                  ggAATCGGTCGCGGCGGC.GGTTTTAAACCACGTCGc                   <  1:180126/36‑1 (MQ=37)
                  ggAATCGATCGTGGCGGC.GGTTTTAAACCACCTCGc                   >  1:7050503/1‑36 (MQ=255)
                  ggAATCGATCGTGGCGGC.GGGTTTAAAACACGTCGc                   >  1:7897570/1‑36 (MQ=38)
                   gAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGCt                  <  1:5618131/36‑1 (MQ=255)
                   gAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGCt                  <  1:288325/36‑1 (MQ=255)
                    aaTCGATCGTGGCGGC.GGTTTTAAACCACGTCGCTg                 <  1:877803/36‑1 (MQ=255)
                    aaTCGATCGTGGCGGC.GGGTTTTAAACACGTCGCTg                 >  1:4200507/1‑36 (MQ=37)
                      tcgatcgTGGCGGC.GGTTTTTAACCACGTCGCTGAt               >  1:208044/1‑36 (MQ=255)
                         atcgTGGCGGA.GGGTTTTAAACACCTCGCTTAtttt            >  1:3629961/1‑34 (MQ=2)
                          tcgTGGCGGC.GGTTTTAAACCACGTCGCTGATTATc           <  1:3825631/36‑1 (MQ=255)
                           cgTGGCGGC.GGTTTTAAACCACGTCGCTGATTATCt          <  1:225544/36‑1 (MQ=255)
                            gTGGCGGC.GGTTTTAAACCACGTCGCTGATTATCTg         <  1:3605415/36‑1 (MQ=255)
                            gTGGCGGC.GGTTTTAAACCACGTCGCTGATTATCTg         <  1:1020737/36‑1 (MQ=255)
                             tGGCGGC.GGTTTTAAACCACGTCGCTGATTTTCtta        >  1:7709266/1‑34 (MQ=37)
                              gacggc.tgTTTTAAACCACGTCGCTTATTATCTGAc       <  1:3426234/34‑1 (MQ=37)
                               gcggc.gGTTTTAAACCACGTCGCTGATTATCTGAcg      <  1:7089420/36‑1 (MQ=255)
                                 ggc.ggTTTTAAACCACGTCGCTGATTATCTTAcgcg    >  1:5689950/1‑36 (MQ=255)
                                  gc.ggTTTTAAACCACGTCGCTGATTATCTGACGCGt   >  1:4811999/1‑36 (MQ=255)
                                   c.ggTTTTAAACCACGTCGCTGATTATCTGACGCGTa  <  1:4589307/36‑1 (MQ=255)
                                    |                                   
CCCGGTAACATCGACAGCGGAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGCTGATTATCTGACGCGTA  >  REL606/2083363‑2083433

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 15 ≤  ATCG  < 28 ≤  ATCG  < 37 ≤  ATCG