Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 474,455 A→G L284P (CTG→CCG)  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,4550AG95.8% 167.4 / ‑2.8 48L284P (CTG→CCG) ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (31/15);  total (31/17)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTCAG  >  REL606/474383‑474525
                                                                        |                                                                      
gTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGc                                                                      >  2:528313/1‑75 (MQ=255)
 tCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCa                                                                     <  2:1093633/75‑1 (MQ=255)
     cAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGc                                                                 >  1:1590312/1‑75 (MQ=255)
       gCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCaa                                                               <  2:1043431/75‑1 (MQ=255)
       gCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCaa                                                               >  1:1507796/1‑75 (MQ=255)
        cACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCaaa                                                              <  2:315468/75‑1 (MQ=255)
        cACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCaaa                                                              <  1:425842/75‑1 (MQ=255)
          cTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAAcg                                                            >  2:740159/1‑75 (MQ=255)
            gCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAAcgcg                                                          <  2:529902/75‑1 (MQ=255)
             ccAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGt                                                         >  1:1526503/1‑75 (MQ=255)
               aGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTcg                                                       >  2:1449439/1‑75 (MQ=255)
                   gTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCa                                                   <  1:894417/75‑1 (MQ=255)
                     tGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAAt                                                 >  1:771673/1‑75 (MQ=255)
                      gCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATg                                                <  2:415966/75‑1 (MQ=255)
                         gAATCAGAATTAACGGATCAAACAACAGCCCGACGGAGACAAAAAACGGCGCCGCAAACGCGTCGCGCAATGGCa                                             <  2:317142/75‑1 (MQ=255)
                          aaTCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAg                                            >  1:120344/1‑75 (MQ=255)
                           aTCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGc                                           <  2:1150584/75‑1 (MQ=255)
                           aTCAGAATTAACCGATAAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGc                                           <  2:726225/75‑1 (MQ=255)
                              aGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTa                                        >  1:857304/1‑75 (MQ=255)
                              aGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTa                                        >  2:361319/1‑75 (MQ=255)
                               gAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTAt                                       >  1:1448442/1‑75 (MQ=255)
                                   tAACGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTg                                   >  1:56206/1‑75 (MQ=255)
                                    aaCGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTgg                                  >  1:1424875/1‑75 (MQ=255)
                                     aCGGATCAAACTACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTggg                                 >  2:497203/1‑75 (MQ=255)
                                     aCGGATCAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTggg                                 >  1:445325/1‑75 (MQ=255)
                                           cAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGggcggcgg                           >  1:37574/1‑73 (MQ=255)
                                           cAAACAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCAc                           <  2:64943/75‑1 (MQ=255)
                                            aaaCAACATCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACg                          >  1:1450402/1‑75 (MQ=255)
                                             aacaacATCCCGACGGAGACAAAACACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACgg                         >  2:803940/1‑75 (MQ=255)
                                                  caTCCCGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCCCAGGGACt                    >  1:1603160/1‑75 (MQ=255)
                                                   aTCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTc                   <  1:217050/75‑1 (MQ=255)
                                                     cccGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAg                 <  1:672435/75‑1 (MQ=255)
                                                     cccGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAg                 <  2:1228629/75‑1 (MQ=255)
                                                      ccGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGt                <  1:765589/75‑1 (MQ=255)
                                                      ccGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGt                <  2:1058731/75‑1 (MQ=255)
                                                      ccGACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGGATCGTGGGCGGCACGGGGACTCAgg                >  1:399389/1‑74 (MQ=255)
                                                       cGACGGAGACAAAAAACGGTCCCGCAAACGCGGCGCGCAAAGGCTGCGTATCGTGGGCGGCACGGTGACTCAGtt               <  2:826698/75‑1 (MQ=255)
                                                        gACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTc              >  2:270291/1‑75 (MQ=255)
                                                        gACGGAGACAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGCGACTCAGTTc              >  1:1584496/1‑75 (MQ=255)
                                                             agaCAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGTGCGGCACGGGGACTCAGTTCAGACt         >  1:771552/1‑75 (MQ=255)
                                                             agaCAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACt         >  1:1481054/1‑75 (MQ=255)
                                                              gaCAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTc        >  1:676411/1‑75 (MQ=255)
                                                               aCAAAAAACGGCTCCGCAAGCGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCg       <  1:1535035/75‑1 (MQ=255)
                                                                cAAAAAACGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGt      >  1:1092714/1‑75 (MQ=255)
                                                                 aaaaaaCAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGtt     <  1:59914/75‑1 (MQ=255)
                                                                  aaaaaCGGCGCCGCAACCGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGGCTCAGTTCAGACTCGTTc    <  2:643600/75‑1 (MQ=255)
                                                                  aaaaaCGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTc    <  1:733834/75‑1 (MQ=255)
                                                                  aaaaaCGGCACCGCAAACGCGGCGCGCAATGGCAGCGTATCGGGGGCGGCACGGTGACTCAGTTCAGACTCGTTc    >  1:955728/1‑75 (MQ=255)
                                                                   aaaaCGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTCa   >  1:1598968/1‑75 (MQ=255)
                                                                   aaaaCGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTCa   >  1:715214/1‑75 (MQ=255)
                                                                   aaaaCGGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTCa   >  1:1246688/1‑75 (MQ=255)
                                                                    aaaCGGCACCGCAAGCGCGTCGAGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTCAg  <  1:1392127/75‑1 (MQ=255)
                                                                    aaaCGGCACCGCAAACGCGTCGCGCAATGGCAGGGTATCGTGGGCGGCACGGGGACTCAGTTCAGAGTCGTTCAg  >  2:1607444/1‑75 (MQ=255)
                                                                        |                                                                      
GTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCCGACGGAGACAAAAAACAGCACCGCAAACGCGTCGCGCAATGGCAGCGTATCGTGGGCGGCACGGTGACTCAGTTCAGACTCGTTCAG  >  REL606/474383‑474525

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