Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,616,538 A→C 94.2% Y337S (TAC→TCC)  nadR → nicotinamide‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,616,5380A→C94.2% 13.1 ‑0.3 ‑0.5182Y337S (TAC→TCC) nadRnicotinamide‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  new base (84/87):  ref base (6/5):  total (90/92)
Fisher's exact test strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.84e-01

Only 100 of 183 total aligned reads displayed.

GCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAAC  >  REL606/4616506‑4616569
                                |                               
gCATCCGTTCGTACAGGCGTTGATTGATGAATCcc                               >  1:14543107/1‑35 (MQ=255)
gCATCCGTTCGTACAGGCGTTGATTGATGAATCcc                               <  1:12255550/35‑1 (MQ=255)
 cATCCGTTCGTACAGTCGTTGATTGATGAATCCCg                              >  1:1541508/1‑35 (MQ=38)
 cATCCGTTCGTACAGGCGTTGATTGATGAATCCCg                              <  1:12802035/35‑1 (MQ=255)
 cATCCGTTCGTACAGGCGTTGATTGATGAATCCCg                              >  1:12548944/1‑35 (MQ=255)
  aTCCGTTCGTACAGGCGTTGATTGATGAATCCCGt                             <  1:14396070/35‑1 (MQ=255)
  aTCCGTTCGTACAGGCGTTGATTGATGAATCCCGt                             <  1:1599789/35‑1 (MQ=255)
  aTCCGTTCGTACAGGCGTTGATTGATGAATCCCGt                             >  1:18583314/1‑35 (MQ=255)
  aTCCGTTCGTACAGGCGTTGATTGATGAATCCCGt                             >  1:19490857/1‑35 (MQ=255)
  aTCCGTTCGTACAGGCGTTGATTGATGAATCCCGt                             >  1:7935974/1‑35 (MQ=255)
   tCCGTTCGTACAGGCGTTTATTGATGAATCCCGtt                            >  1:987457/1‑35 (MQ=37)
    ccGTTCGTACAGGCGTTGATTGATGAATCCCGttt                           <  1:7737656/35‑1 (MQ=255)
     cGTTCGTACAGGCGTTGATTGATGAATCCCGTTTc                          <  1:3344161/35‑1 (MQ=255)
     cGTTCGTACAGGCGTTGATTGATGAATCCCGTTTc                          >  1:787177/1‑35 (MQ=255)
     cGTTCGTACAGGCGTTGATTGATGAATCCCGTTTc                          <  1:9303661/35‑1 (MQ=255)
      gTTCGTACAGGCGTTGATTGATGAATCCCGTTTCg                         >  1:6650389/1‑35 (MQ=255)
      gTTCGTACAGGCGTTGATTGATGAATCCCGTTTCg                         <  1:8751744/35‑1 (MQ=255)
      gTTCGTACAGGCGTTGATTGATGAATCCCGTTTCg                         <  1:16719788/35‑1 (MQ=255)
      gTTCGTACAGGCGTTGATTGATGAATCCCGTTTCg                         <  1:17033780/35‑1 (MQ=255)
      gTTCGTACAGGCGTTGATTGATGAATACCGTTTCg                         >  1:1041502/1‑35 (MQ=255)
       ttggTACAGGCGTTGATTGATGAATACCGTTTCGa                        <  1:15686211/32‑1 (MQ=255)
       ttCGTACAGGCGTTGATTGATGAATCCCGTTTCGa                        <  1:9773032/35‑1 (MQ=255)
       ttCGTACAGGCGTTGATTGATGAATCCCGTTTCGa                        <  1:17522495/35‑1 (MQ=255)
       ttCGTACAGGCGTTGATTGATGAATCCCGTTTCGa                        >  1:12595867/1‑35 (MQ=255)
       ttCGTACAGGCGTTGATTGATGAATCCCGTTTCGa                        >  1:2916400/1‑35 (MQ=255)
       ttCGTACAGGCGTTGATTGATGAATCCCGTTTCGa                        >  1:11865796/1‑35 (MQ=255)
        tCGTACAGGCGTTGATTGATGAATCCCGTTTCGAt                       <  1:18228505/35‑1 (MQ=255)
        tCGTACAGGCGTTGATTGATGAATCCCGTTTCGAt                       >  1:5016028/1‑35 (MQ=255)
        tCGTACAGGCGTTGATTGATGAATCCCGTTTCGAt                       >  1:10260038/1‑35 (MQ=255)
        tCGTACAGGCGTTGATTGATGAATCCCGTTTCGAt                       >  1:10907948/1‑35 (MQ=255)
         cGTACAGGCGTTGATTGATGAATCCCGTTTCGATc                      <  1:9688178/35‑1 (MQ=255)
         cGTACAGGCGTTGATTGATGAATCCCGTTTCGATc                      >  1:4545736/1‑35 (MQ=255)
         cGTACAGGCGTTGATTGATGAATCCCGTTTCGATc                      <  1:19910220/35‑1 (MQ=255)
         cGTACAGGCGTTGATTGATGAATCCCGTTTCGATc                      <  1:6504969/35‑1 (MQ=255)
         cGTACAGGCGTTGATTGATGAATCCCGTTTCGATc                      <  1:16253376/35‑1 (MQ=255)
         cGTACAGGCGTTGATTGATGAATCCAGTTTCGATc                      <  1:2703947/35‑1 (MQ=37)
          gTACAGGCGTTGATTGATGAATCCCGTTTCGATCt                     <  1:7616333/35‑1 (MQ=255)
           tACAGGCGTTGATTGATGAATCCCGTTTCGATCTg                    <  1:12041388/35‑1 (MQ=255)
           tACAGGCGTTGATTGATGAATCCCGTTTCGATCTg                    >  1:2111421/1‑35 (MQ=255)
           tACAGGCGTTGATTGATGAATCCCGTTTCGATCTg                    >  1:20398986/1‑35 (MQ=255)
           tACAGGCGTTGATTGATGAATCCCGTTTCGATCTg                    <  1:5040718/35‑1 (MQ=255)
           tACAGGCGTTGATTGATGAATCCCGTTTCGATCTg                    >  1:14099792/1‑35 (MQ=255)
            aCAGGCGTTGATTGATGAATCCCGTTTCGATCTgg                   <  1:16719095/35‑1 (MQ=255)
             cAGGCGTTGATTGATGAATCCCGTTTCGATCTGGt                  >  1:12071957/1‑35 (MQ=255)
              aGGCGTTGATTGATGAATCCCGTTTCGATCTGGTg                 >  1:3676039/1‑35 (MQ=255)
              aGGCGTTGATTGATGAATCCCGTTTCGATCTGGTg                 >  1:425955/1‑35 (MQ=255)
              aGGCGTTGATTGATGAATCCCGTTTCGATCTGGTg                 <  1:15169639/35‑1 (MQ=255)
              aGGCGTTGATTGATGAATCCCGTTTCGATCTGGTg                 <  1:19860996/35‑1 (MQ=255)
              aGGCGTTGATTGATGAATCCCGTTTCGATCTGGTg                 <  1:12235519/35‑1 (MQ=255)
               ggCGTTGATTGATGAATCCCGTTTCGATCTGGTGa                <  1:10767903/35‑1 (MQ=255)
               ggCGTTGATTGATGAATCCCGTTTCGATCTGGTGa                >  1:13114951/1‑35 (MQ=255)
                gCGTTGATTGATGAATCCCGTTTCGATCTGGTGAt               <  1:2606008/35‑1 (MQ=255)
                gCGTTGATTGATGAATACCGTTTCGATCTGGTGAt               >  1:12751287/1‑35 (MQ=255)
                 cGTTGATTGATGAATCCCGTTTCGATCTGTTGATc              >  1:12959333/1‑35 (MQ=37)
                 cGTTGATTGATGAATCCCGTTTCGATCTGGTGATc              >  1:19184432/1‑35 (MQ=255)
                 cGTTGATTGATGAATCCCGTTTCGATCTGGTGATc              >  1:12138263/1‑35 (MQ=255)
                 cGTTGATTGATGAATCCCGTTTCGATCTGGTGATc              <  1:1967990/35‑1 (MQ=255)
                 cGTTGATTGATGAATCCCGTTTCGATCTGGTGATc              >  1:14687671/1‑35 (MQ=255)
                 cGTTGATTGATGAATCCCGTTTCGATCTGGTGATc              <  1:2180736/35‑1 (MQ=255)
                 cGTTGATTGATGAATCCCGTTTCGATCTGGTGATc              >  1:16587938/1‑35 (MQ=255)
                 cGTTGATTGATGAATCCCGTTTAGATCTGGTGATc              <  1:815053/35‑1 (MQ=38)
                  gTTGATTGATGAATCCCGTTTCGATCTGGTGATcc             <  1:8540831/35‑1 (MQ=255)
                  gTTGATTGATGAATCCCGTTTCGATCTGGTGATcc             >  1:2984338/1‑35 (MQ=255)
                  gTTGATTGATGAATCCCGTTTCGATCTGGTGATcc             >  1:15391858/1‑35 (MQ=255)
                   ttgattgaTGAATCCCGTTTCGATCTGGTGATCct            >  1:17164385/1‑35 (MQ=255)
                   ttgattgaTGAATCCCGTTTCGATCTGGTGATCct            <  1:9190602/35‑1 (MQ=255)
                   ttgattgaTGAATCCCGTTTCGATCTGGTGATCct            <  1:7337197/35‑1 (MQ=255)
                   ttgattgaTGAATCCCGTTTCGATCTGGTGATCct            >  1:7300607/1‑35 (MQ=255)
                    ttattgatGAATCCCGTTTCGATCTGGTGATCctg           <  1:16507873/33‑1 (MQ=38)
                    tgattgatGAATCCCGTTTCGATCTGGTGATCctg           >  1:9143626/1‑35 (MQ=255)
                    tgattgatGAATCCCGTTTCGATCTGGTGATCctg           <  1:7166226/35‑1 (MQ=255)
                    tgattgatGAATCCCGTTTCGATCTGGTGATCctg           <  1:2540188/35‑1 (MQ=255)
                     gattgatGAATCCCGTTTCGATCTGGTGATCctgc          >  1:13696955/1‑35 (MQ=255)
                     gattgatGAATCCCGTTTCGATCTGGTGATCctgc          >  1:11497288/1‑35 (MQ=255)
                     gattgatGAATCCCGTTTCGATCTGGTGATCctgc          >  1:18039345/1‑35 (MQ=255)
                     gattgatGAATCCCGTTTCGATCTGGTGATCctgc          <  1:18938557/35‑1 (MQ=255)
                     gattgatGAATCCCGTTTCGATCTGGTGATCctgc          <  1:5949636/35‑1 (MQ=255)
                     gatcgatGAATCCCGTTTCGATCTGGTGATCctgc          >  1:17935054/1‑35 (MQ=37)
                      attgatGAATCCCGTTTCGATCTGGTGATCctgct         >  1:14308615/1‑35 (MQ=255)
                      attgatGAATCCCGTTTCGATCTGGTGATCctgct         <  1:13533742/35‑1 (MQ=255)
                      attgatGAATCCCGTTTCGATCTGGTGATCctgct         >  1:5333232/1‑35 (MQ=255)
                      attgatGAATACCGTTTCGATCTGGTGATCctgct         >  1:387071/1‑35 (MQ=255)
                       ttgatTAATCCCGTTTCGATCTGGTGATCctgctg        >  1:1853348/1‑35 (MQ=37)
                       ttgatGAATCCCGTTTCGATCTGGTGATCctgctg        <  1:6517953/35‑1 (MQ=255)
                       ttgatGAATCCCGTTTCGATCTGGTGATCctgctg        <  1:11713404/35‑1 (MQ=255)
                       ttgatGAATCCCGTTTCGATCTGGTGATCctgctg        <  1:1892525/35‑1 (MQ=255)
                        tgatgaATCCCGTTTCGATCTGGTGATCCTGCTgg       >  1:3069346/1‑35 (MQ=255)
                        tgatgaATCCCGTTTCGATCTGGTGATCCTGCTgg       >  1:17619655/1‑35 (MQ=255)
                        tgatgaATCCCGTTTCGATCTGGTGATCCTGCTgg       >  1:11965761/1‑35 (MQ=255)
                         gatgaATCCCGTTTCGATCTGGTGATCCTGCTGga      >  1:11627000/1‑35 (MQ=255)
                         gatgaATCCCGTTTCGATCTGGTGATCCTGCTGga      >  1:13045269/1‑35 (MQ=255)
                         gatgaATCCCGTTTCGATCTGGTGATCCTGCTGga      >  1:1547247/1‑35 (MQ=255)
                         gatgaATCCCGTTTCGATCTGGTGATCCTGCTGga      >  1:14860938/1‑35 (MQ=255)
                          atgaATCCCGTTTCGATCTGGTGATCCTGCTGgag     >  1:2296320/1‑35 (MQ=255)
                          atgaATCCCGTTTCGATCTGGTGATCCTGCTGgag     <  1:13930193/35‑1 (MQ=255)
                          atgaATCCCGTTTCGATCTGGTGATCCTGCTGgag     <  1:1057959/35‑1 (MQ=255)
                           tgaATCCCGTTTCGATCTGGTGATCCTGCTGgaga    <  1:10699568/35‑1 (MQ=255)
                            gaATCCCGTTTCGATCTGGTGATCCTGCTGGAGaa   >  1:19288501/1‑35 (MQ=255)
                             aaTCCCGTTTCGATCTGGTGATCCTGCTGGAGaac  >  1:3435218/1‑35 (MQ=255)
                             aaTCCCGTTTCGATCTGGTGATCCTGCTGGAGaac  >  1:8106661/1‑35 (MQ=255)
                                |                               
GCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAAC  >  REL606/4616506‑4616569

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.