Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,126,975 T→C L515P (CTG→CCG)  mdoH → glucosyltransferase MdoH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,126,9750TC97.8% 172.7 / ‑3.7 45L515P (CTG→CCG) mdoHglucosyltransferase MdoH
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/19);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

TGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGT  >  REL606/1126902‑1127047
                                                                         |                                                                        
tGCCACGGTATCCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGa                                                                        <  2:382059/76‑1 (MQ=255)
tGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGa                                                                        <  2:166792/76‑1 (MQ=255)
   cACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACgg                                                                     >  1:442197/1‑76 (MQ=255)
    aCGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTTTTCCCGCCggt                                                                    >  2:1191226/1‑75 (MQ=255)
    aCGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACggg                                                                    <  2:546577/76‑1 (MQ=255)
       gTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGt                                                                 <  2:173824/76‑1 (MQ=255)
               gtgaACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCtt                                                         <  1:960105/75‑1 (MQ=255)
               gtgaACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCtt                                                         <  2:518736/75‑1 (MQ=255)
               atgaACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCtt                                                         <  2:1266031/76‑1 (MQ=255)
               atgaACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCtt                                                         <  1:943392/76‑1 (MQ=255)
               atgaACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCtt                                                         >  1:1516917/1‑76 (MQ=255)
               atgaACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATATCtt                                                         >  2:1333759/1‑76 (MQ=255)
                tgaTCTTCCGTCTGTTCCTGGTGGAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTa                                                        <  2:230072/76‑1 (MQ=255)
                tgaACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTa                                                        <  2:750402/76‑1 (MQ=255)
                    cTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTAtctc                                                    <  1:547597/76‑1 (MQ=255)
                      tCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTCAtctctc                                                  >  2:611721/1‑76 (MQ=255)
                       ccGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTctcc                                                 >  2:1208793/1‑76 (MQ=255)
                        cGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTctccg                                                >  1:1198615/1‑76 (MQ=255)
                          tCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGTTGTTCCCGACGGGCGTGATGTCTTATCTctccgct                                              <  2:1391979/76‑1 (MQ=255)
                          tCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTctccgct                                              >  1:407853/1‑76 (MQ=255)
                          tCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTctccgct                                              >  1:757479/1‑76 (MQ=255)
                              ttCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCtccgctccgc                                          >  1:238164/1‑76 (MQ=255)
                               tCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTccgctccgct                                         >  1:1294604/1‑76 (MQ=255)
                                ccTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCtg                                        >  2:1283338/1‑76 (MQ=255)
                                ccTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCtg                                        >  2:339588/1‑76 (MQ=255)
                                 cTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCtgn                                       <  1:1121396/76‑2 (MQ=255)
                                     tGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTgggt                                   >  2:935775/1‑74 (MQ=255)
                                       aaGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTa                                 >  2:809197/1‑76 (MQ=255)
                                        aGGGTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTAt                                >  1:358315/1‑76 (MQ=255)
                                          ggTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGt                              >  1:1401475/1‑76 (MQ=255)
                                          ggTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGt                              >  1:1200284/1‑76 (MQ=255)
                                           gTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGtt                             >  2:1331130/1‑76 (MQ=255)
                                           gTATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGtt                             >  1:1276070/1‑76 (MQ=255)
                                            tATGCACCCGGTTCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTc                            >  2:721016/1‑76 (MQ=255)
                                                       ttCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGctctc                 >  2:271915/1‑76 (MQ=255)
                                                        tCACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGctctct                >  2:105517/1‑76 (MQ=255)
                                                         cACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTa               <  1:200536/76‑1 (MQ=255)
                                                         cACCGTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTa               <  1:1136578/76‑1 (MQ=255)
                                                         cACCGTGCGGTGCTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGATCCTCGCGCTCTCTa               >  2:311672/1‑76 (MQ=255)
                                                             gTGCGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGc           <  1:570793/76‑1 (MQ=255)
                                                                cGGTGTTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCAtt        >  1:1323775/1‑76 (MQ=255)
                                                                    gtTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAg    >  1:1047565/1‑76 (MQ=255)
                                                                    gtTCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAg    <  2:9838/76‑1 (MQ=255)
                                                                     ttCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAgg   >  2:1300468/1‑76 (MQ=255)
                                                                      tCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGt  <  2:269208/76‑1 (MQ=255)
                                                                      tCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGt  >  2:1193062/1‑76 (MQ=255)
                                                                      tCCCGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGt  <  1:455853/76‑1 (MQ=255)
                                                                         |                                                                        
TGCCACGGTAACCTGATGAACTTCCGTCTGTTCCTGGTGAAGGGTATGCACCCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTCTCCGCTCCGCTGTGGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGT  >  REL606/1126902‑1127047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑