Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,370,027 T→A 59.4% K117M (AAG→ATG)  rpsM ← 30S ribosomal protein S13

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,370,0270T→A59.4% 72.0 ‑3.3 ‑1.759350rpsM30S ribosomal protein S13
Reads supporting (aligned to +/- strand):  new base (12/23):  ref base (15/9):  total (27/32)
Fisher's exact test strand distribution p-value = 3.84e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.02e-02

TGCCTTTGCCATTATTCAATCACCCCGATTATTTCTTGATCGGTTTGCGCGGACCCTTACGGGTACGTGC  >  REL606/3369992‑3370061
                                   |                                  
tGCCTTTGCCATTATTCAATCACCCCGATTATTTCt                                    <  1:6943290/36‑1 (MQ=255)
  ccTTTGCCATTATTCAATCACCCCGATTATTTCTTg                                  >  1:6444648/1‑36 (MQ=255)
     ttGCCATTATTCAATCACCCCGATTATTTCTTGATc                               <  1:6399127/36‑1 (MQ=255)
     ttGCCATTATTCAATCACCCCGATTATTTCATGATc                               <  1:5835610/36‑1 (MQ=255)
     ttGCCATTATTCAATCACCCCGATTATTTCATGATc                               >  1:7463842/1‑36 (MQ=255)
     ttGACATTATTCAATCACCCCAATTATGTCTTGATc                               >  1:3700368/1‑36 (MQ=25)
      tGCCATTATTCAATCACCCCGATTATTTCTTGATCg                              <  1:6833518/36‑1 (MQ=255)
      tGCCATTATTCAATCACCCCGATTATTTCATGATCg                              <  1:6186491/36‑1 (MQ=255)
      tGCCATTATTCAATCACCCCGATTATTTCATGATCg                              >  1:5110204/1‑36 (MQ=255)
      tGCCATTATTCAATCACCCCGATTATTTCATGATCg                              <  1:3623832/36‑1 (MQ=255)
       gCCATTATTCAATCACCCCGATTATTTCTTGATCga                             >  1:7666160/1‑35 (MQ=255)
       gCCATTATTCAATCACCCCGATTATTTCATGATCgg                             >  1:3442310/1‑36 (MQ=255)
       gCCATTATTCAATCACCCCGATTATTTCATGATCgg                             <  1:2357775/36‑1 (MQ=255)
       gCCATTATTCAATCACCCCGATTATTTCATGATCgg                             <  1:7853014/36‑1 (MQ=255)
          attattCAATCACCCCGATTATTTCTTGATCGGttt                          <  1:474646/36‑1 (MQ=255)
          attattCAATCACCCCGATTATTTCTTGATCGGttt                          >  1:5314353/1‑36 (MQ=255)
          attattCAATCACCCCGATTATTTCTTGATCGGttt                          >  1:5568911/1‑36 (MQ=255)
          attattCAATCACCCCGATTATTTCTTGATCGGttt                          >  1:6859182/1‑36 (MQ=255)
          attattCAATCACCCCGATTATTTCATGATCGGttt                          >  1:4400269/1‑36 (MQ=255)
           ttattCAATCACCCCGATTATTTCATGATCGGTTTg                         <  1:2510371/36‑1 (MQ=255)
            tattCAATCACCCCGATTATTTCATGATCGGTTTgc                        >  1:7065325/1‑36 (MQ=255)
            tattCAATCACCCCGATTATTTCATGATCGGTTTgc                        <  1:840665/36‑1 (MQ=255)
             attCAATCACCCCGATTATTTCTTGATCGGTTTgcg                       >  1:6070432/1‑36 (MQ=255)
             attCAATCACCCCGATTATTTCATGATCGGTTTgcg                       >  1:7323290/1‑36 (MQ=255)
              ttCAATCACCCCGATTATTTCATGATCGGTTTgcgc                      <  1:1297444/36‑1 (MQ=255)
              ttCAATCACCCCGATTATTTCATGATCGGTTTgcgc                      <  1:1146150/36‑1 (MQ=255)
               tCAATCACCCCGATTATTTCATGATCGGTTTGcgcg                     <  1:5388850/36‑1 (MQ=255)
                cAATCACCCCGATTATTTCATGATCGGTTTGCGCgg                    >  1:5608888/1‑36 (MQ=255)
                 aaTCACCCCGATTATTTCTTGATCGGTTTGCGCGGa                   >  1:2539692/1‑36 (MQ=255)
                 aaTCACCCCGATTATTTCATGATCGGTTTGCGCGGa                   >  1:1209726/1‑36 (MQ=255)
                  aTCACCCCGATTATTTCATGATCGGTTTGCGCGGAc                  <  1:7665363/36‑1 (MQ=255)
                   tCACCCCGATTATTTCATGATCGGTTTGCGCGGAcc                 <  1:247466/36‑1 (MQ=255)
                    cACCCCGATTATTTCTTGATCGGTTTGCGCGGAccc                >  1:2289122/1‑36 (MQ=255)
                    cACCCCGATTATTTCATGATCGGTTTGCGCGGAccc                >  1:5614157/1‑36 (MQ=255)
                     aCCCCGATTATTTCTTGATCGGTTTGCGCGGACCCt               >  1:4586323/1‑36 (MQ=255)
                     aCCCCGATTATTTCATGATCGGTTTGCGCGGACCCt               >  1:7247048/1‑36 (MQ=255)
                     aCCCCGATTATTTCATGATCGGTTTGCGAGGACCCt               >  1:3418367/1‑36 (MQ=38)
                      ccccGATTATTTCATGATCGGTTTGCGCGGACCCtt              <  1:988219/36‑1 (MQ=255)
                      ccccGATTATTTCATGATCGGTTTGCGCGGACCCtt              <  1:4738751/36‑1 (MQ=255)
                       cccGATTATTTCATGATCGGTTTGCGCGGACCCTTa             <  1:1782724/36‑1 (MQ=255)
                       cccGATTATTTCATGATCGGTTTGCGCGGACCCTTa             <  1:337605/36‑1 (MQ=255)
                       cccGATTATTTCATGATCGGTTTGCGCGGACCCTTa             <  1:6504899/36‑1 (MQ=255)
                        ccGATTATTTCTTGATCGGTTTGCGCGGACCCTTAc            >  1:2516016/1‑36 (MQ=255)
                         cGATTATTTCTTGATCGGTTTGCGCGGACCCTTACg           >  1:1263215/1‑36 (MQ=255)
                          gATTATTTCATGATCGGTTTGCGCGGACCCTTACgg          <  1:1500932/36‑1 (MQ=255)
                            ttattTCTTGATCGGTTTGCGCGGACCCTTACGGGt        <  1:671199/36‑1 (MQ=255)
                             tattTCATGATCGGTTTGCGCGGACCCTTACGGGTa       <  1:5689861/36‑1 (MQ=255)
                             tattTCATGATCGGTTTGCGCGGACCCTTACGGGTa       >  1:408849/1‑36 (MQ=255)
                              attTCTTGATCGGTTTGCGCGGACCCTTACGGGTAc      <  1:1613029/36‑1 (MQ=255)
                              attTCTTGATCGGTTTGCGCGGACCCTTACGGGTAc      <  1:6154927/36‑1 (MQ=255)
                              attTCTTGATCGGTTTGCGCGAACCCTTACGGGTAc      >  1:4669167/1‑36 (MQ=255)
                              attTCATGATCGGTTTGCGCGGACCCTTACGGGTAc      <  1:4019990/36‑1 (MQ=255)
                              attTCATGATCGGTTTGCGCGGACCCTTACGGGTAc      <  1:6316641/36‑1 (MQ=255)
                               tttCATGATCGGTTTGCGCGGACCCTTACGGGTACg     <  1:7475126/36‑1 (MQ=255)
                                ttCATGATCGGTTTGCGCGGACCCTTACGGGTACGt    <  1:3340917/36‑1 (MQ=255)
                                 tCTTGATCGGTTTGCGCGGACCCTTACGGGTACGTg   <  1:1110869/36‑1 (MQ=255)
                                 tCTTGATCGGTTTGCGCGGACCCTTACGGGTACGTg   <  1:7000553/36‑1 (MQ=255)
                                  cTTGATCGGTTTGCGCGGACCCTTACGGGTACGTGc  <  1:2563932/36‑1 (MQ=255)
                                  cTTGATCGGTTTGCGCGGACCCTTACGGGTACGTGc  >  1:3562137/1‑36 (MQ=255)
                                  cTTGATCGGTTTGCGCGGACCATTGCGGGTACGTGc  >  1:6274588/1‑36 (MQ=37)
                                   |                                  
TGCCTTTGCCATTATTCAATCACCCCGATTATTTCTTGATCGGTTTGCGCGGACCCTTACGGGTACGTGC  >  REL606/3369992‑3370061

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG