Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,111,342 C→T 48.3% G467G (GGC→GGT metL → bifunctional aspartate kinase II/homoserine dehydrogenase II

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,111,3420C→T48.3% 204.6 ‑4.2 ‑0.0153G467G (GGC→GGTmetLbifunctional aspartate kinase II/homoserine dehydrogenase II
Reads supporting (aligned to +/- strand):  new base (27/47):  ref base (45/34):  total (72/81)
Fisher's exact test strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.42e-01

Only 100 of 162 total aligned reads displayed.

AAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCC  >  REL606/4111308‑4111376
                                  |                                  
aaaaaCGCATCGGCCTGGTATTGTTCGGTAAGGGc                                    >  1:2889684/1‑35 (MQ=255)
aaaaaCGCATCGGCCTGGTATTGTTCGGTAAGGGc                                    >  1:935309/1‑35 (MQ=255)
 aaaaCGCATCGGCCTGGTATTGTTCGGTAAGGGCa                                   >  1:183911/1‑35 (MQ=255)
 aaaaCGCATCGGCCTGGTATTGTTCGGTAAGGGCa                                   <  1:16691570/35‑1 (MQ=255)
 aaaaCGCATCGGCCTGGTATTGTTCGGTAAGGGCa                                   >  1:16028834/1‑35 (MQ=255)
 aaaaCGCATCGGCCTGGTATTGTTCGGTAAGGGCa                                   >  1:9842713/1‑35 (MQ=255)
  aaaCGCATCGGCCTGGTATTGTTCGGTAAGGGTaa                                  >  1:5468693/1‑35 (MQ=255)
  aaaCGCATCGGCCTGGTATTGTTCGGTAAGGGCaa                                  >  1:20338726/1‑35 (MQ=255)
  aaaCGCATCAGCCTGGTATTGTTCGGTAAGGGTaa                                  <  1:18537677/35‑1 (MQ=37)
   aaCGCATCGGCCTGGTATTGTTCGGTAAGGGTAat                                 <  1:13432580/35‑1 (MQ=255)
   aaCGCATCGGCCTGGTATTGTTCGGTAAGGGCAat                                 >  1:10924808/1‑35 (MQ=255)
   aaCGCATCGGCCTGGTATTGTTCGGTAAGGGCAat                                 >  1:17440205/1‑35 (MQ=255)
   aaCGCATCGGCCTGGTATTGTTCGGTAAGGGCAat                                 <  1:3935495/35‑1 (MQ=255)
   aaCGCATCGGCCTGGTATTGTTCGGTAAGGGCAat                                 >  1:8071610/1‑35 (MQ=255)
   aaCGCATCGGCCTGGTATTGTTCGGTAAGGGCAat                                 >  1:19911111/1‑35 (MQ=255)
    aCGCATCGGCCTGGTATTGTTCGGTAAGGGTAata                                <  1:7967420/35‑1 (MQ=255)
    aCGCATCGGCCTGGTATTGTTCGGTAAGGGTAata                                <  1:16925746/35‑1 (MQ=255)
    aCGCATCGGCCTGGTATTGTTCGGTAAGGGCAata                                <  1:9552378/35‑1 (MQ=255)
     cGCATCGGCCTGGTATTGTTCGGTAAGGGTAatat                               <  1:14476877/35‑1 (MQ=255)
     cGCATCGGCCTGGTATTGTTCGGTAAGGGTAatat                               >  1:8221359/1‑35 (MQ=255)
     cGCATCGGCCTGGTATTGTTCGGTAAGGGTAatat                               <  1:6896350/35‑1 (MQ=255)
     cGCATCGGCCTGGTATTGTTCGGTAAGGGCAatat                               <  1:4659422/35‑1 (MQ=255)
     cGCATCGGCCTGGTATTGTTCGGTAAGGGCAatat                               >  1:13696231/1‑35 (MQ=255)
     cGCATCGGCCTGGTATTGTTCGGTAAGGGCAatat                               <  1:4323799/35‑1 (MQ=255)
      gCATCGGCCTGGTATTGTTCGGTAAGGGTAATATc                              >  1:11554008/1‑35 (MQ=255)
      gCATCGGCCTGGTATTGTTCGGTAAGGGTAATATc                              <  1:6012806/35‑1 (MQ=255)
      gCATCGGCCTGGTATTGTTCGGTAAGGGCAATATc                              >  1:12774994/1‑35 (MQ=255)
      gCATCGGCCTGGTATTGTTCGGTAAGGGCAATATc                              >  1:17881406/1‑35 (MQ=255)
      gCATCGGCCTGGTATTGTTCGGTAAGGGCAATATc                              >  1:5016106/1‑35 (MQ=255)
       cATCGGCCTGGTATTGTTCGGTAAGGGTAATATCg                             <  1:2955974/35‑1 (MQ=255)
       cATCGGCCTGGTATTGTTCGGTAAGGGTAATATCg                             <  1:11181320/35‑1 (MQ=255)
       cATCGGCCTGGTATTGTTCGGTAAGGGCAATATCg                             <  1:6971226/35‑1 (MQ=255)
       cATCGGCCTGGTATTGTTCGGTAAGGGCAATATCg                             <  1:11881875/35‑1 (MQ=255)
        aTCGGCCTGGTATTGTTCGGTAAGGGTAATATCgg                            <  1:14349273/35‑1 (MQ=255)
        aTCGGCCTGGTATTGTTCGGTAAGGGCAATATCgg                            >  1:17605118/1‑35 (MQ=255)
         tCGGCCTGGTATTGTTCGGTAAGGGTAATATCggg                           >  1:18860792/1‑34 (MQ=255)
         tCGGCCTGGTATTGTTCGGTAAGGGTAATATCggc                           >  1:16010333/1‑34 (MQ=255)
         tCGGCCTGGTATTGTTCGGTAAGGGTAATATCGGt                           <  1:18447989/35‑1 (MQ=255)
         tCGGCCTGGTATTGTTCGGTAAGGGTAATATCGGt                           <  1:27465/35‑1 (MQ=255)
         tCGGCCTGGTATTGTTCGGTAAGGGTAATATCGGt                           >  1:10761465/1‑35 (MQ=255)
         tCGGCCTGGTATTGTTCGGTAAGGGCAATATCggc                           >  1:609701/1‑34 (MQ=255)
         tCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGt                           <  1:15368268/35‑1 (MQ=255)
           ggCCTGGTATTGTTCGGTAAGGGTAATATCGGTTc                         <  1:4505510/35‑1 (MQ=255)
           ggCCTGGTATTGTTCGGTAAGGGCAATATCGGTTc                         >  1:824189/1‑35 (MQ=255)
           ggCCTGGTATTGTTCGGTAAGGGCAATATCGGTTc                         <  1:12470455/35‑1 (MQ=255)
           ggCCTGGTATTGTTCGGTAAGGGCAATATCGGTTc                         <  1:15453355/35‑1 (MQ=255)
           ggCCTGGTATTGTTCGGTAAGGGCAATATCGGTTc                         <  1:3094262/35‑1 (MQ=255)
            gCCTGGTATTGTTCGGTAAGGGTAATATCGGTTcc                        <  1:6442531/35‑1 (MQ=255)
            gCCTGGTATTGTTCGGTAAGGGTAATATCGGTTcc                        <  1:5582746/35‑1 (MQ=255)
            gCCTGGTATTGTTCGGTAAGGGTAATATCGGTTcc                        >  1:17395973/1‑35 (MQ=255)
             ccTGGTATTGTTCGGTAAGGGTAATATCGGTTccc                       <  1:14098389/35‑1 (MQ=255)
             ccTGGTATTGTTCGGTAAGGGTAATATCGGTTccc                       >  1:15053343/1‑35 (MQ=255)
             ccTGGTATTGTTCGGTAAGGGCAATATCGGTTccc                       >  1:13627994/1‑35 (MQ=255)
             ccTGGTATTGTTCGGTAAGGGCAATATCGGTTccc                       >  1:12744152/1‑35 (MQ=255)
              cTGGTATTGTTCGGTAAGGGTAATATCGGTTCCCg                      <  1:18149585/35‑1 (MQ=255)
              cTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCg                      <  1:19312004/35‑1 (MQ=255)
              cTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCg                      >  1:19315177/1‑35 (MQ=255)
              cTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCg                      <  1:7901921/35‑1 (MQ=255)
              cTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCg                      >  1:9073051/1‑35 (MQ=255)
               tGGTATTGTTCGGTAAGGGTAATATCGGTTcccct                     >  1:1087609/1‑33 (MQ=38)
               tGGTATTGTTCGGTAAGGGTAATATCGGGTCCCGt                     >  1:4154780/1‑35 (MQ=37)
                ggTATTGTTCGGTAAGGGTAATATCGGTTCCCGtt                    <  1:7217479/35‑1 (MQ=255)
                ggTATTGTTCGGTAAGGGCAATATCGGTTCCCGtt                    <  1:6312317/35‑1 (MQ=255)
                 gTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTg                   >  1:10163267/1‑35 (MQ=255)
                 gTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTg                   >  1:12268415/1‑35 (MQ=255)
                  tATTGTTCGGTAAGGGTAATATCGGTTCCCGTTgg                  <  1:19580779/35‑1 (MQ=255)
                  tATTGTTCGGTAAGGGCAATATCGGTTCCCTTTgg                  >  1:5300679/1‑35 (MQ=255)
                   aTTGTTCGGTAAGGGTAATATCGGTTCCCGTTGGc                 <  1:10643048/35‑1 (MQ=255)
                   aTTGTTCGGTAAGGGTAATATCGGTTCCCGTTGGc                 >  1:7122099/1‑35 (MQ=255)
                   aTTGTTCGGTAAGGGTAATATCGGTTCCCGTTGGc                 <  1:2458157/35‑1 (MQ=255)
                   aTTGTTCGGTAAGGGTAATATCGGTTCCCGTTGGc                 >  1:19721947/1‑35 (MQ=255)
                   aTTGTGCGGTAAGGGCAATATCGGTTCCCGTTGGc                 <  1:14992790/35‑1 (MQ=255)
                    ttGTTCGGTAAGGGTAATATCGGTTCCCGTTGGCt                >  1:19783360/1‑35 (MQ=255)
                    ttGTTCGGTAAGGGTAATATCGGTTCCCGTTGGCt                <  1:4484063/35‑1 (MQ=255)
                    ttGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCt                <  1:3944275/35‑1 (MQ=255)
                     tGTTCTGTAAGGGTAATATCGGTTCCCGGTggggg               >  1:13896356/1‑32 (MQ=255)
                     tGTTCGGTAAGGGTAATATCGGTTCCCGTTGGCTg               <  1:7639334/35‑1 (MQ=255)
                      gTTCGGTAAGGGTAATATCGGTTCCCCTTGGCTgg              >  1:11897790/1‑35 (MQ=255)
                      gTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTgg              <  1:10244821/35‑1 (MQ=255)
                       ttCGGTAAGGGTAATATCGGTTCCCGTTGGCTGGa             <  1:14734954/35‑1 (MQ=255)
                       ttCGGTAAGGGCAATATCGGGTCCCGTTGGGTggg             >  1:16277534/1‑34 (MQ=37)
                        tCGGTAAGGGTAATATCGGTTCCCGTTGGCTGGaa            <  1:9538040/35‑1 (MQ=255)
                        tCGGTAAGGGTAATATCGGTTCCCGTTGGCTGGaa            <  1:5905733/35‑1 (MQ=255)
                        tCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGaa            <  1:14464690/35‑1 (MQ=255)
                         cGGTAAGGGTAATATCGGTTCCCGTTGGCTGGAAc           <  1:1197846/35‑1 (MQ=255)
                         cGGTAAGGGTAATATCGGTTCACGTTGGCTGGAAc           <  1:13581788/35‑1 (MQ=38)
                          ggTAAGGGCAATATCGGTTCCCGTTGGCTGGGACt          >  1:14778416/1‑35 (MQ=255)
                          ggTAAGGGCAATATCGGTTCCCGTTGGCTGGAACt          <  1:13081339/35‑1 (MQ=255)
                          ggTAAGGGCAATATCGGTGCCCGTTGGCTGGAACt          <  1:20615331/35‑1 (MQ=255)
                           gTAAGGGCAATATCGGTTCCCGTTGGCTGGGACTg         >  1:13439280/1‑35 (MQ=255)
                           gTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTg         <  1:17866812/35‑1 (MQ=255)
                             aaGGGTAATATCGGTTCCCGTTGGCTGGAACTGtt       <  1:1886070/35‑1 (MQ=255)
                             aaGGGCAATATCGGTTCCCGTTGGCTGGGACTGtt       >  1:10605146/1‑35 (MQ=255)
                              aGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTc      <  1:7574453/35‑1 (MQ=255)
                              aGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTc      <  1:14005445/35‑1 (MQ=255)
                               gggTAATATCGGTTCCCGTTGGCTGGAACTGTTCg     <  1:7290125/35‑1 (MQ=255)
                               gggTAATATCGGTTCCCGTTGGCTGGAACTGTTCg     <  1:14284682/35‑1 (MQ=255)
                                ggCAATATCTGTTCCCGTTGGCTGGAACTGTTCGc    >  1:7740221/1‑35 (MQ=255)
                                ggCAATATCGGTTCCCGTTGGCTGGAACTGTTCGc    <  1:10185748/35‑1 (MQ=255)
                                  cAATATCGGTTCCCGTTGGCTGGGACTGTTCGccc  >  1:9147126/1‑35 (MQ=255)
                                  |                                  
AAAAACGCATCGGCCTGGTATTGTTCGGTAAGGGCAATATCGGTTCCCGTTGGCTGGAACTGTTCGCCC  >  REL606/4111308‑4111376

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