Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,125,281 G→T 13.8% G80C (GGT→TGT)  yqiJ → predicted inner membrane protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,125,2810G→T13.8% 3.1 ‑3.6 0.057238yqiJpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  new base (8/1):  ref base (22/26):  total (30/27)
Fisher's exact test strand distribution p-value = 2.71e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CTCTCGTCGTTCTCTGCTTACTGGCTGGTTTTTTTGGTCTTATTGGCATCCTGTTGCAGCATGCCTGCATC  >  REL606/3125246‑3125316
                                   |                                   
ctctCGTCGTTCTCTGCTTACTGGCTGGTTTTTTTg                                     <  1:1915679/36‑1 (MQ=255)
   tcgtcgTTCTCTGCTTACTGGCTGGTTTTTTTGGTc                                  <  1:2124735/36‑1 (MQ=255)
    cgtcgtTCTCTGCTTCCTGGCTGGTTTTTTTGGTCt                                 <  1:5920552/36‑1 (MQ=255)
     gtcgtTCTCTGCTTACTGGCTGGTTTTTTTTGTTtt                                >  1:2825017/1‑36 (MQ=25)
     gtcgtTCTCTGCTTACTGGCTGGTTTTTTTTGTCtt                                >  1:4169639/1‑36 (MQ=255)
     gtcgtTCTCTGCTTACTGGCTGGTTTTTTTGGTCtt                                <  1:842607/36‑1 (MQ=255)
     gtcgtTCTCTGCTTACTGGCTGGTTTTTTTGGTCtt                                <  1:1762577/36‑1 (MQ=255)
     gtcgtTCTCTGCTTACTGGCTGGTTTTTTTGGTCtt                                <  1:2302879/36‑1 (MQ=255)
     gacgtTCTCTGCTTACTGGCTGGCTTTTTTTGTTtt                                >  1:5388147/3‑36 (MQ=21)
      tcgtTCCCTGCTTACTGGCTGGTTTTTTTGGTCTTa                               <  1:5407041/36‑1 (MQ=255)
       cgtTCTCTGCTTACTGGCTGGTTTTTTTGGTCTTAt                              <  1:5711313/36‑1 (MQ=255)
        gtTCTCTGCTTACTGGCTGGTTTTTTTTGTCTTAtt                             >  1:2106672/1‑36 (MQ=255)
        gtTCTCTGCTTACTGGCTGGTTTTTTTGGTCTTAtt                             <  1:1392784/36‑1 (MQ=255)
        gtTCTCTGCTTACTGGCTGGGTTTTTTTGTCTTTtt                             >  1:2272826/1‑36 (MQ=25)
          tctcTGCTTACTGGCTGGTTTTTTTTGTTTTAtttt                           >  1:4892071/1‑34 (MQ=21)
           ctctGCTTACTGGCTGGTTTTTTTGGTCTTATTGGc                          <  1:2705604/36‑1 (MQ=255)
            tctGCTTACTGGCTGGTTTTTTTTGTCTTATTTGcc                         >  1:2090690/1‑35 (MQ=25)
              tGCTTACTGGCTGGTTTTTTTGGTCTTATTTGCATc                       >  1:2273981/1‑36 (MQ=255)
              tGCTTACTGGCTGGTTTTTTTGGTCTTATTGGCATc                       <  1:4935921/36‑1 (MQ=255)
               gCTTACTGGCTGGTTTTTTTGGTCTTATTTTCATcc                      >  1:3018713/1‑36 (MQ=37)
                cTTACTGGCTGGTTTTTTTTGTCTTATTTGCATCCt                     >  1:316425/1‑36 (MQ=37)
                cTTACTGGCTGGTTTTTTTGGTCTTATTGGCATCCt                     <  1:4642493/36‑1 (MQ=255)
                  tACTGGCTGGTTTTTTTGGTCTTATTGGCATCCTGt                   <  1:5433314/36‑1 (MQ=255)
                  tACTGGCTGGTTTTTTTGGTCTTATTGGCATCCTGt                   <  1:2377081/36‑1 (MQ=255)
                    ctggctggTTTTTTTGGTCTTATTGGCATCCTGTTg                 <  1:251438/36‑1 (MQ=255)
                     tggctggTTTTTTTGGTCTTATTTGCATCCTGTTgc                >  1:5608752/1‑36 (MQ=255)
                      ggctggTTTTTTTGGTCTTATTGGCATCCTTTTtca               >  1:4224586/1‑36 (MQ=37)
                      ggctggTTTTTTTGGTCTTATTGGCATCCTGTTgca               >  1:4574149/1‑36 (MQ=255)
                      ggctggTTTTTTTGGTCTTATTGGCATCCTGTTgca               >  1:2052592/1‑36 (MQ=255)
                      ggctggTTTTTTTGGTCTTATTGGCATCCTGTTgca               <  1:3847431/36‑1 (MQ=255)
                      ggcgggTTTTTTTGGTCTTATTGGCATCCTTTTtca               >  1:2819624/1‑36 (MQ=25)
                       gctggTTTTTTTTGTCTTCTTGGCATCCTGTTgcag              <  1:2286069/36‑1 (MQ=37)
                       gctggTTTTTTTGGTCTTATTGGCATCCTGTTgcag              <  1:752863/36‑1 (MQ=255)
                       gctggTTTTTTTGGTCTTATTGGCATCCTGTTgcag              >  1:4218677/1‑36 (MQ=255)
                       gctggTTTTTTTGGTCTTATTGGCATCCTGTTgcag              >  1:3543307/1‑36 (MQ=255)
                       gctggTTTTTTTGGTCTTATTGGCATCCTGTTgcag              <  1:5301969/36‑1 (MQ=255)
                       gctggTTTTTTTGGTCTTATTGGCATCCTGTTgcag              >  1:3432524/1‑36 (MQ=255)
                        ctggTTTTTTTGGTCTTATTGGCATCCTGTTgctgc             >  1:946894/1‑36 (MQ=255)
                        ctggTTTTTTTGGTCTTATTGGCATCCTGTTgcagc             >  1:3412093/1‑36 (MQ=255)
                        ctggTTTTTTTGGTCTTATTGGCATCCTGTTgcagc             <  1:271590/36‑1 (MQ=255)
                         tggTTTTTTTGGTCTTATTGGCATCTTGTTgcagca            >  1:513390/1‑36 (MQ=255)
                         tggTTTTTTTGGTCTTATTGGCATCCTGTTgcagca            >  1:3453308/1‑36 (MQ=255)
                         tggTTTTTTTGGTCTTATTGGCATCCTGTTgcagca            >  1:3133230/1‑36 (MQ=255)
                          ggTTTTTTTGGTCTTATTGGCATCCTGTTGCAGCAt           >  1:847225/1‑36 (MQ=255)
                              tttttGGTCTTATTGGCATCCTGTTGCAGCATGCCt       >  1:1095016/1‑36 (MQ=255)
                              tttttGGTCTTATTGGCATCCTGTTGCAGCATCCCt       <  1:1808462/36‑1 (MQ=255)
                               ttttGGTCTTATTGGCATCCTGTTGCAGCATGCCTg      >  1:4904161/1‑36 (MQ=255)
                               ttttGGTCTTATTGGCATCCTGTTGCAGCATGCCTg      >  1:1276704/1‑36 (MQ=255)
                               ttttGGTCTTATTGGCATCCTGTTGCAGCATGCCTg      <  1:1776067/36‑1 (MQ=255)
                                tttGGTCTTATTGGCATCCTGTTGCAGCATTCCTGc     <  1:5584418/36‑1 (MQ=255)
                                tttGGTCTTATTGGCATCCTGTTGCAGCATGCCTGc     <  1:1298/36‑1 (MQ=255)
                                tttGGTCTTATTGGCATCCTGTTGCAGCATGCCTGc     >  1:5810955/1‑36 (MQ=255)
                                tttGGTCTTATTGGCATCCTGTTGCAGCATGCCTGc     <  1:5041398/36‑1 (MQ=255)
                                tttGGTCTTATTGGCATCCTGTTGCAGCATGCCTGc     >  1:4314939/1‑36 (MQ=255)
                                tttGGTCTTATTGGCATCCTGTTGCAGCATGCCTGc     <  1:4404134/36‑1 (MQ=255)
                                 ttGGTCTTATTGGCATCCTGTTGCAGCATGCCTGca    <  1:3516881/36‑1 (MQ=255)
                                  tGGTCTTATTGGCATCCTGTTGCAGCATGCCTGcat   >  1:5524262/1‑36 (MQ=255)
                                  tGGTCTTATTGGCATCCTGTTGCAGCATGCCTGcat   <  1:5486934/36‑1 (MQ=255)
                                   ggTCTTATTGGCATCCTGTTGCAGCATGCCTGcatc  <  1:1880885/36‑1 (MQ=255)
                                   ggTCTTATTGGCATCCTGTTGCAGCATGCCTGcatc  <  1:2306130/36‑1 (MQ=255)
                                   ggTCTTATTGGCATCCTGTTGCAGCATGCCTGcatc  <  1:1492838/36‑1 (MQ=255)
                                   |                                   
CTCTCGTCGTTCTCTGCTTACTGGCTGGTTTTTTTGGTCTTATTGGCATCCTGTTGCAGCATGCCTGCATC  >  REL606/3125246‑3125316

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