Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,747,028 A→C 19.5% intergenic (+112/+97) srmB → / ← yfiE ATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,747,0280AC19.5% 154.9 / 37.0 77intergenic (+112/+97)srmB/yfiEATP‑dependent RNA helicase/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base C (7/8);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTA  >  Exported/1746959‑1747096
                                                                     |                                                                    
aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGa                                                                      >  1:3973708/1‑70 (MQ=255)
aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                                                     >  1:4283624/1‑71 (MQ=255)
aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                                                     >  1:1474693/1‑71 (MQ=255)
aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                                                     >  1:1613768/1‑71 (MQ=255)
 cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATc                                                                    <  1:947765/71‑1 (MQ=255)
 cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATc                                                                    <  1:5576898/71‑1 (MQ=255)
  gCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCt                                                                   <  1:3997055/71‑1 (MQ=255)
    aaaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                                                 <  1:5417556/71‑1 (MQ=255)
     aaatTGCCTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGt                                                                <  1:3822201‑M2/67‑1 (MQ=255)
     aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGt                                                                <  1:5319875/71‑1 (MQ=255)
     aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGt                                                                <  1:1747399/71‑1 (MQ=255)
     aaaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATAAATTGTATTTGCACGATCTTGt                                                                <  1:2638081/71‑1 (MQ=255)
      aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTg                                                                 >  1:1255450/1‑69 (MQ=255)
         tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGa                                                                      <  1:4109525/61‑1 (MQ=255)
         tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCAAGa                                                                      <  1:1051562/61‑1 (MQ=255)
           ccTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCg                                                          <  1:3962906/71‑1 (MQ=255)
           ccTGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCg                                                          <  1:337807/71‑1 (MQ=14)
            cTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCg                                                          >  1:3841755/1‑70 (MQ=255)
             tGATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTg                                                                 <  1:5583573/62‑1 (MQ=14)
              gATGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGa                                                        >  1:2652614/1‑70 (MQ=14)
               aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                                       <  1:4543787/70‑1 (MQ=255)
               aTGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGAt                                                       <  1:4336819/70‑1 (MQ=255)
               aTGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGAt                                                       <  1:2012413/70‑1 (MQ=14)
                tGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATa                                                      <  1:5656641/70‑1 (MQ=255)
                tGCGCTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATaa                                                     <  1:3724856/71‑1 (MQ=14)
                  cgcTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAg                                                    <  1:537265/70‑1 (MQ=255)
                  cgcTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTg                                                                 <  1:2101676/57‑1 (MQ=14)
                  cgcTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTg                                                                 <  1:1601060/57‑1 (MQ=14)
                   gcTACGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGc                                                  <  1:1839796/71‑1 (MQ=14)
                    cTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCg                                                 <  1:211605/71‑1 (MQ=255)
                    cTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCg                                                 <  1:3888979/71‑1 (MQ=255)
                      aCGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATaa                                                     >  1:4356994/1‑65 (MQ=14)
                       cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGc                                                            >  1:2093406/1‑57 (MQ=14)
                       cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGc                                                            >  1:3622613/1‑57 (MQ=14)
                       cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGtt                                               >  1:3952435/1‑70 (MQ=14)
                       cGCTTATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGtt                                               >  1:3799248/1‑70 (MQ=14)
                         cTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAc                                            <  1:116695/71‑1 (MQ=255)
                          ttATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGtt                                               <  1:1500793/67‑1 (MQ=255)
                           tATCAGGCCTACATGGCCCATGCAATGTATTGAATTTGCACGCTTTTGTAGGCCGGATAAGGCGTTTAcgc                                          >  1:2397133/1‑71 (MQ=11)
                              cAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgccg                                        >  1:2219218/1‑70 (MQ=255)
                               aGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCa                                      <  1:3075374/71‑1 (MQ=255)
                                  ccTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAgg                                                   <  1:877077/55‑1 (MQ=255)
                                  ccTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAgg                                                   <  1:5315591/55‑1 (MQ=255)
                                   cTACGTGGCTCATGCAATATATTGAATTTGCACGATc                                                                    >  1:2937576/1‑37 (MQ=255)
                                      cGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGc                                                  >  1:3182663/1‑52 (MQ=255)
                                       gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                                                     >  1:4465721/1‑48 (MQ=255)
                                       gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                                                     >  1:4243018/1‑48 (MQ=255)
                                       gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                     <  1:4993722/64‑1 (MQ=255)
                                       gTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                     <  1:3191174/64‑1 (MQ=255)
                                            tCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAt                          >  1:4811377/1‑70 (MQ=255)
                                            tCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAt                          >  1:2646128/1‑70 (MQ=255)
                                              aTGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATa                         >  1:952675/1‑69 (MQ=255)
                                               tGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg                                           <  1:2940977/50‑1 (MQ=255)
                                                gCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg                                           <  1:2895432/49‑1 (MQ=255)
                                                 cAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                                                     >  1:5574907/1‑38 (MQ=255)
                                                  aaTATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg                                           <  1:1655501/47‑1 (MQ=255)
                                                    tatatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                           <  1:2996414/61‑1 (MQ=255)
                                                     atatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                            >  1:3934492/1‑59 (MQ=255)
                                                      tatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGtt                                               <  1:4305556/39‑1 (MQ=40)
                                                      tatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAgcg                    >  1:3753945/1‑66 (MQ=255)
                                                      tatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAg                      >  1:2210189/1‑64 (MQ=255)
                                                      tatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACtt               >  1:2994516/1‑71 (MQ=255)
                                                      tatTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACtt               >  1:3855468/1‑71 (MQ=255)
                                                        ttGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTc                                              >  1:1032134/1‑38 (MQ=40)
                                                        ttGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCa                                             >  1:5161856/1‑39 (MQ=40)
                                                        ttGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCa                                             >  1:5242872/1‑39 (MQ=40)
                                                        ttGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCa                                             >  1:703815/1‑39 (MQ=40)
                                                        ttGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACttt              <  1:1318810/70‑1 (MQ=255)
                                                         tGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgc                                          >  1:5058831/1‑41 (MQ=38)
                                                         tGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgc                                          >  1:1350355/1‑41 (MQ=38)
                                                           aaTTTGCACGATCTTGTAGGCCGGATAAGGCGTTCa                                             >  1:4226580/1‑36 (MQ=40)
                                                           aaTTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTc           <  1:909282/70‑1 (MQ=255)
                                                             tttGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGAc        <  1:3447179/71‑1 (MQ=255)
                                                              ttGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGAc        >  1:1737715/1‑70 (MQ=255)
                                                               tGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACa       <  1:2079870/70‑1 (MQ=255)
                                                                  aCGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGtt   <  1:5129410/71‑1 (MQ=255)
                                                                   cGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACtt               >  1:3331635/1‑58 (MQ=255)
                                                                   cGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa  <  1:695727/71‑1 (MQ=255)
                                                                   cGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa  <  1:3360371/71‑1 (MQ=255)
                                                                    gATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTa  <  1:909579/70‑1 (MQ=255)
                                                                     |                                                                    
ACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAAAGCGCACTTTGTCGACAAGTTA  >  Exported/1746959‑1747096

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