Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,152,406 A→C 21.7% K41N (AAA→AAC ygjG → putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,152,4060A→C21.7% 118.3 ‑0.3 ‑0.2171K41N (AAA→AACygjGputrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  new base (16/21):  ref base (63/71):  total (79/92)
Fisher's exact test strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.88e-01

Only 100 of 185 total aligned reads displayed.

TTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAG  >  REL606/3152372‑3152440
                                  |                                  
 tGAGAAGCGAACGCTGGATCATGAGGAGATGAAAg                                   >  1:20099457/1‑35 (MQ=255)
 tGAGAAGCGAACGCTGGATCATGAGGAGATGAAAg                                   <  1:1334289/35‑1 (MQ=255)
  gagaAGCGAACGCTGGATCATGAGGAGATGAACGc                                  <  1:8995197/35‑1 (MQ=255)
  gagaAGCGAACGCTGGATCATGAGGAGATGAAAGc                                  >  1:3866883/1‑35 (MQ=255)
  gagaAGCGAACGCTGGATCATGAGGAGATGAAAGc                                  <  1:19816295/35‑1 (MQ=255)
  gagaAGCGAACGCTGGATCATGAGGAGATGAAAGc                                  <  1:8193061/35‑1 (MQ=255)
   agaAGCGAACGCTGGATCATGAGGAGATGAACGCa                                 <  1:18682840/35‑1 (MQ=255)
   agaAGCGAACGCTGGATCATGAGGAGATGAAAGCa                                 >  1:6937796/1‑35 (MQ=255)
   agaAGCGAACGCTGGATCATGAGGAGATGAAAGCa                                 <  1:18628238/35‑1 (MQ=255)
   agaAGCGAACGCTGGATCATGAGGAGATGAAAGCa                                 <  1:1781979/35‑1 (MQ=255)
    gaAGCGAACGCTGGATCATGAGGAGATGAAAGCAc                                <  1:6337725/35‑1 (MQ=255)
    gaAGCGAACGCTGGATCATGAGGAGATGAAAGCAc                                >  1:17965843/1‑35 (MQ=255)
     aaGCGAACGCTGGATCATGAGGAGATGAACGCACt                               >  1:11003341/1‑35 (MQ=255)
      aGCGAACGCTGGATCATGAGGAGATGAAAGCACtt                              <  1:2177121/35‑1 (MQ=255)
      aGCGAACGCTGGATCATGAGGAGATGAAAGCACtt                              <  1:8748191/35‑1 (MQ=255)
      aGCGAACGCTGGATCATGAGGAGATGAAAGCACtt                              <  1:11013591/35‑1 (MQ=255)
      aGCGAACGCTGGATCATGAGGAGATGAAAGCACtt                              <  1:19154066/35‑1 (MQ=255)
       gCGAACGCTGGATCATGAGGAGATGAAAGCACTTa                             >  1:8712088/1‑35 (MQ=255)
        cGAACGCTGGATCATGAGGAGATGAAAGCACTTaa                            >  1:10366821/1‑35 (MQ=255)
        cGAACGCTGGATCATGAGGAGATGAAAGCACTTaa                            <  1:8445481/35‑1 (MQ=255)
        cGAACGCTGGATCATGAGGAGATGAAAGCACTTaa                            <  1:18218699/35‑1 (MQ=255)
        cGAACGCTGGATCATGAGGAGATGAAAGCACTTaa                            <  1:14008781/35‑1 (MQ=255)
         gAACGCTGGATCATGAGGAGATGAAAGCACTTAAc                           <  1:19467639/35‑1 (MQ=255)
          aaCGCTGGATCATGAGGAGATGAACGCACTTAAcc                          >  1:9187481/1‑35 (MQ=255)
           aCGCTGGATCATTAGGAGATGAAAGCACTTAACCg                         >  1:5272270/1‑35 (MQ=255)
           aCGCTGGATCATGAGGAGATGAACGCACTTAACCg                         >  1:11037351/1‑35 (MQ=255)
           aCGCTGGATCATGAGGAGATGAAAGCACTTAACCg                         <  1:19394041/35‑1 (MQ=255)
           aCGCTGGATCATGAGGAGATGAAAGCACTTAACCg                         >  1:13037650/1‑35 (MQ=255)
            cGCTGGATCATGAGGAGATGAAAGCACTTAACCga                        <  1:6997052/35‑1 (MQ=255)
            cGCTGGATCATGAGGAGATGAAAGCACTTAACCga                        <  1:17024800/35‑1 (MQ=255)
            cGCTGGATCATGAGGAGATGAAAGCACTTAACCga                        <  1:9674233/35‑1 (MQ=255)
             gCTGGATCATGAGGAGATGAACGCACTTAACCgag                       <  1:1645349/35‑1 (MQ=255)
             gCTGGATCATGAGGAGATGAAAGCACTTAACCgag                       <  1:657431/35‑1 (MQ=255)
             gCTGGATCATGAGGAGATGAAAGCACTTAACCgag                       >  1:10674813/1‑35 (MQ=255)
             gCTGGATCATGAGGAGATGAAAGCACTTAACCgag                       >  1:10217163/1‑35 (MQ=255)
             gCTGGATCATGAGGAGATGAAAGCACTTAACCgag                       >  1:4050183/1‑35 (MQ=255)
             gCTGGATCATGAGGAGATGAAAGCACTTAACCgag                       >  1:2330273/1‑35 (MQ=255)
              cTGGATCATGAGGAGATGAAAGCACTTAACCgaga                      <  1:3807351/35‑1 (MQ=255)
              cTGGATCATGAGGAGATGAAAGCACTTAACCgaga                      <  1:2150064/35‑1 (MQ=255)
               tGGATCATGAGGAGATGAAAGCACTTAACCGagag                     >  1:13922976/1‑35 (MQ=255)
               tGGATCATGAGGAGATGAAAGCACTTAACCGagag                     >  1:7665540/1‑35 (MQ=255)
                ggATCATGAGGAGATGAAAGCACTTAACCGAGAgg                    >  1:5786728/1‑35 (MQ=255)
                ggATCATGAGGAGATGAAAGCACTTAACCGAGAgg                    >  1:11328258/1‑35 (MQ=255)
                ggATCATGAGGAGATGAAAGCACTTAACCGAGAgg                    <  1:16008952/35‑1 (MQ=255)
                ggATCATGAGGAGATGAAAGCACTTAACCGAGAgg                    >  1:16096211/1‑35 (MQ=255)
                 gATCATGAGGAGATGAACGCACTTAACCGAGAGGt                   <  1:14216623/35‑1 (MQ=255)
                 gATCATGAGGAGATGAACGCACTTAACCGAGAGGt                   >  1:9019225/1‑35 (MQ=255)
                 gATCATGAGGAGATGAACGCACTTAACCGAGAGGt                   >  1:13279668/1‑35 (MQ=255)
                 gATCATGAGGAGATGAAAGCACTTAACCGAGAGGt                   >  1:15337964/1‑35 (MQ=255)
                  aTCATGAGGAGATGAAAGCACTTAACCGAGAGGTg                  >  1:150555/1‑35 (MQ=255)
                  aTCATGAGGAGATGAAAGCACTTAACCGAGAGGTg                  <  1:7710883/35‑1 (MQ=255)
                   tCATGAGGAGATGAAAGCACTTAACCGAGAGGTGa                 <  1:16225624/35‑1 (MQ=255)
                    cATGAGGAGATGAAAGCACTTAACCGAGAGGTTAt                >  1:14799251/1‑35 (MQ=255)
                    cATGAGGAGATGAAAGCACTTAACCGAGAGGTGAt                <  1:9738300/35‑1 (MQ=255)
                    cATGAGGAGATGAAAGCACTTAACCGAGAGGTGAt                >  1:20479187/1‑35 (MQ=255)
                    cATGAGGAGATGAAAGCACTTAACCGAGAGGTGAt                >  1:4215041/1‑35 (MQ=255)
                    cATGAGGAGATGAAAGCACTTAACCGAGAGGTGAt                >  1:2792708/1‑35 (MQ=255)
                     aTGAGGAGATGAACGCACTTAACCGAGAGGTGAtt               <  1:17483087/35‑1 (MQ=255)
                     aTGAGGAGATGAACGCACTTAACCGAGAGGTGAtt               <  1:2436164/35‑1 (MQ=255)
                     aTGAGGAGATGAACGCACTTAACCGAGAGGTGAtt               <  1:4338440/35‑1 (MQ=255)
                     aTGAGGAGATGAACGCACTTAACCGAGAGGTGAtt               <  1:17098952/35‑1 (MQ=255)
                     aTGAGGAGATGAAAGCACTTAACCGAGAGGTGAtt               <  1:6002433/35‑1 (MQ=255)
                     aTGAGGAGATGAAAGCACTTAACCGAGAGGTGAtt               <  1:15173662/35‑1 (MQ=255)
                      tGAGGAGATGAACGCACTTAACCGAGAGGTGATTg              <  1:12457515/35‑1 (MQ=255)
                      tGAGGAGATGAAAGCACTTAACCGAGATGTGATTg              >  1:4372375/1‑35 (MQ=255)
                      tGAGGAGATGAAAGCACTTAACCGAGAGGTGATTg              <  1:16646011/35‑1 (MQ=255)
                       gaggagATGAACGCACTTAACCGAGAGGTGATTGa             >  1:4761066/1‑35 (MQ=255)
                       gaggagATGAAAGCACTTAACCGAGAGGTGATTGa             <  1:19096914/35‑1 (MQ=255)
                       gaggagATGAAAGCACTTAACCGAGAGGTGATTGa             <  1:12214134/35‑1 (MQ=255)
                       gaggagATGAAAGCACTTAACCGAGAGGTGATTGa             <  1:2659549/35‑1 (MQ=255)
                        aggagATGAAAGCACTTAACCGAGAGGTGATTGaa            <  1:6403629/35‑1 (MQ=255)
                        aggagATGAAAGCACTTAACCGAGAGGTGATTGaa            <  1:19902381/35‑1 (MQ=255)
                         ggagATGAAAGCACTTAACCGAGAGGTGATTGAAt           <  1:14509307/35‑1 (MQ=255)
                         ggagATGAAAGCACTTAACCGAGAGGTGATTGAAt           <  1:17166712/35‑1 (MQ=255)
                         ggagATGAAAGCACTTAACCGAGAGGTGATTGAAt           >  1:11512589/1‑35 (MQ=255)
                          gagaTGAACGCACTTAACCGAGAGGTGATTGAATa          <  1:3414001/35‑1 (MQ=255)
                          gagaTGAACGCACTTAACCGAGAGGTGATTGAATa          <  1:857097/35‑1 (MQ=255)
                          gagaTGAAAGCACTTAACCGAGAGGTGATTGAATa          >  1:17056612/1‑35 (MQ=255)
                           agaTGAACGCACTTAACCGAGAGGTGATTGAATAc         >  1:3483703/1‑35 (MQ=255)
                           agaTGAACGCACTTAACCGAGAGGTGATTGAATAc         <  1:4975838/35‑1 (MQ=255)
                            gaTGAAAGCACTTAACCGAGAGGTGATTGAATACt        >  1:12964379/1‑35 (MQ=255)
                            gaTGAAAGCACTTAACCGAGAGGTGATTGAATACt        <  1:17300109/35‑1 (MQ=255)
                            gaTGAAAGCACTTAACCGAGAGGTGATTGAATACt        <  1:18117984/35‑1 (MQ=255)
                             aTGAACGCACTTAACCGAGAGGTGATTGAATACtt       >  1:12098091/1‑35 (MQ=255)
                             aTGAAAGCACTTAACCGAGAGGTGATTGAATACtt       >  1:5574378/1‑35 (MQ=255)
                             aTGAAAGCACTTAACCGAGAGGTGATTGAATACtt       <  1:1886965/35‑1 (MQ=255)
                             aTGAAAGCACTTAACCGAGAGGTGATTGAATACtt       <  1:17902617/35‑1 (MQ=255)
                             aTGAAAGCACTTAAACGAGAGGTGATTGAATACtt       <  1:3669890/35‑1 (MQ=255)
                              tGAAAGCACTTAACCGAGAGGTGATTGAATACTTc      <  1:3042874/35‑1 (MQ=255)
                              tGAAAGCACTTAACCGAGAGGTGATTGAATACTTc      <  1:14827431/35‑1 (MQ=255)
                               gAAAGCACTTAACCGAGAGGTGATTGAATACTTCa     <  1:5006801/35‑1 (MQ=255)
                               gAAAGCACTTAACCGAGAGGTGATTGAATACTTCa     >  1:18216735/1‑35 (MQ=255)
                               gAAAGCACTTAACCGAGAGGTGATTGAATACTTCa     >  1:17462022/1‑35 (MQ=255)
                               gAAAGCACTTAACCGAGAGGTGATTGAATACTTCa     >  1:14249282/1‑35 (MQ=255)
                               gAAAGCACTTAACCGAGAGGTGATTGAATACTTCa     <  1:11776259/35‑1 (MQ=255)
                               gAAAGCACTTAACCGAGAGGTGATTGAATACTTCa     <  1:10610035/35‑1 (MQ=255)
                                aaaGCACTTAACCGAGAGGTGATTGAATACTTCaa    >  1:5317306/1‑35 (MQ=255)
                                aaaGCACTTAACCGAGAGGTGATTGAATACTTCaa    >  1:15681077/1‑35 (MQ=255)
                                 aaGCACTTAACCGAGAGGTGATTGAATACTTCaaa   >  1:10544017/1‑35 (MQ=255)
                                 aaGCACTTAACCGAGAGGTGATTGAATACTTCaaa   <  1:10215607/35‑1 (MQ=255)
                                  |                                  
TTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAG  >  REL606/3152372‑3152440

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