Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,367,537 C→T 42.1% intergenic (‑85/‑105) pepQ ← / → fadB proline dipeptidase/fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,367,5370CT42.1% 26.8 / 68.5 57intergenic (‑85/‑105)pepQ/fadBproline dipeptidase/fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase and 3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (14/19);  new base T (16/8);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

TTTGACCTTAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTGG  >  Exported/2367471‑2367605
                                                                  |                                                                    
tttGACCTTAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAg                                                                    <  1:1630908/69‑1 (MQ=255)
        tAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCAttt                                                          >  1:1509256/1‑71 (MQ=255)
          aGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCt                                                        <  1:2869113/71‑1 (MQ=255)
           gTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCtt                                                       <  1:1934330/71‑1 (MQ=255)
            tctcCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCttt                                                      <  1:3408843/71‑1 (MQ=255)
                cGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAAt                                                   <  1:3147531/70‑1 (MQ=255)
                 gCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCt                                                 >  1:1153358/1‑71 (MQ=255)
                      ttAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGTCATTTCTTTAATCtttt                                              >  1:1787656/1‑69 (MQ=255)
                      ttAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCtttt                                              >  1:411875/1‑69 (MQ=255)
                      ttAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCt                                                 >  1:2158012/1‑66 (MQ=255)
                      ttAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCt                                                 >  1:720704/1‑66 (MQ=255)
                        aGCACATCGTTCATCCAGAGGGTGATTTCTGCCGAGCGTGATCAGATCGGCAt                                                            <  1:1555531/53‑1 (MQ=255)
                        aGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCAt                                                            <  1:3261401/53‑1 (MQ=255)
                        aGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCAt                                                            <  1:3448063/53‑1 (MQ=255)
                        aGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTtttgt                                            <  1:3091992/69‑1 (MQ=255)
                        aGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTtttgt                                            <  1:1661034/69‑1 (MQ=255)
                            caTCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCttt                                                      <  1:2660/55‑1 (MQ=255)
                             aTCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCg                                                                <  1:2381800/44‑1 (MQ=255)
                             aTCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCg                                                                <  1:2259267/44‑1 (MQ=255)
                             aTCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCg                                                                <  1:458920/44‑1 (MQ=255)
                              tCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTtttgttt                                          >  1:2305522/1‑65 (MQ=255)
                              tCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTtttgttt                                          >  1:1451512/1‑65 (MQ=255)
                                 ttCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATAtttt                                 <  1:1986811/71‑1 (MQ=255)
                                  tCATCCAGAGCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATAtttt                                 >  1:1464828/1‑70 (MQ=255)
                                  tCATCCAGAGAGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCat                                      >  1:1454836/1‑65 (MQ=255)
                                         gagCGTGATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACaa                         <  1:2069181/71‑1 (MQ=255)
                                            cGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACaaaa                       >  1:2980576/1‑70 (MQ=255)
                                             gTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATa                     <  1:1907395/71‑1 (MQ=255)
                                             gTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATa                     <  1:1698230/71‑1 (MQ=255)
                                             gTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATa                     <  1:1682041/71‑1 (MQ=255)
                                               gATTTCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATaca                   >  1:2407305/1‑71 (MQ=255)
                                               gATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACaaa                        >  1:2685032/1‑66 (MQ=255)
                                                aTTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCtt                                                >  1:262226/1‑41 (MQ=255)
                                                  ttCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCat                                      >  1:87565/1‑49 (MQ=255)
                                                   tCTGCCGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATacagac                <  1:3040400/70‑1 (MQ=255)
                                                       ccGAGCGTGATTAGATCGGCATTTCTTTAATCTtttgtttg                                         >  1:183865/1‑41 (MQ=255)
                                                       ccGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACtt              <  1:998389/68‑1 (MQ=255)
                                                       ccGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGa           >  1:2236831/1‑71 (MQ=255)
                                                        cGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTaa                              >  1:949488/1‑51 (MQ=255)
                                                        cGAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTaa                              >  1:1322392/1‑51 (MQ=255)
                                                         gAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGAc          >  1:3113855/1‑70 (MQ=255)
                                                         gAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACt         >  1:3430043/1‑71 (MQ=255)
                                                         gAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACt         >  1:961128/1‑71 (MQ=255)
                                                         gAGCGTGATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACt         >  1:496643/1‑71 (MQ=255)
                                                           gCGTGATCAGATCGGCATTTCTTTAATCTtttgtttg                                         >  1:3599461/1‑37 (MQ=255)
                                                           gCGTGATCAGATCGGCATTTCTTTAATCTtttgtttg                                         >  1:1224096/1‑37 (MQ=255)
                                                           gCGTGATCAGATCGGCATTTCTTTAATCTtttgttt                                          >  1:164664/1‑36 (MQ=255)
                                                           gCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCa       <  1:3311538/71‑1 (MQ=255)
                                                           gCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCa       <  1:3307186/71‑1 (MQ=255)
                                                               gATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTg   >  1:3335343/1‑71 (MQ=255)
                                                               gATTAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTg   >  1:2712308/1‑71 (MQ=255)
                                                                aTCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTgg  >  1:2874440/1‑71 (MQ=255)
                                                                aTCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTgg  >  1:1073861/1‑71 (MQ=255)
                                                                aTCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTg   <  1:3682802/70‑1 (MQ=255)
                                                                aTCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTg   <  1:1199817/70‑1 (MQ=255)
                                                                aTCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTg   <  1:2197024/70‑1 (MQ=255)
                                                                aTCAGATCGGCATTTCTTTAATCTTATGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTg   <  1:3277797/70‑1 (MQ=255)
                                                                  |                                                                    
TTTGACCTTAAGTCTCCGCATCTTAGCACATCGTTCATCCAGAGCGTGATTTCTGCCGAGCGTGATCAGATCGGCATTTCTTTAATCTTTTGTTTGCATATTTTTAACACAAAATACACACTTCGACTCATCTGG  >  Exported/2367471‑2367605

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: