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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,083,3981.→G11.4% 14.4 ‑3.6 ‑0.535coding (768/1941 nt)yegIhypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (11/20):  total (15/20)
Fisher's exact test strand distribution p-value = 2.61e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.88e-01

CCCGGTAACATCGACAGCGGAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGCTGATTATCTGACGCGTA  >  REL606/2083363‑2083433
                                    |                                   
cccGGTAACATCGACAGCGGAATCGATCGTggcggc                                      >  1:6601218/1‑36 (MQ=255)
cccGGTAACATCGACAGCGGAATCGATCGTggcggc                                      <  1:4176231/36‑1 (MQ=255)
cccGGTAACATCGACAGCGGAATCGATCGTggcggc                                      <  1:2907011/36‑1 (MQ=255)
 ccGGTAACATCGACAGCGGAATCGATCGTGgcggc.g                                    <  1:4649820/36‑1 (MQ=255)
 ccGGTAACATCGACAGCGGAATCGATCGTGgcggc.g                                    >  1:7296124/1‑36 (MQ=255)
 ccGGTAACATCGACAGCGGAATCGATCGTGgcggc.g                                    >  1:314472/1‑36 (MQ=255)
 ccGGTAACATCGACAGCGGAATCGATCGTGgcggc.g                                    <  1:8389977/36‑1 (MQ=255)
  cGGTAACATCGACAGCGGAATCGATCGTGGcggc.gg                                   <  1:4678726/36‑1 (MQ=255)
   ggcaaCAGCGACAGCGGAATCGATCGTGGCGGC.GGt                                  <  1:5105222/33‑1 (MQ=37)
     tAACATCGACAGCGGAATCGATCGTGGCGGC.GGgtt                                >  1:6573309/1‑36 (MQ=255)
     tAAAATCGACAGCGGAATCGATCGTGGcggg.ggggt                                >  1:6357058/1‑33 (MQ=25)
      aaCATCGACAGCGGAATCGATCGTGGCGGC.GGtttt                               >  1:3204732/1‑36 (MQ=255)
      aaCATCGACAGCGGAATCGATCGTGGCGGC.GGtttt                               <  1:4073282/36‑1 (MQ=255)
       aCATCGACAGCGGAATCGATCGTGGCGGC.GGgtttt                              >  1:4537863/1‑35 (MQ=255)
       aCATCGACAGCGGAATCGATCGTGGCGGC.GGTTTTa                              <  1:5979194/36‑1 (MQ=255)
         aTCGACAGCGGAATCGATCGTGGCGGC.GGTTTTaaa                            <  1:703714/36‑1 (MQ=255)
          tCGACAGCGGAATCGATCGTGGCGGCGGGTTTTaaa                            >  1:4368741/1‑36 (MQ=38)
           cGACAGCGGAATCGATCGTGGCGGCGGGTTTTAAAc                           >  1:6019132/1‑36 (MQ=38)
           cGACAGCGGAATCGATCGTGGCGGC.GGTTTTAAAcc                          <  1:8010885/36‑1 (MQ=255)
            gACAGCGGAATCGATCGTGGCGGCGGGTTTTAAAca                          >  1:7013495/1‑35 (MQ=37)
            gACAGCGGAATCGATCGTGGCGGC.GGGTTTaaaacc                         >  1:7843365/1‑33 (MQ=38)
             aCAGCGGAATCGATCGTGGCGGC.GGGTTTAAACCAc                        >  1:5043089/1‑36 (MQ=255)
               aGCGGAATCGATCGTGGCGGC.GGTTTTAAACCACGt                      <  1:6224266/36‑1 (MQ=255)
                gCGGAATCGATCGTGGCGGC.GGTTTTAAACCACGTc                     <  1:7269665/36‑1 (MQ=255)
                gCGGAATCGATCGTGGCGGC.GGTTTTAAACCACGTc                     <  1:4886378/36‑1 (MQ=255)
                gCGGAATCGATCGTGGCGGC.GGTTTTAAACCACGTc                     <  1:4511629/36‑1 (MQ=255)
                 cGGAATCGATCGTGGCGGC.GGGTTTTAACCACGGCg                    >  1:2676383/1‑36 (MQ=25)
                 cGGAATCGATCGTGGCGGC.GGGTTTTAAACACGTCg                    >  1:565879/1‑36 (MQ=37)
                 cGGAATCGATCGTGGCGGC.GGGTTTTAAACACGGCg                    >  1:656100/1‑36 (MQ=19)
                  ggACTCGATCGTGGCGGC.GGTTTTAAACCACTTCGc                   <  1:4741841/36‑1 (MQ=37)
                  ggAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGc                   <  1:8358797/36‑1 (MQ=255)
                  ggAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGc                   <  1:3394965/36‑1 (MQ=255)
                   gAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGCt                  <  1:20040/36‑1 (MQ=255)
                   gAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGCt                  <  1:266814/36‑1 (MQ=255)
                     atcgatcgTGGCGGCGGGTTTTAAACC.CGTCGCtcg                >  1:7830989/1‑34 (MQ=19)
                       cgatcgTGGCGGC.GGTTTTAAACCACGTCGCTGAtt              <  1:6559128/36‑1 (MQ=255)
                        gatcgTGGCGGC.GGTTTTTAACCACGTCGCTGAttt             >  1:2056179/1‑35 (MQ=255)
                        gatcgTGGCGGC.GGGTTTTAAACACGTCGCTGGttt             >  1:4709002/1‑35 (MQ=19)
                         atcgTGGCGGC.GGTTTTAAACCACGTCGCTGATTAt            <  1:5050659/36‑1 (MQ=255)
                         atcgTGGCGGC.GGGTTTTAACCACGTCGCTGATTAt            >  1:6651152/1‑36 (MQ=25)
                          tcgTGGCGGC.GGTTTTAAACCACGTCGCTGATTATc           <  1:4615676/36‑1 (MQ=255)
                           cgTGGCGGC.GGTTTTAAACCACGTCGCTGATTATCt          <  1:8002616/36‑1 (MQ=255)
                             tGGCGGC.GGTTTTAAACCACGTCGCTGATTATCTGa        <  1:2412489/36‑1 (MQ=255)
                             tGGCGGC.GGTTTTAAACCACGTCGCTGATTATCTGa        <  1:1033624/36‑1 (MQ=255)
                              ggcggc.GGTTTTAAACCACGTCGCTGATTATCTGAc       <  1:7280164/36‑1 (MQ=255)
                               gcggc.gGTTTTAAACCACGTCGCTGATTATCTGAcg      <  1:4032831/36‑1 (MQ=255)
                                 ggc.ggTTTTAAACCACGTCGCTGGTTATCTGAcgcg    >  1:3794995/1‑36 (MQ=255)
                                 ggc.ggTTTTAAACCACGTCGCTGATTATCTGAcgcg    <  1:6973996/36‑1 (MQ=255)
                                 ggc.ggTTTTAAACCACGTCGCTGATTATCTGAcgcg    >  1:3155655/1‑36 (MQ=255)
                                  gc.ggTTTTAAACCACGTCGCTGATTATCTGACGCGt   <  1:7064390/36‑1 (MQ=255)
                                   c.ggTTTTAAACCACGTCGCTGATTATCTGACGCGTa  >  1:7618060/1‑36 (MQ=255)
                                    |                                   
CCCGGTAACATCGACAGCGGAATCGATCGTGGCGGC.GGTTTTAAACCACGTCGCTGATTATCTGACGCGTA  >  REL606/2083363‑2083433

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG