Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,369,432 T→C D50G (GAC→GGC)  rpsD 30S ribosomal protein S4

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,369,4320T→C100.0% 206.4 63D50G (GAC→GGC) rpsD30S ribosomal protein S4
Reads supporting (aligned to +/- strand):  new base (44/18):  ref base (0/0):  total (44/19)

CTTTTTGCTTTTCACGCAACTGCACACCATAGTCAGACAGACGCGGTTTACGCGCACCGTGCTGG  >  REL606/3369400‑3369464
                                |                                
cTTTTTGCTTTTCACGCAACTGCACACCATAGCcag                               <  1:3734486/36‑1 (MQ=255)
cTTTTTGCTTTTCACGCAACTGCACACCATAGCcag                               >  1:2827041/1‑36 (MQ=255)
cTTTTTGCTTTTCACCCAACTGCACACCATAGCcag                               <  1:6554935/36‑1 (MQ=25)
  ttttGCTTTTCACGCAACTGCACACCATAGCcagtc                             >  1:1855894/1‑34 (MQ=25)
  ttttGCTTTTCACGCAACTGCACACCATAGCcagac                             >  1:3167795/1‑36 (MQ=255)
   tttGCTTTTCACGCAACTGCACACCATAGCcagaca                            >  1:1450586/1‑36 (MQ=255)
    ttGCTTTTCACGCAACTGCACACCATAGCcagacag                           >  1:7415941/1‑36 (MQ=255)
    ttGCTTTTCACGCAACTGCACACCATAGCcagacag                           >  1:369809/1‑36 (MQ=255)
    ttGCTTTTCACGCAACTGCACACCATAGCcagacag                           >  1:3355888/1‑36 (MQ=255)
     tGCTTTTCACGCAACTGCACACCATAGCcagacaga                          <  1:5926815/36‑1 (MQ=255)
      gCTTTTCACGCAACTGCACACCATAGCCagacagac                         >  1:494902/1‑36 (MQ=255)
      gCTTTTCACGCAACTGCACACCATAGCCagacagac                         <  1:5834242/36‑1 (MQ=255)
      gCTTTTCACGCAACTGCACACCATAGCCagacagac                         <  1:7439414/36‑1 (MQ=255)
       cTTTTCACGCAACTGCACACCATAGCCAGACAGAcg                        >  1:5515402/1‑36 (MQ=255)
       cTTTTCACGCAACTGCACACCATAGCCAGACAGAcg                        >  1:2936609/1‑36 (MQ=255)
        ttttCACGCAACTGCACACCATACCCACACACAcgc                       >  1:2886382/1‑36 (MQ=18)
         tttCACGCAACTGCACACCATAGCCAGACAGAcgcg                      >  1:1555998/1‑36 (MQ=255)
         tttCACGCAACTGCACACCATAGCCAGACAGAcgcg                      >  1:2133170/1‑36 (MQ=255)
          ttCACGCAACTGCACACCATAGCCAGACAGACGCgg                     >  1:4876034/1‑36 (MQ=255)
           tCACGCAACTGCACACCATAGCCAGACAGAcgcgct                    >  1:1899743/1‑34 (MQ=25)
           tCACGCAACTGCACACCATAGCCAGACAGACGCGGt                    >  1:4807284/1‑36 (MQ=255)
            cACGCAACTGCACACCATAGCCAGACAGACGAGGtt                   >  1:6766915/1‑36 (MQ=25)
             aCGCAACTGCACACCATAGCCAGACAGACGCGGttt                  <  1:5009754/36‑1 (MQ=255)
             aCGCAACTGCACACCATAGCCAGACAGACGCGGttt                  <  1:5529244/36‑1 (MQ=255)
             aCGCAACTGCACACCATAGCCAGACAGACGCGGttt                  >  1:3360607/1‑36 (MQ=255)
             aCGCAACTGCACACCATAGCCAGACAGACGCGGttt                  >  1:858144/1‑36 (MQ=255)
              cGCAACTGCACACCATAGCCAGACAGACGCGGTTTa                 >  1:4651140/1‑36 (MQ=255)
               gCAACTGCACACCATAGCCAGACAGACGCGGTTTAc                >  1:1970939/1‑36 (MQ=255)
               gCAACTGCACACCATAGCCAGACAGACGCGGTTTAc                >  1:1979882/1‑36 (MQ=255)
               aaaaCTACACACCATAGCCAGACAGACGCGGTTTAc                <  1:1141489/34‑1 (MQ=25)
                cAACTGCACACCATAGCCAGACAGACGCGGTTTAcg               >  1:1604366/1‑36 (MQ=255)
                cAACTGCACACCATAGCCAGACAGACGCGGTTTAcg               <  1:2415627/36‑1 (MQ=255)
                 aaCTGCACACCATAGCCAGACAGACGCGGTTTAcgc              >  1:5501233/1‑36 (MQ=255)
                 aaCTGCACACCATAGCCAGACAGACGCGGTTTAcgc              >  1:2876524/1‑36 (MQ=255)
                 aaCTGCACACCATAGCCAGACAGACGCGGTTTAcgc              <  1:4563514/36‑1 (MQ=255)
                  aCTGCACAGCATAGCCGGACAGACGCGGTTTAcgcg             <  1:2549300/36‑1 (MQ=21)
                  aCTGCACACCATTGCCAGTCAGACGCGGTTTAcgcg             >  1:4363932/1‑36 (MQ=21)
                  aCTGCACACCATAGCCAGACAGACGCGGTTTAcgcg             >  1:374614/1‑36 (MQ=255)
                  aCTGCACACCATAGCCAGACAGACGCGGTTTAcgcg             >  1:5086901/1‑36 (MQ=255)
                   cTGCACACCATAGCCAGACAGACGCGGTTTACgcgc            >  1:5469750/1‑36 (MQ=255)
                   cTGCACACCATAGCCAGACAGACGCGGTTTACgcgc            >  1:5981601/1‑36 (MQ=255)
                   cTGCACACCATAGCCAGACAGACGCGGTTTACgcgc            <  1:2973007/36‑1 (MQ=255)
                   cTGCACACCATAGCCAGACAGACGCGGTTTACgcgc            >  1:7325542/1‑36 (MQ=255)
                   cTGCACACCATAGCCAGACAGACGCGGTTTACgcgc            >  1:1668126/1‑36 (MQ=255)
                   cTGCACACCATAGCCAGACAGACGCGGTTTACacgc            >  1:4356502/1‑36 (MQ=25)
                    tGCACACCATAGCCAGACAGACGCGGTTTTCGCGCa           >  1:1120981/1‑36 (MQ=21)
                    tGCACACCATAGCCAGACAGACGCGGTTTACGCGCa           >  1:7378104/1‑36 (MQ=255)
                    tGCACACCATAGCCAGACAGACGCGGTTTACGCGCa           >  1:3825907/1‑36 (MQ=255)
                    tGCACACCATAGCCAGACAGACGCGGTTTACACGCa           >  1:2239932/1‑36 (MQ=25)
                     gCACACCATAGCCAGACAGACGCGGTTTACGCGCAc          >  1:122656/1‑36 (MQ=255)
                     gCACACCATAGCCAGACAGACGCGGTTTACGCGCAc          <  1:2633350/36‑1 (MQ=255)
                     gCACACCATAGCCAGACAGACGCGGTTTACGCGCAc          <  1:3744109/36‑1 (MQ=255)
                     gCACACCATAGCCAGACAGACGCGGTTTACGCGCAc          <  1:5197057/36‑1 (MQ=255)
                      cacaCCATAGCCAGACAGACGCGGTTTATGCGTAca         >  1:6442356/1‑35 (MQ=21)
                       acacCATAGCCAGACAGACGCGGTTTACGCGCACCg        <  1:7265844/36‑1 (MQ=255)
                       acacCAGAGCCAGACAGACGCGGTTTACGCGCACCg        <  1:42763/36‑1 (MQ=21)
                        cacCATAGCCAGACAGACGCGGTTTACGCGCACCGt       <  1:5674173/36‑1 (MQ=255)
                         acCATAGCCAGACAGACGCGGTTTACGCGCACCGTg      >  1:6784021/1‑36 (MQ=255)
                          ccATAGCCAGACAGACGCGGTTTACGCGCACCGTGc     >  1:3737860/1‑36 (MQ=255)
                          ccATAGCCAGACAGACGCGGTTTACGCGCACCGTGc     >  1:2280308/1‑36 (MQ=255)
                          ccATAGCCAGACAGACGCGGTTTACGCGCACCGTGc     >  1:1501667/1‑36 (MQ=255)
                           cATAGCCAGACAGACGCGGTTTACGCGCACCGTGCt    >  1:5166114/1‑36 (MQ=255)
                             tAGACAGACAGACGCGGTTTACGCGCACCGTGCTgg  <  1:6205993/36‑1 (MQ=255)
                                |                                
CTTTTTGCTTTTCACGCAACTGCACACCATAGTCAGACAGACGCGGTTTACGCGCACCGTGCTGG  >  REL606/3369400‑3369464

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG