Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,252,016 C→T G249S (GGT→AGT)  nusA ← transcription elongation factor NusA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,252,0160CT100.0% 164.0 / NA 42G249S (GGT→AGT) nusAtranscription elongation factor NusA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/24);  total (18/24)

TATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTG  >  REL606/3251944‑3252087
                                                                        |                                                                       
tatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTAc                                                                      <  1:1150173/76‑1 (MQ=255)
tatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTAc                                                                      >  1:1280889/1‑76 (MQ=255)
tatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTAc                                                                      >  1:1767097/1‑76 (MQ=255)
  tCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCg                                                                    <  1:1526992/76‑1 (MQ=255)
   cGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCgg                                                                   <  1:1040703/76‑1 (MQ=255)
   cGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCgg                                                                   >  1:589202/1‑76 (MQ=255)
     aTACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGAt                                                                 <  2:1796385/76‑1 (MQ=255)
      tACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATc                                                                >  1:618733/1‑76 (MQ=255)
          cTCGCCCCCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATa                                                            <  1:1281170/76‑1 (MQ=255)
          cTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATa                                                            <  1:23664/76‑1 (MQ=255)
           tCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATAc                                                           <  1:295848/76‑1 (MQ=255)
           tCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATAc                                                           >  2:39545/1‑76 (MQ=255)
           tCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATAc                                                           >  1:469222/1‑76 (MQ=255)
           tCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATAc                                                           <  1:1423490/76‑1 (MQ=255)
               cACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGttt                                                       <  1:328371/76‑1 (MQ=255)
               cACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGttt                                                       <  2:1785676/76‑1 (MQ=255)
               cACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGttt                                                       <  1:805678/76‑1 (MQ=255)
               cACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGttt                                                       >  2:1528949/1‑76 (MQ=255)
                aCCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTa                                                      <  1:646461/76‑1 (MQ=255)
                 cccAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTAt                                                     <  1:428934/76‑1 (MQ=255)
                   cAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCg                                                   <  2:864752/76‑1 (MQ=255)
                   cAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCg                                                   >  1:1228011/1‑76 (MQ=255)
                    aGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGt                                                  >  2:1868819/1‑76 (MQ=255)
                           tAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGtttt                                           <  2:1835131/76‑1 (MQ=255)
                           tAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGtttt                                           <  2:351875/76‑1 (MQ=255)
                           tAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGtttt                                           >  2:674244/1‑76 (MQ=255)
                            aGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTc                                          <  2:1768789/76‑1 (MQ=255)
                                     cTGAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                                 <  1:1138029/76‑1 (MQ=255)
                                       gAACACGCGCGCCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtttt                               <  2:1127730/76‑1 (MQ=255)
                                                cgcCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACgaga                      <  2:122749/76‑1 (MQ=255)
                                                 gcCACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGaa                     <  1:538044/76‑1 (MQ=255)
                                                  ccACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAc                    >  2:138832/1‑76 (MQ=255)
                                                   cACGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAcc                   <  2:35071/76‑1 (MQ=255)
                                                    aCGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAccc                  <  2:876386/76‑1 (MQ=255)
                                                     cGCATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCg                 >  1:921747/1‑76 (MQ=255)
                                                       cATACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGa               >  1:1406635/1‑76 (MQ=255)
                                                        aTACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGAt              >  2:385289/1‑76 (MQ=255)
                                                         tACCTACGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATc             <  2:1557665/76‑1 (MQ=255)
                                                              aCGCAAGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAg        >  1:1323360/1‑76 (MQ=255)
                                                                  aaGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgc    >  2:1724667/1‑76 (MQ=255)
                                                                  aaGCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgc    >  1:625042/1‑76 (MQ=255)
                                                                    gCACTTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgctg  >  1:676021/1‑76 (MQ=255)
                                                                        |                                                                       
TATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTG  >  REL606/3251944‑3252087

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