Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 473,901 +C 46.7% coding (1405/1677 nt) ybaL ← predicted transporter with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606473,9011.→C46.7% 0.0 ‑2.4 ‑0.7169coding (1405/1677 nt)ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  new base (37/41):  ref base (55/35):  total (92/77)
Fisher's exact test strand distribution p-value = 8.83e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.01e-01

Only 100 of 194 total aligned reads displayed.

TCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CGCATTGCCCAATACTGCGCGGACCCCGCGCTCT  >  REL606/473867‑473935
                                   |                                   
tCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGC                                      >  1:16964333‑M1/1‑35 (MQ=255)
tCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGC                                      <  1:15132258‑M1/35‑1 (MQ=255)
tCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGC                                      >  1:3246228‑M1/1‑35 (MQ=255)
tCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGC                                      <  1:6260556‑M1/35‑1 (MQ=255)
 ccAGATGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑C                                   <  1:10746977‑M1/35‑1 (MQ=255)
 ccAGATGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑C                                   >  1:5626243‑M1/1‑35 (MQ=255)
 ccAGATGCGCCAGTTGCATAATTTCTTCGTTAGC‑‑C                                   <  1:7530476‑M1/35‑1 (MQ=255)
  cAGATGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CG                                  >  1:14350651‑M1/1‑35 (MQ=255)
  cAGATGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CG                                  <  1:5555546‑M1/35‑1 (MQ=255)
  cAGATGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CG                                  >  1:2265118‑M1/1‑35 (MQ=255)
  cAGATGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CG                                  >  1:12862137‑M1/1‑35 (MQ=255)
   aaaTGCGCCAGTTTCATAATTTCTTCTTTcgc‑‑cgc                                 <  1:8216419/33‑1 (MQ=25)
   aGATGCGCCAGTTTCATAATTTCTTCGTTCGC‑‑CGC                                 >  1:8126251‑M1/1‑35 (MQ=255)
     aTGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgc                               >  1:16704037‑M1/1‑33 (MQ=255)
     aTGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgc                               >  1:4184225‑M1/1‑33 (MQ=255)
     aTGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgc                               >  1:20561706‑M1/1‑33 (MQ=255)
     aTGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgc                               >  1:11376272‑M1/1‑33 (MQ=255)
     aTGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgc                               >  1:10245524‑M1/1‑33 (MQ=255)
     aTGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CGCAt                               >  1:4457490/1‑35 (MQ=255)
       gcgcCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcat                             >  1:16187358‑M1/1‑31 (MQ=255)
       gcgcCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcat                             <  1:2536497‑M1/35‑5 (MQ=255)
       gcgcCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcat                             >  1:5908500‑M1/1‑31 (MQ=255)
       gcgcCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcat                             <  1:13948649‑M1/35‑5 (MQ=255)
        cgcCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcatt                            <  1:10150244‑M1/35‑6 (MQ=255)
        cgcCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcatt                            >  1:1680785‑M1/1‑30 (MQ=255)
        cgcCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcatt                            >  1:11295395‑M1/1‑30 (MQ=255)
         gcCAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattg                           >  1:12601666‑M1/1‑29 (MQ=255)
          ccAGTTTCATAATTTCTTCGTTCGC‑‑CGCgcattgc                          >  1:18923369‑M1/1‑28 (MQ=255)
          ccAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgc                          <  1:6301406‑M1/35‑8 (MQ=255)
          ccAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgc                          >  1:11666950‑M1/1‑28 (MQ=255)
          ccAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgc                          >  1:11870647‑M1/1‑28 (MQ=255)
           cAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgcc                         >  1:16739866‑M1/1‑27 (MQ=255)
           cAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgcc                         >  1:20040315‑M1/1‑27 (MQ=255)
           cAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgcc                         <  1:7805302‑M1/35‑9 (MQ=255)
           cAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgcc                         <  1:948228‑M1/35‑9 (MQ=255)
           cAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgcc                         >  1:13491868‑M1/1‑27 (MQ=255)
           cAGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgcc                         <  1:13012124‑M1/35‑9 (MQ=255)
            aGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgccc                        <  1:16044783‑M1/35‑10 (MQ=255)
            aGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgccc                        <  1:3328955‑M1/35‑10 (MQ=255)
            aGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgccc                        <  1:17589770‑M1/35‑10 (MQ=255)
            aGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgccc                        >  1:20097143‑M1/1‑26 (MQ=255)
            aGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgccc                        >  1:11044705‑M1/1‑26 (MQ=255)
            aGTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgccc                        >  1:6800492‑M1/1‑26 (MQ=255)
            aGTTGCATAATTTCTTCGTTCGC‑‑CGCATTGCCCaa                        >  1:3985835/1‑35 (MQ=255)
             gTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgccca                       >  1:3564444‑M1/1‑25 (MQ=255)
             gTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgccca                       >  1:5328282‑M1/1‑25 (MQ=255)
             gTTGCATAATTTCTTCGTTCGC‑‑CGCgcattgccca                       <  1:7676113‑M1/35‑11 (MQ=255)
              ttGCATAATTTCTTCGTTCGC‑‑CGCgcattgcccaa                      >  1:611723‑M1/1‑24 (MQ=255)
              ttGCATAATTTCTTCGTTCGC‑‑CGCgcattgcccaa                      >  1:5813837‑M1/1‑24 (MQ=255)
              ttGCATAATTTCTTCGTTCGC‑‑CGCgcattgcccaa                      >  1:15335595‑M1/1‑24 (MQ=255)
               tGCATAATTTCTTCGTTCGCCGCGCATTGCCCAAt                       <  1:1248537/35‑1 (MQ=21)
               tGCATAATTTCTTCGTTCGCCGCGCATTGCCCAAt                       >  1:11383328/1‑35 (MQ=21)
               tGCATAATTTCTTCGTTCGCCGCGCATTGCCCAAt                       >  1:17051941/1‑35 (MQ=21)
               tGCATAATTTCTTCGTTCGCCGCGCATTGCCCAAt                       >  1:4586943/1‑35 (MQ=21)
               tGCATAATTTCTTCGTTCGCCGCGCATTGCCCAAt                       <  1:640921/35‑1 (MQ=21)
                gCATAATTTCTTCGTTCGCCGCGCATTGCCCAATa                      <  1:3911167/35‑1 (MQ=21)
                gCATAATTTCTTCGTTCGCCGCGCATTGCCCAATa                      >  1:8660606/1‑35 (MQ=21)
                gCATAATTTCTTCGTTCGCCGCGCATTGCCCAATa                      >  1:8407940/1‑35 (MQ=21)
                 cATAATTTCTTCTTTCGCCGCGCATTGCCCAATAc                     <  1:692669/35‑1 (MQ=12)
                 cATAATTTCTTCGTTCGCCGCGCATTGCCCAATAc                     >  1:2642220/1‑35 (MQ=21)
                 cATAATTTCTTCGTTCGCCGCGCATTGCCCAATAc                     >  1:6528868/1‑35 (MQ=21)
                 cATAATTTCTTCGTTCGCCGCGCATTGCACAATan                     <  1:20319554/35‑2 (MQ=16)
                  aTAATTTCTTCGTTCGCCGCGCATTGCCCAATACt                    <  1:19976529/35‑1 (MQ=21)
                   tAATTTCTTCGTTCGG‑‑CGCATTGCCCAATACTgcg                 >  1:9672242/1‑35 (MQ=255)
                   tAATTTCTTCGTTCGCCGCGCATTGCCCAATACTg                   >  1:5683721/1‑35 (MQ=21)
                   tAATTTCTTCGTTCGCCGCGCATTGCCCAATACTg                   <  1:19594005/35‑1 (MQ=21)
                    aaTTTCTTCGTTCGCCGCGCATTGCCCAATACTgc                  >  1:9465597/1‑35 (MQ=21)
                    aaTTTCTTCGTTCGCCGCGCATTGCCCAATACTgc                  >  1:285157/1‑35 (MQ=21)
                    aaTTTCTTCGTTCGCCGCGCATTGCCCAATACTgc                  <  1:2285936/35‑1 (MQ=21)
                    aaTTTCTTCGTTCGCCGCGCATTGCCCAATACTgc                  <  1:13881239/35‑1 (MQ=21)
                    aaTTTCTTCGTTCGCCGCGCATTGCCCAATACTgc                  >  1:7297378/1‑35 (MQ=21)
                    aaTTTCTTCGTTCGCCGCGCATTGCCCAATACTgc                  >  1:19452066/1‑35 (MQ=21)
                    aaTTTCTTCGTTCGCCGCGCATTGCCAAATACTgc                  <  1:20247175/35‑1 (MQ=12)
                     aTTTCTTCGTTCGCCGCGCATTGCCCAATACTgcg                 <  1:10293009/35‑1 (MQ=21)
                     aTTTCTTCGTTCGCCGCGCATTGCCCAATACTgcg                 >  1:9128955/1‑35 (MQ=21)
                     aTTTCTTCGTTCGCCGCGCATTGCCCAATACTgcg                 >  1:15608167/1‑35 (MQ=21)
                     aTTTCTTCGTTCGCCGCGCATTGCCCAATACTgcg                 >  1:14903721/1‑35 (MQ=21)
                     aTTTCTTCGTTCGCCGCGCATTGCCCAATACTgcg                 >  1:17935425/1‑35 (MQ=21)
                      tttCTTCGTTCGCCGCGCATTGCCCAATACTgcgc                <  1:4191422/35‑1 (MQ=21)
                      tttCTTCGTTCGCCGCGCATTGCCCAATACTgcgc                >  1:17249795/1‑35 (MQ=21)
                       ttcttcGTTCGCCGCGCATTGCCCAATACTgcgcn               <  1:15873307/35‑2 (MQ=21)
                       ttcttcGTTCGCCGCGCATTGCCCAATACTGcgcg               <  1:8195988/35‑1 (MQ=21)
                        tcttcGTTCGCCGCGCATTGCCCAATACTGCGCgg              <  1:13726239/35‑1 (MQ=21)
                        tcttcGTTCGCCGCGCATTGCCCAATACTGCGCgg              <  1:16556428/35‑1 (MQ=21)
                         cttcGTTCGCCGCGCATTGCCCAATACTGCGCGGa             <  1:6219449/35‑1 (MQ=21)
                         cttcGTTCGCCGCGCATTGCCCAATACTGCGCGGa             <  1:14071038/35‑1 (MQ=21)
                         cttcGTTCGCCGCGCATTGCCCAATACTGCGCGGa             <  1:18183627/35‑1 (MQ=21)
                         cttcGTTCGCCGCGCATTGCCCAATACTGCGCGGa             >  1:19782683/1‑35 (MQ=21)
                          ttcgttcgCCGCGCATTGCCCAATACTGCGCGGAc            >  1:9805429/1‑35 (MQ=21)
                          ttcgttcgCCGCGCATTGCCCAATACTGCGCGGAc            >  1:5433473/1‑35 (MQ=21)
                          ttcgttcgCCGCGCATTGCCCAATACTGCGCGGAc            <  1:10028700/35‑1 (MQ=21)
                              ttcgC‑‑CGCATTGCCCAATACTGCGCGGACCCcgcg      <  1:5068881/35‑1 (MQ=255)
                              ttcgC‑‑CGCATTGCCCAATACTGCGCGGACCCcgcg      >  1:16934997/1‑35 (MQ=255)
                                cgc‑‑cgcgcattgccacccgtctgctgcgtcgagtg    <  1:17749964‑M1/35‑30 (MQ=255)
                                cgc‑‑cgcgcattgccaaccgtctgctgcgtcgagtg    >  1:2036906‑M1/1‑6 (MQ=255)
                                cgc‑‑cgcgcattgccaaccgtctgctgcgtcgagtg    >  1:19216754‑M1/1‑6 (MQ=255)
                                cgc‑‑cgcgcattgccaaccgtctgctgcgtcgagtg    >  1:10176853‑M1/1‑6 (MQ=255)
                                 gc‑‑cgcgcattgccaaccgtctgctgcgtcgagtgc   <  1:16615536‑M1/35‑31 (MQ=255)
                                 gc‑‑cgcgcattgccaaccgtctgctgcgtcgagtgc   >  1:1353857‑M1/1‑5 (MQ=255)
                                  c‑‑cgcgcattgccaaccgtctgctgcgtcgagtgcg  <  1:16092172‑M1/35‑32 (MQ=255)
                                   |                                   
TCCAGATGCGCCAGTTGCATAATTTCTTCGTTCGC‑‑CGCATTGCCCAATACTGCGCGGACCCCGCGCTCT  >  REL606/473867‑473935

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ 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.