Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,228,283 G→T P96Q (CCG→CAG)  yraL ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,228,2830GT100.0% 180.1 / NA 45P96Q (CCG→CAG) yraLpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (20/25);  total (20/25)

GTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCGGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGTTTCAGCTTTT  >  REL606/3228212‑3228354
                                                                       |                                                                       
ntAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGt                                                                     >  1:1038790/2‑76 (MQ=255)
 tAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGtt                                                                    <  2:912509/76‑1 (MQ=255)
  aGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTc                                                                   <  1:1379647/76‑1 (MQ=255)
     ggCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCgg                                                                >  2:1826636/1‑76 (MQ=255)
        accacGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCAt                                                             <  2:1921720/76‑1 (MQ=255)
          cacGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCg                                                           <  1:903882/76‑1 (MQ=255)
           acGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCgg                                                          <  1:1546084/76‑1 (MQ=255)
            cgcgGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTACGGCATCGGa                                                         >  2:1393628/1‑76 (MQ=255)
             gcgGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGaa                                                        <  2:421524/76‑1 (MQ=255)
             gcgGATCCCCGCTTCACGGCAGGTACGCACCAGATGGGAGCCAGGATCGTTAATTAGCTGCGGTCCGGCATCGGaa                                                        >  2:1198975/1‑76 (MQ=255)
              cgGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGaan                                                       <  1:1915768/76‑2 (MQ=255)
              cgGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGaaa                                                       <  2:1508991/76‑1 (MQ=255)
               ggATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAAc                                                      <  2:278679/76‑1 (MQ=255)
               ggATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAAc                                                      <  1:1634273/76‑1 (MQ=255)
               ggATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAAc                                                      <  1:998925/76‑1 (MQ=255)
                  tCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAg                                                   >  2:442264/1‑76 (MQ=255)
                     ccGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAgcgc                                                <  1:424718/76‑1 (MQ=255)
                     ccGCTTCACGGCAGGGACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAgcgc                                                >  2:790645/1‑76 (MQ=255)
                      cGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCa                                               <  1:956017/76‑1 (MQ=255)
                      cGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCa                                               >  2:1268078/1‑76 (MQ=255)
                        ctccACGGCAGGTACGCGCCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAAt                                             <  2:1530013/73‑1 (MQ=255)
                        cTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAAt                                             >  2:1403079/1‑76 (MQ=255)
                         ttCACGGCAGGTCCGCCCCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCaaag                                            <  2:1115543/76‑3 (MQ=255)
                             cGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGtttt                                        >  2:1578419/1‑76 (MQ=255)
                             cGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGtttt                                        <  1:1215832/76‑1 (MQ=255)
                               gCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGc                                      <  1:1053912/76‑1 (MQ=255)
                                 aGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGccc                                    >  1:918806/1‑76 (MQ=255)
                                  ggTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGcccc                                   <  2:519130/76‑1 (MQ=255)
                                      cGCACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCtt                               >  2:786487/1‑76 (MQ=255)
                                        cACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGc                             <  2:150524/76‑1 (MQ=255)
                                        cACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGc                             >  1:821242/1‑76 (MQ=255)
                                        cACCAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGc                             >  2:167239/1‑76 (MQ=255)
                                           cAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGc                          >  1:500619/1‑76 (MQ=255)
                                           cAGATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGc                          >  1:347385/1‑76 (MQ=255)
                                             gATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCtt                        >  1:1883621/1‑76 (MQ=255)
                                             gATGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCtt                        <  2:517233/76‑1 (MQ=255)
                                              aTGGTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTc                       <  2:84593/76‑1 (MQ=255)
                                                 gTAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGcc                    >  2:1538571/1‑76 (MQ=255)
                                                  tAGCCAGGATCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCca                   >  1:1455590/1‑76 (MQ=255)
                                                   aGCCAGGATCGTTACTTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCcag                  >  1:1061617/1‑76 (MQ=255)
                                                          aTCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGtt           <  2:535124/76‑1 (MQ=255)
                                                          aTCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGtt           >  2:1296103/1‑76 (MQ=255)
                                                           tCGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGttt          <  2:1279248/76‑1 (MQ=255)
                                                            cGTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGTTTc         >  2:182608/1‑76 (MQ=255)
                                                             gTTAATTAGCTGCGTTCCGGCATCGGAAACCGGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGTTTCa        <  2:1560142/76‑1 (MQ=255)
                                                             gTTAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGTTTCa        >  2:985538/1‑76 (MQ=255)
                                                               tAATTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGTTTCAGc      >  1:482438/1‑76 (MQ=255)
                                                                 aTTAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGTTTCAGCtt    >  2:852164/1‑76 (MQ=255)
                                                                   tAGCTGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGTTTCAGCtttt  <  2:1038790/76‑1 (MQ=255)
                                                                       |                                                                       
GTAGCGGCACCACGCGGATCCCCGCTTCACGGCAGGTACGCACCAGATGGTAGCCAGGATCGTTAATTAGCGGCGTTCCGGCATCGGAAACCAGCGCAATGTTTTGCCCCTCTTGCAGCTTCGCCAGCAGCGTTTCAGCTTTT  >  REL606/3228212‑3228354

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: