Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 863,793 G→T 8.0% E32* (GAA→TAA)  yciQ → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported863,7930GT 8.0% 252.5 / 5.7 75E32* (GAA→TAA) yciQpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (28/41);  new base T (3/3);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

TTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAA  >  Exported/863724‑863861
                                                                     |                                                                    
ttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                      >  1:1769590/1‑70 (MQ=255)
ttGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                      <  1:4531384/70‑1 (MQ=255)
 tGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                      >  1:215396/1‑69 (MQ=255)
  gCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                      >  1:1225798/1‑68 (MQ=255)
   cTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATgaag                                                                   <  1:285349/70‑1 (MQ=255)
      ttGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAgg                                                               >  1:2065495/1‑71 (MQ=255)
      ttGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAgg                                                               >  1:1696408/1‑71 (MQ=255)
          tAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtc                                                           >  1:1590703/1‑71 (MQ=255)
            gCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtctc                                                         >  1:2046183/1‑71 (MQ=255)
             cAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtctc                                                         >  1:220302/1‑70 (MQ=255)
             cAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCg                                                        >  1:4117445/1‑71 (MQ=255)
              aGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCg                                                        <  1:2007772/70‑1 (MQ=255)
              aGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCg                                                        <  1:1889471/70‑1 (MQ=255)
                     tttGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGaagaag                                                                <  1:4417079/55‑1 (MQ=255)
                     tttGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGaagaag                                                                <  1:3138625/55‑1 (MQ=255)
                     tttGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTg                                                 <  1:1280187/70‑1 (MQ=255)
                       tGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATg                                              >  1:3393488/1‑71 (MQ=255)
                         tATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGtt                                            <  1:306970/71‑1 (MQ=255)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAgg                                                               <  1:1781082/50‑1 (MQ=25)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAgg                                                               <  1:3494380/50‑1 (MQ=25)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTg                                           <  1:4247866/70‑1 (MQ=255)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTg                                           <  1:574056/70‑1 (MQ=255)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:3635176/71‑1 (MQ=255)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:1603042/71‑1 (MQ=255)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:3781643/71‑1 (MQ=255)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:4410732/71‑1 (MQ=255)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:699198/71‑1 (MQ=255)
                           tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:1019086/71‑1 (MQ=255)
                             atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGa                                                             >  1:270836/1‑50 (MQ=21)
                             atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGa                                                             >  1:3087692/1‑50 (MQ=21)
                             atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGa                                                             >  1:829399/1‑50 (MQ=21)
                             atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAg                                        >  1:4423694/1‑71 (MQ=255)
                             atcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAg                                        >  1:1666873/1‑71 (MQ=255)
                                  tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCa                                     <  1:4261111/69‑1 (MQ=255)
                                  tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCa                                     <  1:2587104/69‑1 (MQ=255)
                                  tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCa                                     <  1:2632304/69‑1 (MQ=255)
                                  tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCa                                     <  1:1884001/69‑1 (MQ=255)
                                  tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAct                                   >  1:2234809/1‑71 (MQ=255)
                                    aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGa                                                >  1:2750890/1‑56 (MQ=255)
                                    aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGa                                                >  1:2995220/1‑56 (MQ=255)
                                    aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCActct                                 <  1:4695120/71‑1 (MQ=255)
                                    aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAct                                   <  1:640971/69‑1 (MQ=255)
                                    aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCAActct                                 >  1:2307132/1‑71 (MQ=255)
                                     gTTATGCAGAAAACACGGAGATCCCTTCTTATgaaga                                                                  <  1:4407412/37‑1 (MQ=255)
                                      ttATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTg                                >  1:2616481/1‑70 (MQ=255)
                                       tATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGc                               >  1:3133113/1‑70 (MQ=255)
                                       tATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGc                               >  1:1385211/1‑70 (MQ=255)
                                        aTGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCaa                             >  1:4431598/1‑71 (MQ=255)
                                        aTGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCGTTGATGTTGAAGCCACTCTGCaa                             >  1:1155152/1‑71 (MQ=255)
                                           cAGAAAACACGGAGATCCCTTCTTATGTAGAAGGGATCTCGCTCTTTGATGTTGa                                          >  1:2020278/1‑55 (MQ=255)
                                           cAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:1108434/71‑1 (MQ=255)
                                            aGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAg                         <  1:3976209/71‑1 (MQ=255)
                                             gAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAg                         <  1:4383889/70‑1 (MQ=255)
                                                aaCACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                     >  1:2178086/1‑71 (MQ=255)
                                                  cacGGAGATCCCTTCTTATTAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:1272243/48‑1 (MQ=255)
                                                  cacGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                          <  1:1640396/48‑1 (MQ=255)
                                                      gagaTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAct                                   <  1:1833316/51‑1 (MQ=255)
                                                         aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:1289934/57‑1 (MQ=255)
                                                         aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:3627075/57‑1 (MQ=255)
                                                         aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:1506010/57‑1 (MQ=255)
                                                          tCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGAc           <  1:369722/71‑1 (MQ=255)
                                                          tCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGAc           <  1:2320173/71‑1 (MQ=255)
                                                           cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:2210202/55‑1 (MQ=255)
                                                           cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:2176745/55‑1 (MQ=255)
                                                           cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:28689/55‑1 (MQ=255)
                                                            ccTTCTTATGAAGAAGTTATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                          <  1:262598/54‑1 (MQ=255)
                                                            ccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAg                         <  1:3163966/55‑1 (MQ=255)
                                                              ttcttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATggg                    <  1:1653596/58‑1 (MQ=255)
                                                              ttcttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATc        >  1:1456205/1‑70 (MQ=255)
                                                                cttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCa       >  1:3547894/1‑69 (MQ=255)
                                                                 ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            >  1:4398449/1‑63 (MQ=255)
                                                                 ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            <  1:3824416/63‑1 (MQ=255)
                                                                 ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            <  1:1652286/63‑1 (MQ=255)
                                                                 ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            <  1:3327521/63‑1 (MQ=255)
                                                                 ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa            >  1:4738369/1‑63 (MQ=255)
                                                                   aTGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATg                      >  1:1463483/1‑51 (MQ=255)
                                                                   aTGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATg                      >  1:4130459/1‑51 (MQ=255)
                                                                   aTGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAg    >  1:1100416/1‑69 (MQ=255)
                                                                   aTGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGa   >  1:3336499/1‑70 (MQ=255)
                                                                    tGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGaa  <  1:518666/70‑1 (MQ=255)
                                                                    tGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGaa  <  1:3871669/70‑1 (MQ=255)
                                                                     |                                                                    
TTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAA  >  Exported/863724‑863861

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: