Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 500,385 T→G D90E (GAT→GAG ylbH → conserved protein, rhs‑like

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606500,3850TG97.4% 809.3 / 5.2 189D90E (GAT→GAGylbHconserved protein, rhs‑like
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base G (78/106);  total (82/107)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 191 total aligned reads displayed.

CTCAGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGG  >  REL606/500312‑500457
                                                                         |                                                                        
cTCAGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGcc                                                                        >  2:3622925/1‑76 (MQ=255)
  cAGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCAt                                                                      >  1:2638101/1‑76 (MQ=255)
  cAGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCAt                                                                      >  2:3828531/1‑76 (MQ=255)
  cAGGACCAGATAGGGTTGGAGGGGGGAGGGAGTCTGTAAGCGTATTCGCTGAATCCGGTGGAAGGTATAGATCCAt                                                                      >  2:893606/1‑76 (MQ=255)
   aGGACCGGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCAtt                                                                     <  2:3332459/76‑1 (MQ=255)
   aGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCAtt                                                                     >  2:4228522/1‑76 (MQ=255)
    ggccccGATAGGGCGGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTa                                                                    <  1:2762435/73‑1 (MQ=37)
    ggACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTa                                                                    <  1:865964/76‑1 (MQ=255)
    ggACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTa                                                                    <  1:4789642/76‑1 (MQ=255)
    ggACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTAAGCCCTTa                                                                    >  2:1299476/1‑76 (MQ=37)
      aCCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAgg                                                                  <  2:5093514/76‑1 (MQ=255)
         cGATAGGGCTGGGGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGtt                                                               <  2:605437/76‑1 (MQ=255)
         cGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGtt                                                               >  2:2471101/1‑76 (MQ=255)
          gATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTa                                                              <  2:3270004/76‑1 (MQ=255)
          gATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTa                                                              >  1:4216700/1‑76 (MQ=255)
           aTAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTaa                                                             <  2:4850763/76‑1 (MQ=255)
             aGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGt                                                           >  2:3406825/1‑76 (MQ=255)
             aGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGt                                                           >  1:4996484/1‑76 (MQ=255)
                 cTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCg                                                       <  1:2992776/76‑1 (MQ=255)
                  tGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGc                                                      <  2:1151486/76‑1 (MQ=255)
                   ggAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCa                                                     <  1:1910602/76‑1 (MQ=255)
                    gAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAg                                                    <  2:2015274/76‑1 (MQ=255)
                     aGGGGGGCTGGAGTCTGTATGCGTATCCGCTGAATCCGGTGACTGGTATTGAGCCATTAGGGTTAAGTCACGCAGa                                                   >  2:4301758/1‑76 (MQ=38)
                     aGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGa                                                   >  2:2356378/1‑76 (MQ=255)
                     aGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGa                                                   >  1:3837815/1‑76 (MQ=255)
                      ggggggATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGAt                                                  <  1:2227422/76‑1 (MQ=255)
                      ggggggATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGAt                                                  <  2:3941398/76‑1 (MQ=255)
                         gggATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGCGCCATTAGGGTTAAGTCCCGCAGATGTa                                               >  2:1282355/1‑76 (MQ=255)
                         gggATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTa                                               >  2:2709302/1‑76 (MQ=255)
                          ggATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAg                                              >  1:2046532/1‑76 (MQ=255)
                          ggATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAg                                              <  1:3828531/76‑1 (MQ=255)
                            aTGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAgcg                                            >  2:1255314/1‑76 (MQ=255)
                            aTGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAgcg                                            <  1:1089214/76‑1 (MQ=255)
                            aTGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAgcg                                            >  1:331334/1‑76 (MQ=255)
                            aTGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAgcg                                            <  1:1109013/76‑1 (MQ=255)
                            aTGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAgcg                                            >  1:1166711/1‑76 (MQ=255)
                             tGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAgcgc                                           >  2:2110791/1‑76 (MQ=255)
                              ggAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCt                                          >  2:4508014/1‑76 (MQ=255)
                               gAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCta                                         >  2:1411490/1‑76 (MQ=255)
                               gAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCta                                         <  1:4077979/76‑1 (MQ=255)
                               gAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCta                                         <  1:3619512/76‑1 (MQ=255)
                                aGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCtaa                                        <  1:2292894/76‑1 (MQ=255)
                                 gTCTGTATGCTTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCtaat                                       <  1:46011/76‑1 (MQ=255)
                                 gTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCtaat                                       <  2:1846186/76‑1 (MQ=255)
                                  tCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCtaata                                      <  1:1307103/76‑1 (MQ=255)
                                  tCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCtaata                                      <  2:4547094/76‑1 (MQ=255)
                                   cTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCtaataa                                     >  2:995601/1‑76 (MQ=255)
                                    tgtgtgCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaag                                    <  2:4154586/76‑1 (MQ=255)
                                    tGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaag                                    >  2:1978439/1‑76 (MQ=255)
                                    tGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaag                                    <  2:3490770/76‑1 (MQ=255)
                                    tGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaag                                    <  2:4347617/76‑1 (MQ=255)
                                     gTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaaga                                   >  1:471464/1‑76 (MQ=255)
                                     gTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaaga                                   >  1:4587889/1‑76 (MQ=255)
                                     gTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaaga                                   >  2:3561913/1‑76 (MQ=255)
                                     gTATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaaga                                   >  2:2826100/1‑76 (MQ=255)
                                      tATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaagaa                                  <  1:2382354/76‑1 (MQ=255)
                                      tATGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaagaa                                  <  1:4970773/76‑1 (MQ=255)
                                       aTGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaagaag                                 >  1:5013790/1‑76 (MQ=255)
                                       aTGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaagaag                                 >  2:3310032/1‑76 (MQ=255)
                                       aTGCGTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATaagaag                                 <  1:381121/76‑1 (MQ=255)
                                           gTATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGaaaa                             <  2:3150692/76‑1 (MQ=255)
                                            tATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAg                            <  1:3069779/76‑1 (MQ=255)
                                            tATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAg                            >  1:963644/1‑76 (MQ=255)
                                            tATCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAg                            <  2:1750228/76‑1 (MQ=255)
                                             aTCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGa                           <  1:257140/76‑1 (MQ=255)
                                             aTCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGa                           <  2:5117130/76‑1 (MQ=255)
                                              tCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGAt                          <  1:1242569/76‑1 (MQ=255)
                                              tCCGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGAt                          <  1:801288/76‑1 (MQ=255)
                                                cGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCa                        >  2:4961328/1‑76 (MQ=255)
                                                cGCTGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCa                        <  2:2638101/76‑1 (MQ=255)
                                                   tGAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACt                     <  1:3689216/76‑1 (MQ=255)
                                                    gAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTa                    >  2:2396433/1‑76 (MQ=255)
                                                    gAATCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCGCCGGTGGAGCGGTTATTAGGAGAAAAGATCCACTaa                   >  1:3270672/1‑76 (MQ=255)
                                                     aaTCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTaa                   >  2:1602878/1‑76 (MQ=255)
                                                     aaTCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTaa                   <  1:5108162/76‑1 (MQ=255)
                                                     aaTCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTaa                   <  1:3495735/76‑1 (MQ=255)
                                                      aTCCGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTaaa                  <  2:2366748/76‑1 (MQ=255)
                                                        ccGGTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAAcc                <  1:1006698/76‑1 (MQ=255)
                                                          ggTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCAt              <  1:710247/76‑1 (MQ=255)
                                                           gTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATc             <  2:4996484/76‑1 (MQ=255)
                                                           gTGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATc             <  1:3931322/76‑1 (MQ=255)
                                                            tGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCa            <  1:2303785/76‑1 (MQ=255)
                                                            tGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCa            <  1:435728/76‑1 (MQ=255)
                                                            tGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCa            <  2:708569/76‑1 (MQ=255)
                                                            tGAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCa            <  1:2119734/76‑1 (MQ=255)
                                                             gAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCaa           >  1:1340351/1‑76 (MQ=255)
                                                             gAATGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCaa           <  2:1040020/76‑1 (MQ=255)
                                                              aaTGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCaaa          >  2:589031/1‑76 (MQ=255)
                                                               aTGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAag         >  2:2903601/1‑76 (MQ=255)
                                                               aTGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAag         <  1:1188166/76‑1 (MQ=255)
                                                               aTGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAag         <  1:2730387/76‑1 (MQ=255)
                                                               aTGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAag         <  2:2046532/76‑1 (MQ=255)
                                                                tGGTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAaga        <  2:3130477/76‑1 (MQ=255)
                                                                 ggTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAagag       <  1:2870234/76‑1 (MQ=255)
                                                                 ggTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAagag       <  1:4480815/76‑1 (MQ=255)
                                                                 ggTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAagag       <  1:2452230/76‑1 (MQ=255)
                                                                 ggTATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAagag       <  1:2390984/76‑1 (MQ=255)
                                                                  gTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGc      <  1:3388140/76‑1 (MQ=255)
                                                                   tATTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCa     <  1:2242346/76‑1 (MQ=255)
                                                                    aTTGAGCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCAt    >  1:2365027/1‑76 (MQ=255)
                                                                         |                                                                        
CTCAGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGG  >  REL606/500312‑500457

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