Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 851,606 G→A 18.4% S3F (TCT→TTT)  ybiU ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606851,6060G→A18.4% 8.6 ‑5.5 0.0548ybiUhypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/9):  ref base (30/13):  total (32/22)
Fisher's exact test strand distribution p-value = 4.15e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCG  >  REL606/851571‑851639
                                   |                                 
tgtgATCGGCAGGCAATGTGTCGCTGGTAAAAGTAg                                   <  1:8057130/36‑1 (MQ=255)
tgtgATCGGCAGGCAATGTGTCGCTGGTAAAAGTAg                                   >  1:6692773/1‑36 (MQ=255)
 gtgATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGa                                  >  1:3705218/1‑36 (MQ=255)
 gtgATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGa                                  >  1:1088970/1‑36 (MQ=255)
  tgATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGaa                                 <  1:8045390/36‑1 (MQ=255)
  tgATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGaa                                 <  1:1989202/36‑1 (MQ=255)
  tgATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGaa                                 >  1:1799398/1‑36 (MQ=255)
  tgATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGaa                                 >  1:6608459/1‑36 (MQ=255)
    aTCGGCAGGCAATGTGTCGCTGGTTAAAGTAAAAGc                               >  1:1748946/1‑36 (MQ=37)
    aTCGGCAGGCAATGTGTCGCTGGTAAAAGTAAAAGc                               >  1:1335350/1‑36 (MQ=255)
     tCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGcc                              >  1:6771783/1‑36 (MQ=255)
     tCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGcc                              <  1:2208077/36‑1 (MQ=255)
     tCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGcc                              <  1:8095727/36‑1 (MQ=255)
     tCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGcc                              >  1:6856096/1‑36 (MQ=255)
      cGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCa                             <  1:8283760/36‑1 (MQ=255)
      cGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCa                             <  1:4890593/36‑1 (MQ=255)
      cGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCa                             <  1:7719135/36‑1 (MQ=255)
       ggcaggcaATGTGTCGCTGGTAAAAGTAGAAGcctt                            >  1:1109067/1‑34 (MQ=255)
       ggcaggcaATGTGTCGCTGGTAAAAGTAGAAGCCAt                            <  1:1584901/36‑1 (MQ=255)
        gcaggcaATGTGTCGCTGGTAAAAGTAGAAGCCATa                           >  1:6409898/1‑36 (MQ=255)
        gcaggcaATGTGTCGCTGGTAAAAGTAGAAGCCATa                           <  1:3777390/36‑1 (MQ=255)
        gcaggcaATGTGTCGCTGGTAAAAGTAGAAGCCATa                           >  1:3758004/1‑36 (MQ=255)
        gcaggcaATGTGTCGCTGGTAAAAGTAGAAGCCATa                           >  1:1866470/1‑36 (MQ=255)
           ggcaATGTGTCGCTGGTAAAAGTAGAAGCCATAacc                        >  1:7761912/1‑36 (MQ=255)
           gaaaaTGTGTCGCTGGTAAAAGTAGAAGCCATAacc                        <  1:6584183/33‑1 (MQ=255)
            gcaATGTGTCGCTGGTAAAAGTAGAAGCCATAacca                       >  1:2684301/1‑36 (MQ=255)
            gaaaTGTGTCGCTGGTAAAAGTAGAAGCCATAacca                       <  1:2832226/34‑1 (MQ=255)
              aaTGTGTCGCTGGTAAAAGTAGAAGCCATAaccacc                     >  1:510551/1‑36 (MQ=255)
               aTGTGTCGCTGGTAAAAGTAGAAGCCATAACCACct                    >  1:6208907/1‑36 (MQ=255)
                tgtgTCGCTGGTAAAAGTAGAAGCCATAACCACctc                   >  1:3904022/1‑36 (MQ=255)
                tgtgTCGCTGGTAAAAGTAGAAGCCATAACCACctc                   >  1:8122159/1‑36 (MQ=255)
                tgtgTCGCTGGTAAAAGTAGAAGCCATAACCACctc                   >  1:5375341/1‑36 (MQ=255)
                 gtgtCGCTGGAAAAAGAAAAAGCCATAACCACctct                  <  1:8406166/36‑1 (MQ=21)
                 gggtCGCTGGAAAAAGAAAAAGCCATAACCACctct                  <  1:179738/34‑1 (MQ=18)
                  tgtCGCTGGTAAAAGTAGAAGCCATAACCACCtctc                 >  1:5484290/1‑36 (MQ=255)
                  tgtCGCTGGTAAAAGTAGAAGCCATAACCACCtctc                 >  1:5284108/1‑36 (MQ=255)
                  tgtCGCTGGAAAAAGTAAAAGCCATAACCACCtctc                 <  1:6657888/36‑1 (MQ=37)
                  tgtCGCTGATAAAAGTAAAAGCCATAACCACCtctc                 <  1:6240034/36‑1 (MQ=37)
                   gtCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCa                >  1:1850212/1‑36 (MQ=255)
                    tCGCTGGAAAAAGAAAAAGCCATAACCACCTCTCAt               <  1:5965267/36‑1 (MQ=21)
                     cacTGGAAAAAGTAGAAGCCATAACCACCTCTCAtt              <  1:955928/34‑1 (MQ=37)
                      gCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTa             >  1:5209601/1‑36 (MQ=255)
                      gCTGGAAAAAGTAAAAGCCATAACCACCTCTCATTa             <  1:7445686/36‑1 (MQ=37)
                       cTGGAAAAAGAAAAAGCCATAACCACCTCTCATTaa            <  1:4247839/36‑1 (MQ=21)
                         ggaaaaaGAAAAAGCCATAACCACCTCTCATTAACt          <  1:122325/33‑1 (MQ=37)
                         ggTAAAAGTAGAAGCCATAACCACCTCTCATTAACt          >  1:7679145/1‑36 (MQ=255)
                         ggTAAAAGTAGAAGCCATAACCACCTCTCATTAACt          >  1:8001733/1‑36 (MQ=255)
                          gTAAAAGTAGAAGCCATAACCACCTCTCATTAACtt         >  1:7664829/1‑36 (MQ=255)
                           tAAAAGTAGAAGCCATAACCACCTCTCATTAACttt        >  1:5111728/1‑36 (MQ=255)
                            aaaaGTAGAAGCCATAACCACCTCTCATTAACTTTc       >  1:6638636/1‑36 (MQ=255)
                            aaaaGAAAAAGCCATAACCACCTCTCATTAACTTTc       <  1:1772342/36‑1 (MQ=37)
                             aaaGTAGAAGCCATAACCACCTCTCATTAACTTTCa      >  1:5456964/1‑36 (MQ=255)
                              aattAGAAGCCATAACCACCTCTCATTAACTTTCAt     <  1:2363787/33‑1 (MQ=255)
                              aaGTAGAAGCCATAACCACCTCTCATTAACTTTCAt     >  1:4664208/1‑36 (MQ=255)
                               aGTAGAAGCCATAACCACCTCTCATTAACTTTCAtt    >  1:411095/1‑36 (MQ=255)
                                 tAGAAGCCATAACCACCTCTCATTAACTTTCATTCg  >  1:4681327/1‑36 (MQ=255)
                                   |                                 
TGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCG  >  REL606/851571‑851639

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG