Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,557,166 T→G 12.3% intergenic (+1016/+247) dsdA → / ← emrY D‑serine ammonia‑lyase/predicted multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,557,1660TG12.3% 133.1 / 18.4 65intergenic (+1016/+247)dsdA/emrYD‑serine ammonia‑lyase/predicted multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (22/35);  new base G (4/4);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

CCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAAT  >  Exported/1557101‑1557236
                                                                 |                                                                      
ccAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTc                                                                   <  1:916065/71‑1 (MQ=255)
 cAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTc                                                                   >  1:2715270/1‑70 (MQ=255)
   aTCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGGCAGTCCTg                                                                >  1:769471/1‑71 (MQ=255)
     ccAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCtt                                                                <  1:3033530/69‑2 (MQ=255)
     ccAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTg                                                                <  1:2646681/69‑1 (MQ=255)
      cAGATCTCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGt                                                              <  1:1676965/70‑1 (MQ=255)
      cAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGt                                                              <  1:988705/70‑1 (MQ=255)
      cAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGt                                                              <  1:3112796/70‑1 (MQ=255)
      cAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGt                                                              <  1:2542077/70‑1 (MQ=255)
      cAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGGCAGTCCTGGTa                                                             <  1:157496/71‑1 (MQ=255)
      cAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGGCAGTCCTGGTa                                                             <  1:2706676/71‑1 (MQ=255)
       aGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTg                                                                >  1:2262642/1‑67 (MQ=255)
       aGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTg                                                                >  1:1702838/1‑67 (MQ=255)
            gaagagaaactgagacttcaccggtcaaaatgaatctaccAGCACTAACAGGGTCAGTCCt                                                                 >  1:2520265‑M1/41‑61 (MQ=255)
            gaagagaaactgagacttcaccggtcaaaatgaatctaccAGCACTAACAGGGTCAGTCCTg                                                                >  1:2147982‑M1/41‑62 (MQ=255)
             cgcCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCt                                                                 <  1:2436118/60‑1 (MQ=255)
             cgcCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCt                                                                 <  1:2135999/60‑1 (MQ=255)
             cgcCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCt                                                                 <  1:1985070/60‑1 (MQ=255)
              agagaaactgagacttcaccggtcaaaatgaatctaccAGCACTAACAGGGTCAGTCCTGGt                                                              <  1:1297222‑M1/24‑1 (MQ=255)
               cccTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGt                                                              <  1:758487/61‑1 (MQ=255)
               cccTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGt                                                              <  1:339922/61‑1 (MQ=255)
                ccTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCAtttt                                                   >  1:1338735/1‑71 (MQ=255)
                ccTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCAtttt                                                   >  1:2588183/1‑71 (MQ=255)
                 cTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTg                                                  <  1:2225288/71‑1 (MQ=255)
                 cTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTg                                                  <  1:3230498/71‑1 (MQ=255)
                  tttATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGa                                                 >  1:514893/1‑71 (MQ=255)
                   ttATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGAc                                                >  1:391115/1‑71 (MQ=255)
                    tATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGa                                                 <  1:3406931/69‑1 (MQ=255)
                    tATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGa                                                 <  1:3365512/69‑1 (MQ=255)
                    tATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGa                                                 <  1:3188356/69‑1 (MQ=255)
                        aaGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGt                                            >  1:1141868/1‑70 (MQ=255)
                        aaGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTg                                           <  1:595572/71‑1 (MQ=255)
                         aGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGa                                          <  1:207863/71‑1 (MQ=255)
                          gCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGaa                                         <  1:493100/71‑1 (MQ=255)
                            aTAATTTGCAAGCCACCAACACCGAGCACTAACAGGGGCAGTCCTGGTAGATTCATTTTGACCGGTGAAGt                                       <  1:389127/71‑1 (MQ=255)
                              aaTTTGCAAGCCACCAACACCGAGCACTAACAGGGt                                                                        <  1:2499561/36‑1 (MQ=255)
                                tttGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGtct                                     >  1:3472611/1‑69 (MQ=255)
                                tttGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGtct                                     >  1:1788443/1‑69 (MQ=255)
                                tttGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGtct                                     >  1:1129244/1‑69 (MQ=255)
                                tttGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGtct                                     >  1:105533/1‑69 (MQ=255)
                                tttGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCa                                   <  1:1172205/71‑1 (MQ=255)
                                tttGCAAGCCACCAACACCGAGCACTAACAGGGGCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCa                                   <  1:3557672/71‑1 (MQ=255)
                                 ttGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCAtttt                                                   >  1:52568/1‑54 (MQ=255)
                                   gCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGtt                                >  1:1496330/1‑71 (MQ=255)
                                    cAAGCCACCAACACCGAGCACTAACAGGGGCAGTCCTGGTAGATTc                                                        >  1:153614/1‑46 (MQ=255)
                                    cAAGCCACCAACACCGAGCACTAACAGGGGCAGTCCTGGTAGATTc                                                        >  1:2249074/1‑46 (MQ=255)
                                    cAAGCCACCAACACCGAGCACTAACAGGGGCAGTCCTGGTAGATTc                                                        >  1:2302518/1‑46 (MQ=255)
                                      aGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTtct                             <  1:21760/71‑1 (MQ=255)
                                       gCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGt                                 <  1:874991/66‑1 (MQ=255)
                                       gCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGt                                 <  1:3071209/66‑1 (MQ=255)
                                       gCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTtctc                            >  1:1529184/1‑71 (MQ=255)
                                       gCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTtctc                            >  1:2870503/1‑71 (MQ=255)
                                         caccaACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCtt                          <  1:629004/71‑1 (MQ=255)
                                          accaACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTctct                           <  1:3402445/69‑1 (MQ=255)
                                          accaACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTc                         >  1:2874902/1‑71 (MQ=255)
                                               cacCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTta                    <  1:3231655/71‑1 (MQ=255)
                                                   gAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTtaagt                 >  1:1503987/1‑70 (MQ=255)
                                                    AGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTtaagt                 <  1:672222‑M1/69‑1 (MQ=255)
                                                        cTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGtt           <  1:3233846/71‑1 (MQ=255)
                                                        cTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGt            >  1:2067618/1‑70 (MQ=255)
                                                         tAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGtt                                <  1:1270303/49‑1 (MQ=255)
                                                         tAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGtt                                <  1:3166945/49‑1 (MQ=255)
                                                           aCAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTaa         <  1:2389424/70‑1 (MQ=255)
                                                           aCAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTaa         <  1:2347038/70‑1 (MQ=255)
                                                             aGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCa      >  1:194623/1‑71 (MQ=255)
                                                              ggggCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCAt     <  1:932618/71‑1 (MQ=255)
                                                              gggTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGtt           <  1:1811833/65‑1 (MQ=255)
                                                                 tCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAAt  <  1:2757572/71‑1 (MQ=255)
                                                                 |                                                                      
CCAATCCAGATCGCGCCCTTTATCAAGCATAATTTGCAAGCCACCAACACCGAGCACTAACAGGGTCAGTCCTGGTAGATTCATTTTGACCGGTGAAGTCTCAGTTTCTCTTCCTTTAAGTAAGGTTAAGCATAAT  >  Exported/1557101‑1557236

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