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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,370,0270T→A100.0% 184.6 56350rpsM30S ribosomal protein S13
Reads supporting (aligned to +/- strand):  new base (22/34):  ref base (0/0):  total (22/34)

CCTTTGCCATTATTCAATCACCCCGATTATTTCTTGATCGGTTTGCGCGGACCCTTACGGGTACGT  >  REL606/3369994‑3370059
                                 |                                
cgtttGGCATTATTCAATCACCCCGATTATTTCATg                                <  1:1791763/34‑1 (MQ=25)
ccTTTGCCATTATTCAATCACCCCGATTATTTCATg                                <  1:2734448/36‑1 (MQ=255)
ccTTTGCCATTATTCAATCACCCCGATTATTTCATg                                <  2:2232094/36‑1 (MQ=255)
 cTTTGCCATTATTCAATCACCCCGATTATTTCATGa                               >  1:2991536/1‑36 (MQ=255)
 cTTTGCCATTATTCAATCACCCCGATTATTTCATGa                               >  2:1751892/1‑36 (MQ=255)
 cTTTGCCATTATTCAATCACCCCGATTATTTCATGa                               <  1:1721314/36‑1 (MQ=255)
    tGCCATTATTCAATCACCCCGATTATTTCATGATca                            <  1:2827482/36‑2 (MQ=255)
    tGCCATTATTCAATCACCCCGATTATTTCATGATCg                            >  1:3354466/1‑36 (MQ=255)
    tGCCATTATTCAATCACCCCGATTATTTCATGATCg                            <  2:3012736/36‑1 (MQ=255)
     gCCATTATTCAATCACCCCGATTATTTCATGATCgg                           >  1:2688297/1‑36 (MQ=255)
     gCCATTATTCAATCACCCCGAATATTTCATGATCgg                           <  2:2631381/36‑1 (MQ=37)
      ccATTATTCAATCACCCCGATTATTTCATGATCGGt                          >  2:2833560/1‑36 (MQ=255)
      ccATTATTCAATCACCCCGATTATTTCATGATCGGt                          <  1:16553/36‑1 (MQ=255)
      ccATTATTCAATCACCCCGATTATTTCATGATCGGt                          <  2:464607/36‑1 (MQ=255)
      ccATTATTCAATCACCCCGATTATTTCATGATCGGt                          >  1:3311939/1‑36 (MQ=255)
      ccAGTATTCAATCACCCCGATTATTTCATGATCGGt                          <  2:1614907/36‑1 (MQ=37)
       cATTATTCAATCACCCCGATTATTTCATGATCGGtt                         >  2:1907269/1‑36 (MQ=255)
       cATTATTCAATCACCCCGATTATTTCATGATCGGtt                         >  1:3112723/1‑36 (MQ=255)
        attattCAATCACCCCGATTATTTCATGATCGGttt                        >  1:3732391/1‑36 (MQ=255)
        attattCAATCACCCCGATTATTTCATGATCGGttt                        >  1:427871/1‑36 (MQ=255)
         ttattCAATCACCCCGATTATTTCATGATCGGTTTg                       <  1:1043146/36‑1 (MQ=255)
           attCAATCACCCCGATTATTTCATGATCGGTTTgcg                     <  2:175799/36‑1 (MQ=255)
           attCAATCACCCCGATTATTTCATGATCGGTTTgcg                     <  1:677851/36‑1 (MQ=255)
            ttCAATCACCCCGATTATTTCATGATCGGTTTgcgc                    <  1:2074755/36‑1 (MQ=255)
            ttCAATCACCCCGATTATTTCATGATCGGTTTgcgc                    >  2:3247990/1‑36 (MQ=255)
            ttCAATCACCCCGATTATTTCATGATCGGTTTgcgc                    <  1:698054/36‑1 (MQ=255)
            ttCAATCACCCCGATTATTTCATGATCGGTTTgcgc                    <  2:1170420/36‑1 (MQ=255)
              cAATCACCCCGATTATTTCATGATCGGTTTGCGCgg                  >  2:1684069/1‑36 (MQ=255)
              cAATCACCCCGATTATTTCATGATCGGTTTGCGCgg                  <  2:1543032/36‑1 (MQ=255)
              cAATCACCCCGATTATTTCATGATCGGTTTGCGCgg                  <  2:1294555/36‑1 (MQ=255)
              cAATCACCCCGATTATTTCATGATCGGTTTGCGCgg                  >  1:702598/1‑36 (MQ=255)
              cAATCACCCCGATTATTTCATGATCGGTTTGCGCgg                  <  1:1046515/36‑1 (MQ=255)
                aTCACCCCGATTATTTCATGATCGGTTTGCGCGGAc                <  2:2138193/36‑1 (MQ=255)
                  cACCCCGATTATTTCATGATCGGTTTGCGCGGAccc              <  2:2599077/36‑1 (MQ=255)
                  cACCCCGATTATTTCATGATCGGTTTGCGCGGAccc              >  1:2455065/1‑36 (MQ=255)
                   aCCCCGATTATTTCATGATCGGTTTGCGCGGACCCt             <  2:1649099/36‑1 (MQ=255)
                   aCCCCGATTATTTCATGATCGGTTTGCGCGGACCCt             <  2:1973516/36‑1 (MQ=255)
                    ccccGATTATTTCATGATCGGTTTGCGCGGACCCtt            >  1:1079620/1‑36 (MQ=255)
                    ccccGATTATTTCATGATCGGTTTGCGCGGACCCtt            >  1:2401287/1‑36 (MQ=255)
                      ccGATTATTTCATGATCGGTTTGCGCGGACCCTTAc          <  1:3388816/36‑1 (MQ=255)
                       cGATTATTTCATGATCGGTTTGCGCGGACCCTTACg         <  1:26093/36‑1 (MQ=255)
                       cGATTATTTCATGATCGGTTTGCGCGGACCCTTACg         <  2:740499/36‑1 (MQ=255)
                        gATTATTTCATGATCGGTTTGCGCGGACCCTTACgg        <  1:2319218/36‑1 (MQ=255)
                        gATTATTTCATGATCGGTTTGCGCGGACCCTTACgg        <  2:402388/36‑1 (MQ=255)
                        gATTATTTCATGATCGGTTTGCGCGGACCCTTACgg        >  2:3618481/1‑36 (MQ=255)
                         attattTCATGATCGGTTTGCGCGGACCCTTACggg       <  2:1770780/36‑1 (MQ=255)
                         attattTCATGATCGGTTTGCGCGGACCCTTACggg       >  2:458112/1‑36 (MQ=255)
                         attattTCATGATCGGTTTGCGCGGACCCTTACggg       >  2:1280472/1‑36 (MQ=255)
                           ttattCATGATCGGTTTGCGCGGACCCTTACGGGTa     <  1:1930949/33‑1 (MQ=38)
                           tattTCATGATCGGTTTGCGCGGACCCTTACGGGTa     >  2:946237/1‑36 (MQ=255)
                           tattTCATGATCGGTTTGCGCGGACCCTTACGGGTa     <  1:1106239/36‑1 (MQ=255)
                            attTCATGATCGGTTTGCGCGGACCCTTACGGGTAc    <  2:2955400/36‑1 (MQ=255)
                            attTCATGATCGGTTTGCGCGGACCCTTACGGGTAc    <  1:3408130/36‑1 (MQ=255)
                             tttCATGATCGGTTTGCGCGGACCCTTACGGGTACg   <  1:3593969/36‑1 (MQ=255)
                              ttCATGATCGGTTTGCGCGGACCCTTACGGGTACGt  >  2:2478776/1‑36 (MQ=255)
                              ttCATGATCGGTTTGCGCGGACCCTTACGGGTACGt  >  1:3389483/1‑36 (MQ=255)
                                 |                                
CCTTTGCCATTATTCAATCACCCCGATTATTTCTTGATCGGTTTGCGCGGACCCTTACGGGTACGT  >  REL606/3369994‑3370059

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG