Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 783,391 G→A G292G (GGG→GGA ybhI → predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606783,3910GA97.7% 159.8 / ‑3.8 43G292G (GGG→GGAybhIpredicted transporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (24/18);  total (25/18)
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.95e-01

GCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTACCTGGGAAGAC  >  REL606/783324‑783463
                                                                   |                                                                        
gcgAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTcgatga                                                                  >  2:1611296/1‑76 (MQ=255)
    aaaaaTGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATcc                                                              <  2:1373628/76‑1 (MQ=255)
    aaaaaTGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATcc                                                              <  1:42291/76‑1 (MQ=255)
    aaaaaTGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATcc                                                              <  1:356924/76‑1 (MQ=255)
       aaTGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCAcc                                                           <  1:1621560/76‑1 (MQ=255)
         tgctgcTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGt                                                         <  1:333312/76‑1 (MQ=255)
           ctgctgGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTgg                                                       >  2:1095723/1‑76 (MQ=255)
            tgctCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGc                                                      >  2:1421577/1‑76 (MQ=255)
             gctCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCa                                                     >  2:99919/1‑76 (MQ=255)
                cGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATc                                                  <  2:1496735/76‑1 (MQ=255)
                     tCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCgttgt                                             >  1:11422/1‑76 (MQ=255)
                      cTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCgttgtt                                            <  1:96611/76‑1 (MQ=255)
                         tttgCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATg                                         >  2:1119615/3‑76 (MQ=255)
                              tGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAAc                                    >  2:206756/1‑76 (MQ=255)
                                gcgcTGCTGGGCTGGATTTTCAGTAAGTCTCTGAGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCa                                  <  1:922521/76‑1 (MQ=255)
                                  gctgctGGGCTGGATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCAtg                                <  1:1538093/76‑1 (MQ=255)
                                      ctgggctggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCAtgctgc                            <  1:1057932/76‑1 (MQ=255)
                                      ctgggctggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCAtgctgc                            >  1:1355159/1‑76 (MQ=255)
                                      ctgggctggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCAggctgc                            >  2:355275/1‑76 (MQ=255)
                                       tgggctggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCAtgctgct                           >  1:768349/1‑76 (MQ=255)
                                       tgggctggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCAcgctgct                           <  1:697060/76‑1 (MQ=255)
                                       tgggctggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACAAtgctgct                           >  2:323063/1‑76 (MQ=255)
                                        gggctggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCAtgctgctg                          >  1:12885/1‑76 (MQ=255)
                                           ctggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGctgctgctg                       <  2:669560/76‑1 (MQ=255)
                                             ggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTggg                     >  2:330535/1‑76 (MQ=255)
                                             ggATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTggg                     >  1:815347/1‑76 (MQ=255)
                                              gATTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGt                    >  2:92617/1‑76 (MQ=255)
                                               gttttCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTa                   >  1:855400/2‑76 (MQ=255)
                                               aTTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTa                   >  2:821493/1‑76 (MQ=255)
                                               aTTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTa                   >  1:35724/1‑76 (MQ=255)
                                               aTTTTCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTa                   >  1:1567850/1‑76 (MQ=255)
                                                ttttCAGTAAGTCTCTGGGGGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTAt                  >  1:1271180/1‑76 (MQ=255)
                                                ttttCAGTAAGTCTCTGGGAGTTGATGAATCCCCCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTAt                  <  2:1206832/76‑1 (MQ=255)
                                                  ttCAGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCg                <  1:799437/76‑1 (MQ=255)
                                                     aGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTa             >  1:124954/1‑76 (MQ=255)
                                                     aGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTa             >  2:1159629/1‑76 (MQ=255)
                                                     aGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTa             <  1:925698/76‑1 (MQ=255)
                                                     aGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTa             <  2:996533/76‑1 (MQ=255)
                                                     aGTAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTCATGGCAACCATGCTGCTGCAGACTATCGCTa             >  1:428235/1‑76 (MQ=255)
                                                       tAAGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTAcc           >  1:260965/1‑76 (MQ=255)
                                                        aaGTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTACCt          <  2:73778/76‑1 (MQ=255)
                                                          gTCTCTGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTACCTgg        <  1:361385/76‑1 (MQ=255)
                                                            ctctGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTACCTGGGa      <  2:388408/76‑1 (MQ=255)
                                                            ctctGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTACCTGGGa      >  2:891968/1‑76 (MQ=255)
                                                            ctctGGGAGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGGATCGTTACCTGGGa      >  2:891969/1‑76 (MQ=255)
                                                                gggagTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTACCTGGGAAGAc  <  1:520999/76‑1 (MQ=255)
                                                                   |                                                                        
GCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGGGTTGATGAATCCACCGTGGCAATCGTTGTTATGGCAACCATGCTGCTGCTGGGTATCGTTACCTGGGAAGAC  >  REL606/783324‑783463

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