Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 866,370 G→T 24.2% L196F (TTG→TTT yliE → conserved inner membrane protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606866,3700G→T24.2% 7.6 ‑4.3 0.032588yliEconserved inner membrane protein
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (12/12):  total (20/12)
Fisher's exact test strand distribution p-value = 1.35e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GAAATCTATTTTTATTATCCTGGCGGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGC  >  REL606/866335‑866405
                                   |                                   
gAAATCTATTTTTATTATCCTGGCGGTAATTTATTg                                     <  1:2995908/36‑1 (MQ=255)
  aaTCTATTTTTATTATCCTGGCGGTAATTTATTGaa                                   <  1:2262901/36‑1 (MQ=255)
       aTTTTTATTATCCTGGCGGTAATTTATTGAACGCAt                              <  1:2940661/36‑1 (MQ=255)
        tttttATTATCCTGGCGGTAATTTATTTAACGCATc                             >  1:867072/1‑36 (MQ=255)
        tttttATTATCCTGGCGGTAATTTATTGAACGCATc                             <  1:839820/36‑1 (MQ=255)
        tttttATTATCCTGGCGGTAATTTATTGAACGCATc                             >  1:4458160/1‑36 (MQ=255)
         ttttATTATCCTGGCGGTAATTTATTGAACGCATCa                            <  1:3072479/36‑1 (MQ=255)
             attatCCTGGCGGTAATTTATTTAACGCATCAGATa                        >  1:5894952/1‑36 (MQ=255)
             attatCCTGGCGGTAATTTATTTAACGCATCAGATa                        >  1:5392611/1‑36 (MQ=255)
             attatCCTGGCGGCAATTTATTTAACGCATCAGATa                        >  1:7171707/1‑36 (MQ=37)
              ttatCCTGGCGGTAATTTATTTAACGCATCTGATaa                       >  1:7367170/1‑36 (MQ=37)
               tatCCTGGCGGTAATTTATTGAACGCATCAGATaaa                      >  1:1107563/1‑36 (MQ=255)
                atCCTGGCGGTAATTTATTGAACGCATCAGATAAAc                     <  1:3597974/36‑1 (MQ=255)
                 tCCTGGCGGTAATTTTTTTAAAGCATCAGATaaaat                    >  1:2423047/1‑34 (MQ=18)
                 tCCTGGCGGTAATTTATTTAACGCATCAGATAAACt                    >  1:560692/1‑36 (MQ=255)
                 tCCTGGCGGTAATTTATTGAACGCATCAGATAAACt                    <  1:3819963/36‑1 (MQ=255)
                    tGGCGGTAATTTTTTTAACGCATTAGATAAAAtttt                 >  1:5536490/1‑36 (MQ=16)
                    tGGCGGTAATTTATTGAACGCATCAGATAAACtttt                 <  1:6386154/36‑1 (MQ=255)
                       cGGTAATTTATTGAACGCATCAGATAAACTTTTTTa              >  1:3869438/1‑36 (MQ=255)
                         gTAATTTATTGAACGCATCAGATAAACTTTTTTATc            >  1:7803274/1‑36 (MQ=255)
                         gTAATTTATTGAACGCATCAGATAAACTTTTTTATc            >  1:7189010/1‑36 (MQ=255)
                          tAATTTATTGAACGCATCAGATAAACTTTTTTATCt           <  1:2614626/36‑1 (MQ=255)
                          tAATTTATTGAACGCATCAGATAAACTTTTTTATCt           >  1:4616955/1‑36 (MQ=255)
                            aTTTATTGAACGCATCAGATAAACTTTTTTATCTGc         >  1:7168880/1‑36 (MQ=255)
                            aTTTATTGAACGCATCAGATAAACTTTTTTATCTGc         >  1:3008141/1‑36 (MQ=255)
                            aTTTATTGAACGCATCAGATAAACTTTTTTATCTGc         <  1:7891433/36‑1 (MQ=255)
                             tttATTGAACGCATCAGATAAACTTTTTTATCTGcc        >  1:5538828/1‑36 (MQ=255)
                               tATTGAACGCATCAGATAAACTTTTTTATCTGCCGt      >  1:2181395/1‑36 (MQ=255)
                                aTTGAACGCATCAGATAAACTTTTTTATCTGCCGtt     >  1:6553941/1‑36 (MQ=255)
                                 ttGAACGCATCAGATAAACTTTTTTATCTGCCGttt    <  1:2619385/36‑1 (MQ=255)
                                 ttGAACGCATCAGATAAACTTTTTTATCTGCCGttt    >  1:4210701/1‑36 (MQ=255)
                                 ttGAACGCATCAGATAAACTTTTTTATCTGCCGttt    <  1:5394980/36‑1 (MQ=255)
                                  tGAACGCATCAGATAAACTTTTTTATCTGCCGTTTg   <  1:6991765/36‑1 (MQ=255)
                                   gAACGCATCAGATAAACTTTTTTATCTGCCGTTTGc  <  1:680848/36‑1 (MQ=255)
                                   |                                   
GAAATCTATTTTTATTATCCTGGCGGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGC  >  REL606/866335‑866405

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG