Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,658,770 T→G N200K (AAT→AAG ydgB → short chain dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,658,7700TG95.6% 170.1 / ‑4.3 46N200K (AAT→AAGydgBshort chain dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (26/17);  total (27/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 = 1.00e+00

ATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGC  >  REL606/1658700‑1658842
                                                                      |                                                                        
aTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATc                                                                     <  2:606750/76‑1 (MQ=255)
 ttCTATTGCGCCATNGCTGATCCTGTTTAATGAACATGATGATGCCGTATATCGACAACAGGCGCTTAATNACTCa                                                                    <  2:1343665/76‑1 (MQ=255)
 ttCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCa                                                                    >  1:620491/1‑76 (MQ=255)
 ttCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCa                                                                    >  1:121363/1‑76 (MQ=255)
 ttCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATACCGACAACAGGCGCTGAATAAATCa                                                                    >  1:695296/1‑76 (MQ=255)
  tCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCAc                                                                   <  1:256506/76‑1 (MQ=255)
    tATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACtg                                                                 >  2:1198224/1‑76 (MQ=255)
    tATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACtg                                                                 >  2:1198236/1‑76 (MQ=255)
     aTTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACtga                                                                >  2:285565/1‑76 (MQ=255)
     aTTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACtga                                                                <  1:385497/76‑1 (MQ=255)
      ttGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACtgat                                                               >  1:1219810/1‑76 (MQ=255)
       tGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACtgatg                                                              >  2:179161/1‑76 (MQ=255)
          gcCATCGCTGATCCNGTTTAATGAACATGATGATGCCGAATACCGACAACAGGCGCTGAACAACTCACTGATGaaa                                                           <  1:1144814/76‑1 (MQ=255)
          gcCATCGCTGATACTCTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGNAATCACTGATGaaa                                                           >  1:1151914/1‑76 (MQ=255)
              tCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACcg                                                       >  2:880726/1‑76 (MQ=255)
                gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACcgcg                                                     <  2:988241/76‑1 (MQ=255)
                gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACcgcg                                                     <  1:770030/76‑1 (MQ=255)
                 cTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCgcgc                                                    >  2:576015/1‑76 (MQ=255)
                   gATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCt                                                  <  1:830384/76‑1 (MQ=255)
                    aTCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTg                                                 >  1:1222162/1‑76 (MQ=255)
                     tCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTgg                                                >  1:756500/1‑76 (MQ=255)
                       cTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCg                                              <  2:103102/76‑1 (MQ=255)
                        tGTTTAATGAATATGATGATGCCGAATATCGACAACAGGTGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCga                                             >  1:36704/1‑76 (MQ=255)
                          tttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCgaga                                           >  2:579533/1‑76 (MQ=255)
                           ttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGagaa                                          <  1:701537/76‑1 (MQ=255)
                               tGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGagaaagaa                                      <  1:452629/76‑1 (MQ=255)
                               tGAACATGATGATGCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGagaaagaa                                      >  2:1102386/1‑76 (MQ=255)
                                      gatgatgCCGAATATCGACAACAGGCGCTGAAGGAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCg                               >  1:154917/1‑76 (MQ=255)
                                      gatgatgCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCg                               >  1:1084082/1‑76 (MQ=255)
                                        tgatgCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGAc                             >  2:779790/1‑76 (MQ=255)
                                            gCCGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTgg                         >  2:967570/1‑76 (MQ=255)
                                             ccGAATATCGACAACACGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGt                        <  1:398831/76‑1 (MQ=255)
                                              cGAATATCGACAACAGGCGCTGAAGAAATCCCTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTc                       <  2:617279/76‑1 (MQ=255)
                                              cGAATATCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTc                       >  2:808939/1‑76 (MQ=255)
                                                  tatCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAtt                   >  2:897677/1‑76 (MQ=255)
                                                  taccGACAACAGGCGCTGAAGAAATCACTGATGAAACCCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAtt                   <  1:435594/73‑1 (MQ=255)
                                                   atCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAtta                  >  2:1273472/1‑76 (MQ=255)
                                                    tCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttac                 >  1:550633/1‑76 (MQ=255)
                                                    tCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttac                 <  1:256638/76‑1 (MQ=255)
                                                    tCGACAACAGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttac                 <  2:1000501/76‑1 (MQ=255)
                                                        caacaGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttac             >  2:362527/1‑76 (MQ=255)
                                                         aacaGGCGCTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttacn            <  1:837668/76‑2 (MQ=255)
                                                            aGGCGCTGAAGAAATCACTTATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGAttacttacttac         <  2:395670/76‑1 (MQ=255)
                                                               cgcTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAg      >  1:357613/1‑76 (MQ=255)
                                                                 cTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTTGATGACTTACTTACCAGtt    <  2:117411/76‑1 (MQ=255)
                                                                 cTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGtt    <  2:178802/76‑1 (MQ=255)
                                                                 cTGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGtt    >  1:795619/1‑76 (MQ=255)
                                                                  tGAAGAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTg   <  1:520467/76‑1 (MQ=255)
                                                                   gaagaaATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGc  <  2:304628/76‑1 (MQ=255)
                                                                      |                                                                        
ATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAGAAGTGATCGACCTGGTCGATTACTTACTTACCAGTTGC  >  REL606/1658700‑1658842

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