Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,594,508 T→A 17.0% intergenic (+41/+7) yphA → / ← yphB predicted inner membrane protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,594,5080T→A17.0% 2.5 ‑3.7 0.023intergenic (+41/+7)yphA/yphBpredicted inner membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (13/6):  total (13/10)
Fisher's exact test strand distribution p-value = 2.37e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCT  >  REL606/2594473‑2594543
                                   |                                   
aTGCCGGATGCGGCGTGAACGCCTTATCCGGCAttt                                     >  1:2012612/1‑36 (MQ=35)
aTGCCGGATGCGGCGTGAACGCCTTATCCGGCAttt                                     <  1:3460983/36‑1 (MQ=35)
  gCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                   >  1:1097544/1‑36 (MQ=40)
   ccGGATGCGGCGTGAGCGCCTTATCCGGCATTtaat                                  <  1:4335797/36‑1 (MQ=255)
   ccGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                  <  1:7903699/36‑1 (MQ=255)
   ccGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                  <  1:3979860/36‑1 (MQ=255)
   ccGGATGCGGCGTGAAAGCCTTATCCGGCATTtaat                                  >  1:4637177/1‑36 (MQ=255)
    cGGATGCGGCGTGAACGCCTTATCCGGCATTtaata                                 <  1:7289237/36‑1 (MQ=255)
    cGGATGCGGCGTGAACGCCTTATCCGGCATTtaata                                 <  1:2843212/36‑1 (MQ=255)
    aggATGCGGCGTGAACGCCTTATCCGGCATTtaata                                 <  1:3326323/35‑1 (MQ=255)
     ggATGCGGCGTGAACGCCTTATCCGGCATTtaataa                                <  1:8365591/36‑1 (MQ=255)
     ggATGCGGCGTGAACGCCTTATCCGGCATTtaataa                                >  1:6070593/1‑36 (MQ=255)
     ggATGCGGCGTGAACGCCTTATCCGGCATTtaataa                                >  1:4484193/1‑36 (MQ=255)
      gATGCGGCGTGAACGCCTTATCCGGCATTTAATaaa                               <  1:639075/36‑1 (MQ=255)
       aTGCGGCGTGAACGCCTTATCCGGCATTTAATAAAt                              >  1:442135/1‑36 (MQ=255)
       aTGCGGCGTGAACGCCTTATCCGGCATTTAATAAAt                              <  1:3383299/36‑1 (MQ=255)
         gcggcgTGAACGCCTTATCCGGCATTTAATAAATTa                            >  1:3658003/1‑36 (MQ=255)
         gcggcgTGAACGCCTTATCCGGCATTTAATAAATTa                            >  1:5554971/1‑36 (MQ=255)
         gcggcaTGAACTCCTTATCCTGCATTTAATAAATTa                            <  1:2567976/36‑1 (MQ=21)
          cggcgTGAACGCCTTATCCGGCATTTAATAAATTAc                           >  1:4456705/1‑36 (MQ=255)
           ggcgTGAACGCCTTATCCGGCATTTAATAAATTACa                          <  1:379462/36‑1 (MQ=255)
            gcgTGAACGCCTTATCCGGCATTTAATAAATTACAg                         <  1:3602584/36‑1 (MQ=255)
            gcgTGAACGCCTTATCCGGCATTTAATAAATTACAg                         >  1:5126434/1‑36 (MQ=255)
               tGAACGCCTTATCCGGCATTTAATAAATTACAGCCa                      <  1:7958338/36‑1 (MQ=255)
                gAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                     >  1:7326775/1‑36 (MQ=255)
                    gCCTTATCCGGCATTTAATAAATTACAGCCACTCTa                 >  1:4397600/1‑36 (MQ=255)
                     ccTTATCCGGCATTTAATAAATTACAGCCACTCTAc                >  1:77131/1‑36 (MQ=255)
                     ccTTATCCGGCATTTAATAAATTACAGCCACTCTAc                >  1:7451789/1‑36 (MQ=255)
                       ttATCCGGCATTTAATAAATTACAGCCACTCTAccc              >  1:4631942/1‑36 (MQ=255)
                          cccGGAATTAAAAAAATCACAGCCACTCTACCCGCa           <  1:3755998/35‑1 (MQ=16)
                           ccGGCATTTAATAAATTACAGCCACTCTACCCGCaa          >  1:8457670/1‑36 (MQ=255)
                             gacATTAAAAAAACCACAGCCACTCTACCCGCAACg        <  1:1482577/34‑1 (MQ=12)
                              gCATTTAATAAATTACAGCCACTCTACCCGCAACGa       >  1:6757819/1‑36 (MQ=255)
                               cATTTAATAAATTACAGCCACTCTACCCGCAACGAc      >  1:5761884/1‑36 (MQ=255)
                               cATTAAAAAAATCCCAGCCACTCTACCCGCAACGAc      <  1:8349450/36‑1 (MQ=16)
                               aaTTAAAAAAATACCAGCCACTCTACCCGCAACGAc      <  1:2535132/35‑1 (MQ=16)
                                 tttAATAAATTACAGCCACTCTACCCGCAACGACAt    >  1:2448036/1‑36 (MQ=255)
                                 tttAATAAATTACAGCCACTCTACCCGCAACGACAt    >  1:1370865/1‑36 (MQ=255)
                                 tttAATAAATTACAGCCACTCTACCCGCAACGACAt    >  1:3328692/1‑36 (MQ=255)
                                  ttAATAAATTACAGCCACTCTACCCTCAACGACAtc   >  1:2471031/1‑36 (MQ=255)
                                  ttAATAAATTACAGCCACTCTACCCGCAACGACAtc   >  1:214709/1‑36 (MQ=255)
                                  ttAATAAATTACAGCCACTCTACCCGCAACGACAtc   >  1:1972810/1‑36 (MQ=255)
                                  ttAATAAATTACAGCCACTCTACCCGCAACGACAtc   >  1:7532351/1‑36 (MQ=255)
                                  ttAATAAATTACAGCCACTCTACCCGCAACGACAtc   >  1:4509504/1‑36 (MQ=255)
                                   taataaATTACAGCCACTCTACCCGCAACGACAtct  >  1:7807623/1‑36 (MQ=255)
                                   |                                   
ATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCT  >  REL606/2594473‑2594543

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG