Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,739,943 T→C 46.4% intergenic (‑133/+65) lepA ← / ← rseC GTP binding membrane protein/RseC protein involved in reduction of the SoxR iron‑sulfur cluster

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,739,9430TC46.4% 8.9 / 97.1 69intergenic (‑133/+65)lepA/rseCGTP binding membrane protein/RseC protein involved in reduction of the SoxR iron‑sulfur cluster
Reads supporting (aligned to +/- strand):  ref base T (20/17);  new base C (19/13);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

ACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCA  >  Exported/1739877‑1740010
                                                                  |                                                                   
aCGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATg                                                                  >  1:2456299/1‑70 (MQ=255)
  gTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGAt                                                               >  1:2043834/1‑71 (MQ=255)
  gTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGAt                                                               >  1:2506525/1‑71 (MQ=255)
    cTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATtg                                                             >  1:682010/1‑71 (MQ=255)
       tATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATtg                                                             <  1:2166457/68‑1 (MQ=255)
       tATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATtg                                                             <  1:2650341/68‑1 (MQ=255)
           aaTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATc                                                      >  1:1347626/1‑71 (MQ=255)
            aTAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTAt                                                       <  1:2355933/69‑1 (MQ=255)
             tAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTg                                                    >  1:2271136/1‑71 (MQ=255)
             tAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTg                                                    >  1:1214748/1‑71 (MQ=255)
             tAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTg                                                    >  1:2620524/1‑71 (MQ=255)
               gAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATtgtg                                                           >  1:1700992/1‑62 (MQ=255)
                aaaCGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGaa                                                 <  1:2417980/71‑1 (MQ=255)
                 aaCGCCGCATTTTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGaaa                                                <  1:2077811/71‑1 (MQ=255)
                      cgcATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGaaa                                                <  1:159873/66‑1 (MQ=255)
                        cATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTc                                         >  1:1643407/1‑71 (MQ=255)
                        cATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTc                                         >  1:2246855/1‑71 (MQ=255)
                         aTTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCg                                        >  1:499903/1‑71 (MQ=255)
                          ttCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCg                                        >  1:843347/1‑70 (MQ=255)
                            ctacactacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAt                                     >  1:608872/1‑71 (MQ=255)
                             tacactacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTg                                                    >  1:1729922/1‑55 (MQ=255)
                             tacactacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCAtt                                    >  1:3699777/1‑71 (MQ=255)
                              acactacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGaa                                                 >  1:2032766/1‑57 (MQ=255)
                              acactacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGa                                                  >  1:473451/1‑56 (MQ=255)
                               cactacaACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCAttt                                   <  1:2844408/70‑1 (MQ=255)
                                actacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGt                                 >  1:2019675/1‑71 (MQ=255)
                                actacaACATTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGt                                 >  1:179161/1‑71 (MQ=255)
                                    caacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGa                                                  <  1:22111/50‑1 (MQ=255)
                                    caacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGa                                                  <  1:1326974/50‑1 (MQ=255)
                                    caacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGa                                                  <  1:126447/50‑1 (MQ=255)
                                    caacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGa                                                  <  1:3635540/50‑1 (MQ=255)
                                    caacaTTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGa                                                  <  1:269499/50‑1 (MQ=255)
                                    caacaTTGAGGCGAGGAAATGTTCATACCGCATgg                                                                 <  1:3606183/35‑1 (MQ=255)
                                    caacaTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGtt                                <  1:603277/68‑1 (MQ=255)
                                    caacaTTGAGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGtt                                <  1:891198/68‑1 (MQ=255)
                                       caTTGAGGCGAGGAAATGTTCATACCGCATGTATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCaa                          <  1:742220/71‑1 (MQ=255)
                                         ttGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCaaaa                        <  1:1168039/71‑1 (MQ=255)
                                         ttGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCaaaa                        <  1:1023576/71‑1 (MQ=255)
                                            aGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGc                     <  1:3710701/71‑1 (MQ=255)
                                            aGGCGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGc                     <  1:1482007/71‑1 (MQ=255)
                                               cGAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCaaaa                        <  1:1147108/65‑1 (MQ=255)
                                                gAGGAAATGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGcaaca                 <  1:344042/71‑1 (MQ=255)
                                                  ggAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGcaaca                 >  1:3299943/1‑69 (MQ=255)
                                                  ggAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGcaaca                 >  1:3464841/1‑69 (MQ=255)
                                                     aaTGTTCATACCGTATGGATTGTGGTATCTGGAAACg                                              >  1:1199585/1‑37 (MQ=255)
                                                     aaTGTTCATACCGTATGGATTGTGGTATCTGGAAACg                                              >  1:2652066/1‑37 (MQ=255)
                                                      aTGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACTAAGCCa           <  1:5717/71‑1 (MQ=255)
                                                      aTGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGGCa           <  1:393112/71‑1 (MQ=255)
                                                      aTGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCa           <  1:3488925/71‑1 (MQ=255)
                                                      aTGTTCATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCa           <  1:3316232/71‑1 (MQ=255)
                                                         ttCATACCGTATGGATTGTGGTATCTGGAAACGTCCt                                          >  1:3459918/1‑37 (MQ=255)
                                                         ttCATACCGTATGGATTGTGGTATCTGGAAACGTCCt                                          >  1:751434/1‑37 (MQ=255)
                                                         ttCATACCGTATGGATTGTGGTATCTGGAAACGTCCt                                          >  1:2257632/1‑37 (MQ=255)
                                                         ttCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGt         <  1:344783/70‑1 (MQ=255)
                                                         ttCATACCGCATGGATTGTGGTATCTGGAAACGTCCTc                                         >  1:64051/1‑38 (MQ=255)
                                                         ttCATACCGCATGGATTGTGGTATCTGGAAACGTCCTc                                         >  1:2470072/1‑38 (MQ=255)
                                                         ttCATACCGCATGGATTGTGGTATCTGGAAACGTCCTc                                         >  1:1820967/1‑38 (MQ=255)
                                                         ttCATACCGCATGGATTGTGGTATCTGGAAACGTCCTc                                         >  1:2062744/1‑38 (MQ=255)
                                                         ttCATACCGCATGGATTGTGGTATCTGGAAACGTCCTc                                         >  1:2099467/1‑38 (MQ=255)
                                                          tCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGcc            >  1:3104464/1‑66 (MQ=255)
                                                          tCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGcc            >  1:3503075/1‑66 (MQ=255)
                                                           cATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACaaa               >  1:1566077/1‑62 (MQ=255)
                                                           cATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACaaa               >  1:1440833/1‑62 (MQ=255)
                                                           cATACCGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACaaa               >  1:869006/1‑62 (MQ=255)
                                                            aTACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGaa      <  1:2146995/70‑1 (MQ=255)
                                                               ccGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATc   <  1:3442792/70‑1 (MQ=255)
                                                               ccGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCTACAAAGCCAGTGAAATCa  >  1:3415784/1‑71 (MQ=255)
                                                               ccGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCa  >  1:2134293/1‑71 (MQ=255)
                                                               ccGCATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCa  <  1:161776/71‑1 (MQ=255)
                                                                  |                                                                   
ACGTCTGTATTAATAGAAACGCCGCATTCTACACTACAACATTGAGGCGAGGAAATGTTCATACCGTATGGATTGTGGTATCTGGAAACGTCCTCGCATTTGTTATGCAAAATGCAACAAAGCCAGTGAAATCA  >  Exported/1739877‑1740010

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: