Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 474,424 A→T D294E (GAT→GAA ybaL ← predicted transporter with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606474,4240AT97.9% 166.3 / NA 48D294E (GAT→GAAybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (26/21);  minor base C (1/0);  total (27/21)
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

AACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCG  >  REL606/474353‑474496
                                                                       |                                                                        
aaCGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCaa                                                                      >  2:1860398/1‑76 (MQ=255)
aaCGACATCCCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCaa                                                                      >  1:1899020/1‑76 (MQ=255)
 aCGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCaaa                                                                     >  2:707294/1‑76 (MQ=255)
  cGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAaac                                                                    <  1:1104319/76‑1 (MQ=255)
  cGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAaac                                                                    <  1:1104338/76‑1 (MQ=255)
        aCCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATc                                                              >  1:367638/1‑76 (MQ=255)
          cAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATccc                                                            <  1:1217863/76‑1 (MQ=255)
           aaaCAGANAAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCg                                                           >  2:1642880/1‑76 (MQ=255)
           aaaCAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCg                                                           <  2:588893/76‑1 (MQ=255)
           aaaCAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCg                                                           >  2:1319189/1‑76 (MQ=255)
           aaaCAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATGAACGGGTCAAACAACATCCCg                                                           >  2:1223804/1‑76 (MQ=255)
             aCAGAATAAGCGCCAGCGTCGCCAGCACTGCCAGCAGTTCCTGAATCAGAATTAACGGTTCAAACAACATCCCGAc                                                         <  2:1633054/76‑1 (MQ=255)
             aCAAAATGAGCGCCAGCGGCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGAc                                                         <  2:821121/76‑1 (MQ=255)
               aGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACgg                                                       >  1:1413048/1‑76 (MQ=255)
                gAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGga                                                      <  1:759964/76‑1 (MQ=255)
                 aataatCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGgag                                                     <  1:870839/76‑1 (MQ=255)
                     atCGCCGGCATCGCCAGCGCTGCCAGCGGTTGCTGAATCAGAATTAACGATTCAAACAACANCCNGACGGAGACaa                                                 >  2:311702/1‑76 (MQ=255)
                      tCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACaaa                                                >  2:1464511/1‑76 (MQ=255)
                       cGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACaaaa                                               <  1:694511/76‑1 (MQ=255)
                          cAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAAc                                            <  1:407659/76‑1 (MQ=255)
                          cAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAAc                                            <  1:1395792/76‑1 (MQ=255)
                          cAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAAc                                            <  2:1131335/76‑1 (MQ=255)
                          cAGCGTCGCCAGCACTGCCAGCGGTGGCCGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAAc                                            <  1:319440/76‑1 (MQ=255)
                           aGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCTAAGAACATCCCGACGGAGACAAAAAACa                                           >  1:314542/1‑76 (MQ=255)
                           aGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACa                                           >  1:620407/1‑76 (MQ=255)
                              gTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCa                                        >  1:860524/1‑76 (MQ=255)
                                cGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCAcc                                      >  2:1562590/1‑76 (MQ=255)
                                   cAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCa                                   >  1:655080/1‑76 (MQ=255)
                                    aGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCaa                                  >  2:1087471/1‑76 (MQ=255)
                                      cACTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAAc                                <  1:767095/76‑1 (MQ=255)
                                        cTGCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGAACCGCAAAcgc                              >  1:230326/1‑76 (MQ=255)
                                          gCCAGCGGTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGt                            >  2:1356451/1‑76 (MQ=255)
                                                ggTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCa                      <  2:1437727/76‑1 (MQ=255)
                                                ggTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCa                      >  1:367260/1‑76 (MQ=255)
                                                 gTTGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCaa                     >  1:887864/1‑76 (MQ=255)
                                                   tGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATg                   >  1:1044701/1‑76 (MQ=255)
                                                   tGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATg                   <  1:1703842/76‑1 (MQ=255)
                                                   tGCTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATg                   >  1:719159/1‑76 (MQ=255)
                                                     cTGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGc                 <  2:1044700/76‑1 (MQ=255)
                                                      tGAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCa                <  1:1228712/76‑1 (MQ=255)
                                                       gAATCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAg               >  1:1556567/1‑76 (MQ=255)
                                                        aaTCAGAATTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGc              <  1:1857623/76‑1 (MQ=255)
                                                            aGAATTAACGGNTCNCACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGCCGCGCAATCGCAGCGTAt          <  2:314515/76‑1 (MQ=255)
                                                               aTTAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGt       <  1:604073/76‑1 (MQ=255)
                                                                ttAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCACACCCGTCGCGCAATGGCAGCGTATCGTg      >  1:321865/1‑76 (MQ=255)
                                                                ttAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTg      >  1:1336733/1‑76 (MQ=255)
                                                                 tAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTgg     >  1:1538414/1‑76 (MQ=255)
                                                                 tAACGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTgg     <  2:632803/76‑1 (MQ=255)
                                                                 tAACGGCTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTgg     >  1:1064185/1‑76 (MQ=255)
                                                                    cGGTTCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGggcg  >  2:981980/1‑76 (MQ=255)
                                                                       |                                                                        
AACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCG  >  REL606/474353‑474496

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