Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,557,809 G→T 5.9% A155S (GCC→TCC)  secB → protein export chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,557,8090GT 5.9% 388.7 / 4.6 113A155S (GCC→TCC) secBprotein export chaperone
Reads supporting (aligned to +/- strand):  ref base G (50/56);  new base T (5/2);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

TAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTGCCA  >  Exported/2557739‑2557877
                                                                      |                                                                    
tAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATg                                                                      >  1:1936576/1‑71 (MQ=255)
tAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATg                                                                      >  1:127235/1‑71 (MQ=255)
 aaCTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGc                                                                     >  1:1963380/1‑71 (MQ=255)
    ttCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCt                                                                   >  1:1791758/1‑70 (MQ=255)
     tCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtga                                                                 >  1:1010840/1‑71 (MQ=255)
      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                >  1:629047/1‑71 (MQ=255)
      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                >  1:1631171/1‑71 (MQ=255)
      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                >  1:3119868/1‑71 (MQ=255)
      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                >  1:1500639/1‑71 (MQ=255)
      cGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtg                                                                  <  1:1247075/69‑1 (MQ=255)
       gATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGCATCCtgat                                                                <  1:3737989/70‑1 (MQ=255)
       gATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCtgat                                                                <  1:3464701/70‑1 (MQ=255)
         tGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGaa                                                             >  1:2181313/1‑71 (MQ=255)
           cgcTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAAcc                                                           <  1:1043723/71‑1 (MQ=255)
            gcTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCa                                                          <  1:3713169/71‑1 (MQ=255)
            gcTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCa                                                          <  1:3401574/71‑1 (MQ=255)
            gcTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCa                                                          <  1:3269516/71‑1 (MQ=255)
             cTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCaa                                                         <  1:1368688/71‑1 (MQ=255)
             cTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCaa                                                         <  1:480809/71‑1 (MQ=255)
             cTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCa                                                          <  1:1105298/70‑1 (MQ=255)
                ttCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACg                                                       >  1:1123821/1‑70 (MQ=255)
                   aTGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTaa                                                    <  1:3021589/70‑1 (MQ=255)
                    tGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGaa                                                             >  1:1423354/1‑60 (MQ=255)
                     gAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGc                                                 >  1:3674928/1‑71 (MQ=255)
                        cTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGt                                                      <  1:2575221/63‑1 (MQ=255)
                        cTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGt                                                      <  1:2579604/63‑1 (MQ=255)
                        cTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGt                                                      <  1:746977/63‑1 (MQ=255)
                        cTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTc                                              <  1:916661/71‑1 (MQ=255)
                          aTTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACTTAATGCTTCaa                                            >  1:644038/1‑71 (MQ=255)
                          aTTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTaa                                                    >  1:624539/1‑63 (MQ=255)
                          aTTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTaa                                                    >  1:794151/1‑63 (MQ=255)
                             tGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTAATGAc                                        >  1:674477/1‑71 (MQ=255)
                               cagcagcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGAc                                        <  1:2143029/70‑1 (MQ=255)
                                  cagcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgtg                                    <  1:143893/71‑1 (MQ=255)
                                  cagcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgt                                     >  1:3489215/1‑70 (MQ=255)
                                  cagcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgt                                     >  1:1571835/1‑70 (MQ=255)
                                    gcagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGa                                   >  1:872594/1‑70 (MQ=255)
                                     cagGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTTATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:1744750/71‑1 (MQ=255)
                                      agGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCg                                >  1:2181527/1‑71 (MQ=255)
                                      agGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCg                                >  1:1217168/1‑71 (MQ=255)
                                      agGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCg                                >  1:1559478/1‑71 (MQ=255)
                                      agGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCg                                >  1:507549/1‑71 (MQ=255)
                                        gCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGt                              <  1:1545981/71‑1 (MQ=255)
                                          tGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtg                                      >  1:3169690/1‑61 (MQ=255)
                                          tGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGc                            >  1:1976305/1‑71 (MQ=255)
                                            gCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGc                            <  1:1276473/69‑1 (MQ=255)
                                            gCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCg                          >  1:2141495/1‑71 (MQ=255)
                                              gAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAAAGACTGTGATCGGTGCCgg                         >  1:573549/1‑70 (MQ=255)
                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTc                                              <  1:288534/47‑1 (MQ=255)
                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtg                                      >  1:3634398/1‑55 (MQ=255)
                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtg                                      >  1:3629965/1‑55 (MQ=255)
                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:1569690/1‑71 (MQ=255)
                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:2989645/1‑71 (MQ=255)
                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:2989161/1‑71 (MQ=255)
                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:2829414/1‑71 (MQ=255)
                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:2720194/1‑71 (MQ=255)
                                                aGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      >  1:827047/1‑71 (MQ=255)
                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:2711177/71‑1 (MQ=255)
                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:3406379/71‑1 (MQ=255)
                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:1180309/71‑1 (MQ=255)
                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:3156958/71‑1 (MQ=255)
                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     >  1:714420/1‑71 (MQ=255)
                                                 ggTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCg                     <  1:2941791/71‑1 (MQ=255)
                                                  gTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgt                                     <  1:2430258/54‑1 (MQ=255)
                                                  gTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACtgt                                     <  1:1241235/54‑1 (MQ=255)
                                                  gTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 >  1:3048263/1‑58 (MQ=255)
                                                    aCTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGa                                   <  1:2458824/54‑1 (MQ=255)
                                                    aCTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGa                                   >  1:1879032/1‑54 (MQ=255)
                                                     cTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTAc                  <  1:1771923/70‑1 (MQ=255)
                                                      tGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCa                                             <  1:2996987/42‑1 (MQ=255)
                                                      tGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCa                                             <  1:1668620/42‑1 (MQ=255)
                                                       gaagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCa                                             <  1:2274385/41‑1 (MQ=255)
                                                       gaagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCa                                             <  1:1745229/41‑1 (MQ=255)
                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTg                             >  1:1040093/1‑56 (MQ=255)
                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTg                             >  1:515175/1‑56 (MQ=255)
                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACgg                <  1:1398186/69‑1 (MQ=255)
                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACgg                <  1:2495991/69‑1 (MQ=255)
                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACgg                <  1:855288/69‑1 (MQ=255)
                                                        aagaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCa              >  1:1288020/1‑71 (MQ=255)
                                                         agaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCAc             >  1:1846527/1‑71 (MQ=255)
                                                         agaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCAc             >  1:2494856/1‑71 (MQ=255)
                                                          gaaCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTAc                  >  1:3344252/1‑65 (MQ=255)
                                                            aCATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGc          >  1:2592756/1‑71 (MQ=255)
                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGcc                           >  1:3057337/1‑53 (MQ=38)
                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGc                            >  1:2125016/1‑52 (MQ=38)
                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGc                            >  1:1371687/1‑52 (MQ=38)
                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGt                    >  1:1664242/1‑60 (MQ=255)
                                                             cATCAGGCATCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTa                   >  1:2123739/1‑61 (MQ=255)
                                                                cAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCtt      >  1:2032265/1‑71 (MQ=255)
                                                                 aGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTg     <  1:1352153/71‑1 (MQ=255)
                                                                   gATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      <  1:2909443/52‑1 (MQ=255)
                                                                   gATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTc                      <  1:2317709/52‑1 (MQ=255)
                                                                    aggCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:3663940/38‑1 (MQ=255)
                                                                    aTGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:518848/40‑1 (MQ=255)
                                                                    aTGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:2434030/40‑1 (MQ=255)
                                                                    aTGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:3153104/40‑1 (MQ=255)
                                                                    aTGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:3445742/40‑1 (MQ=255)
                                                                     tGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATc                                 <  1:3715663/39‑1 (MQ=255)
                                                                     tGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTGCCa  <  1:2263256/70‑1 (MQ=255)
                                                                      gCCTGATGAACCAACGTAATGCTTCAATGACtgtg                                    <  1:3215041/35‑1 (MQ=255)
                                                                      |                                                                    
TAACTTCGATGCGCTGTTCATGAACTATTTGCAGCAGCAGGCTGGCGAAGGTACTGAAGAACATCAGGATGCCTGATGAACCAACGTAATGCTTCAATGACTGTGATCGGTGCCGGCTCGTACGGCACCGCTCTTGCCA  >  Exported/2557739‑2557877

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