Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,370,027 T→A 100% 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% 10.8 71350rpsM30S ribosomal protein S13
Reads supporting (aligned to +/- strand):  new base (37/28):  ref base (1/5):  total (38/33)
Fisher's exact test strand distribution p-value = 9.00e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.84e-01

CCTTTGCCATTATTCAATCACCCCGATTATTTCTTGATCGGTTTGCGCGGACCCTTACGGGTACGT  >  REL606/3369994‑3370059
                                 |                                
ccTTTGCCATTATTCAATCACCCCGATTATTTCATg                                >  1:7895408/1‑36 (MQ=255)
ccTTTGCCATTATTCAATCACCCCGATTATTTCATg                                <  1:2497544/36‑1 (MQ=255)
ccTTTGCCATTATTCAATCACCCCGATTATTTCATg                                <  1:947778/36‑1 (MQ=255)
 cTTTGCCATTATTCAATCACCCCGATTATTTCATGa                               >  1:7223333/1‑36 (MQ=255)
 cTTTGCCATTATTCAATCACCCCGATTATTTCATGa                               >  1:1653450/1‑36 (MQ=255)
 cTTTGACATTATTCAATCACCCCGATTATTTCATGa                               >  1:4641127/1‑36 (MQ=37)
  tttGCCATTATTCAATCACCCCGATTATTTCATGAt                              <  1:5608883/36‑1 (MQ=255)
  tttGCCATTATTCAATCACCCCGATTATTTCATGAt                              <  1:3806021/36‑1 (MQ=255)
  tttGCCATTATTCAATCACCCCGATTATTTCATGAt                              >  1:8207582/1‑36 (MQ=255)
   ttGCCATTATTCAATCACCCCGATTATTTCATGATc                             <  1:6834100/36‑1 (MQ=255)
   gtGCCATTATTCAATCACCCCGATTATTTCATGATc                             >  1:5405339/2‑36 (MQ=255)
    tGCCATTATTCAATCACCCCGATTATTTCATGATCg                            >  1:1126312/1‑36 (MQ=255)
    tGCCATTATTCAATCACCCCGATTATTTCATGATCg                            >  1:7049001/1‑36 (MQ=255)
     gCCATTATTCAATCACCCCGATTATTTCATGATCgg                           >  1:3882779/1‑36 (MQ=255)
      ccATTATTCAATCACCCCGATTATTTCTTGATCGGt                          <  1:578939/36‑1 (MQ=255)
      ccATTATTCAATCACCCCGATTATTTCATGATCGGt                          <  1:6301894/36‑1 (MQ=255)
      ccATTATTCAATCACCCCGATTATTTCATGATCGGt                          >  1:6975972/1‑36 (MQ=255)
      ccATTATTCAATCACCCCGATTATTTCATGATCGGt                          >  1:313800/1‑36 (MQ=255)
      ccATTATTCAATCACCCCGATTATTTCATGATCGGt                          <  1:6911908/36‑1 (MQ=255)
      ccATTATTCAATCACCCCGATTATTTCATGATCGGt                          <  1:2436526/36‑1 (MQ=255)
       cATTATTCAATCACCCCGATTATTTCATGATCGGtt                         <  1:1755133/36‑1 (MQ=255)
        attattCAATCACCCCGATTATTTCATGATCGGttt                        <  1:4328954/36‑1 (MQ=255)
        attattCAATCACCCCGATTATTTCATGATCGGttt                        <  1:6133549/36‑1 (MQ=255)
         ttattCAATCACCCCGATTATTTCATGATCGGtttt                       >  1:3318985/1‑35 (MQ=255)
         ttattCAATCACCCCGATTATTTCATGATCGGTTTg                       >  1:5408331/1‑36 (MQ=255)
         ttattCAATCACCCCGATTATTTCATGATCGGTTTg                       >  1:8122604/1‑36 (MQ=255)
          tattCAATCACCCCGATTATTTCATGATCGGttttc                      >  1:824676/1‑34 (MQ=38)
          tattCAATCACCCCGATTATTTCATGATCGGttttc                      >  1:3189987/1‑34 (MQ=38)
          tattCAATCACCCCGATTATTTCATGATCGGttttc                      >  1:6137981/1‑34 (MQ=38)
          tattCAATCACCCCGATTATTTCATGATCGGTTTgc                      >  1:7131411/1‑36 (MQ=255)
           attCAATCACCCCGATTATTTCATGATCGGTTTgcg                     >  1:6123303/1‑36 (MQ=255)
           attCAATCACCCCGATTATTTCATGATCGGTTTgcg                     <  1:3679511/36‑1 (MQ=255)
           attCAATCACCCCGATTATTTCATGATCGGTTTgcg                     >  1:1327000/1‑36 (MQ=255)
            ttCAATCACCCCGATTATTTCTTGATCGGTTTgcgc                    >  1:67752/1‑36 (MQ=255)
             tCAATCACCCCGATTATTTCATGATCGGTTTGcgcg                   <  1:6875564/36‑1 (MQ=255)
              cAATCACCCCGATTATTTCATGATCGGTTTGCGCgg                  <  1:4197345/36‑1 (MQ=255)
               aaTCACCCCGATTATTTCATGATCGGTTTGCGCGGa                 >  1:2959158/1‑36 (MQ=255)
                aTCACCCCGATTATTTCATGATCGGTTTTCGCGGAc                >  1:6907458/1‑36 (MQ=37)
                aTCACCCCGATTATTTCATGATCGGTTTGCGCGGAc                <  1:5978188/36‑1 (MQ=255)
                 tCACCCCGATTATTTCATGATCGGTTTGCGCGGAcc               <  1:7597345/36‑1 (MQ=255)
                 tCACCCCGATTATTTCATGATCGGTTTGCGCGGAcc               <  1:4653213/36‑1 (MQ=255)
                 tCACCCCGATTATTTCATGATCGGTTTGCGCGGAcc               <  1:6159080/36‑1 (MQ=255)
                   aCCCCGATTATTTCATGATCGGTTTGCGCGGACCCt             >  1:4152755/1‑36 (MQ=255)
                    ccccGATTATTTCTTGATCGGTTTGCGCGGACCCtt            <  1:1416639/36‑1 (MQ=255)
                    ccccGATTATTTCTTGATCGGTTTGCGCGGACCCtt            <  1:839314/36‑1 (MQ=255)
                    ccccGATTATTTCATGATCGGTTTGCTCGGACCCtt            >  1:3104800/1‑36 (MQ=38)
                    ccccGATTATTTCATGATCGGTTTGCGCGGACCCtt            >  1:7630011/1‑36 (MQ=255)
                    ccccGATTAGTTCATGATCGGTTTGCGCGGACCCtt            <  1:4441379/36‑1 (MQ=38)
                      ccGATTATTTCATGATCGGTTTGTGCGGACACTTat          >  1:2054878/1‑35 (MQ=21)
                      ccGATTATTTCATGATCGGTTTGCGCGGACCCTTAc          <  1:8261631/36‑1 (MQ=255)
                      ccGATTATTTCATGATCGGTTTGCGCGGACCCTTAc          <  1:1346046/36‑1 (MQ=255)
                      ccGATTATTTCATGATCGGTTTGCGCGGACCCTTAc          <  1:5624391/36‑1 (MQ=255)
                       cGATTATTTCATGATCGGTTTGCGCGGACCCTTACg         <  1:7726250/36‑1 (MQ=255)
                        gATTATTTCATGATCGGTTTTCGCGGACCCTTACgg        >  1:6569091/1‑36 (MQ=37)
                        gATTATTTCATGATCGGTTTGCGCGGACCCTTACgg        >  1:5552601/1‑36 (MQ=255)
                          ttattTCTTGATCGGTTTGCGCGGACCCTTACGGGt      <  1:5434058/36‑1 (MQ=255)
                            attTCTTGATCGGTTTGCGCGGACCCTTACGGGTAc    <  1:4432492/36‑1 (MQ=255)
                            attTCATGATCGGTTTGCGCGGACCCTTGCGGGTAc    >  1:5908162/1‑36 (MQ=38)
                            attTCATGATCGGTTTGCGCGGACCCTTACGGGTAc    >  1:1165084/1‑36 (MQ=255)
                            attTCATGATCGGTTTGCGCGGACCCTTACGGGTAc    <  1:4277817/36‑1 (MQ=255)
                            attTCATGATCGGTTTGCGCGGACCCTTACGGGTAc    >  1:500097/1‑36 (MQ=255)
                            attTCATGATCGGTTTGCGCGGACCCTTACGGGTAc    <  1:8444599/36‑1 (MQ=255)
                            attTCATGATCGGTTTGCGCGGACCCTTACGGGTAc    <  1:6928662/36‑1 (MQ=255)
                            attTCATGATCGGGTTGGGAGGATCATTACGGGTAc    >  1:6395048/1‑36 (MQ=12)
                             tttCATGATCGGTTTGCGCGGACCCTTACGGGTACg   >  1:7698585/1‑36 (MQ=255)
                             tttCATGATCGGTTTGCGCGGACCCTTACGGGTACg   >  1:7147437/1‑36 (MQ=255)
                              ttCATGATCGGTTTGCGCGGACCCTTACGGGTACGt  <  1:1066401/36‑1 (MQ=255)
                              ttCATGATCGGTTTGCGCGGACCCTTACGGGTACGt  <  1:3287127/36‑1 (MQ=255)
                              ttCATGATCGGTTTGCGCGGACCCTTACGGGTACGt  >  1:3591943/1‑36 (MQ=255)
                              ttCATGATCGGTTTGCGCGGACCCTTACGGGTACGt  >  1:7538746/1‑36 (MQ=255)
                              ttCATGATCGGTTTGCGCGGACCCTTACGGGTACGt  >  1:6605189/1‑36 (MQ=255)
                                 |                                
CCTTTGCCATTATTCAATCACCCCGATTATTTCTTGATCGGTTTGCGCGGACCCTTACGGGTACGT  >  REL606/3369994‑3370059

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG