Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,572,770 A→G 19.7% intergenic (+1538/‑241) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,572,7700AG19.7% 242.1 / 53.2 122intergenic (+1538/‑241)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/52);  new base G (14/10);  total (60/62)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-02

Only 100 of 131 total aligned reads displayed.

CCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTCTG  >  Exported/2572702‑2572839
                                                                    |                                                                     
ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa                                                                      >  1:4211728/1‑70 (MQ=255)
 cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                                                                     <  1:2640986/70‑1 (MQ=255)
   cctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt                                                                  <  1:2116581/71‑1 (MQ=255)
      acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                                                                >  1:2439551/1‑70 (MQ=255)
       cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTTCGGTTAATGTTACCGATAATAGTta                                                                <  1:3063471/69‑1 (MQ=255)
       cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                                                                <  1:5310379/69‑1 (MQ=255)
        ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                                                             <  1:5513053/71‑1 (MQ=255)
          tCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                                                             >  1:1695969/1‑69 (MQ=255)
           ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTTCTGCGGTTAATGTTACCGATAATAGTTATAGt                                                            <  1:5113808/69‑1 (MQ=255)
           ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                                                           <  1:5045414/70‑1 (MQ=255)
            cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                                                           >  1:2144679/1‑69 (MQ=255)
            cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTAc                                                          <  1:2124118/70‑1 (MQ=255)
               ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                      >  1:2273757/1‑71 (MQ=255)
               ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                      <  1:1764035/71‑1 (MQ=255)
               cccAATCGGGTCCGGGGTGATATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGcc                                                      >  1:64164/1‑71 (MQ=255)
                cccaTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGccc                                                     >  1:1898110/4‑71 (MQ=255)
                ccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGccc                                                     >  1:2761371/1‑71 (MQ=255)
                 cAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                                    <  1:28268/71‑1 (MQ=255)
                  caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg                                                   <  1:2543541/70‑1 (MQ=255)
                  caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg                                                   <  1:3753246/70‑1 (MQ=255)
                   aTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGa                                                  >  1:3363095/1‑71 (MQ=255)
                    tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt                                                 >  1:2551937/1‑71 (MQ=255)
                    tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt                                                 >  1:4835131/1‑71 (MQ=255)
                      ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata                                                                    >  1:4172927/1‑50 (MQ=255)
                      ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt                                                >  1:5074099/1‑70 (MQ=255)
                      ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtc                                               <  1:3677966/71‑1 (MQ=255)
                          ccTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatca                                           >  1:1458927/1‑71 (MQ=255)
                           cTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          <  1:1215386/71‑1 (MQ=255)
                            tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGt                                                            <  1:1550200/52‑1 (MQ=255)
                            tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTAc                                                          >  1:3606836/1‑54 (MQ=255)
                            tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          >  1:4255053/1‑70 (MQ=255)
                            tGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctcat                                          >  1:4848809/1‑70 (MQ=255)
                             gggtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACtcctcatca                                        <  1:3320751/69‑1 (MQ=255)
                             gAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAg                                                                   <  1:5317765/44‑1 (MQ=255)
                              ggtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACTcctcatcat                                       <  1:3500579/70‑1 (MQ=255)
                              aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:4024081/70‑1 (MQ=255)
                              aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        >  1:4364238/1‑70 (MQ=255)
                              aGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTcatcatcat                                       <  1:9140/71‑1 (MQ=255)
                               gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt                                                >  1:1850727/1‑61 (MQ=255)
                               gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt                                                >  1:3075723/1‑61 (MQ=255)
                               gtgaTATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACtcctcat                                          <  1:4888930/67‑1 (MQ=255)
                                taatataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctcat                                          <  1:1756528/63‑1 (MQ=255)
                                 gtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          <  1:3750460/65‑1 (MQ=255)
                                 gtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          <  1:3995348/65‑1 (MQ=255)
                                  ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                                    <  1:4608315/54‑1 (MQ=255)
                                  ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                                    <  1:5120786/54‑1 (MQ=255)
                                  ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTcatcatcat                                       <  1:1729571/67‑1 (MQ=255)
                                  ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                    <  1:704595/70‑1 (MQ=255)
                                  ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                    <  1:1963142/70‑1 (MQ=255)
                                  atataCCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGcccc                                                    >  1:1408611/2‑53 (MQ=255)
                                  atataCCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGCCCCGACTCCTCATCAtgctg                                   <  1:1260426/70‑2 (MQ=255)
                                  atataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                            >  1:2377207/2‑62 (MQ=255)
                                  atataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCAtgctg                                   >  1:5647503/2‑70 (MQ=255)
                                    ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTtata                                                              >  1:1599657/1‑42 (MQ=255)
                                    ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTtata                                                              >  1:3483156/1‑42 (MQ=255)
                                    ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTtat                                                               >  1:4843679/1‑41 (MQ=255)
                                    ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATagtact                                 >  1:311533/1‑71 (MQ=255)
                                     taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                      >  1:3598182/1‑49 (MQ=255)
                                     taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                                      >  1:3381715/1‑49 (MQ=255)
                                     taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCATGctgctg                                >  1:1084867/1‑71 (MQ=255)
                                      aCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          <  1:3216827/60‑1 (MQ=255)
                                       ccGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTActgct                              <  1:3679169/71‑1 (MQ=255)
                                            ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGa                                                  >  1:4672802/1‑46 (MQ=255)
                                            ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGa                                                  >  1:707945/1‑46 (MQ=255)
                                            ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                    <  1:4925133/60‑1 (MQ=255)
                                            ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATactactact                              <  1:2410494/66‑1 (MQ=255)
                                            ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCgg                          >  1:1158939/1‑70 (MQ=255)
                                            ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                         <  1:527883/71‑1 (MQ=255)
                                             ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:2079703/55‑1 (MQ=255)
                                             ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:878015/55‑1 (MQ=255)
                                             ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:1131489/55‑1 (MQ=255)
                                             ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:3793395/55‑1 (MQ=255)
                                              tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:4503456/54‑1 (MQ=255)
                                              tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:303116/54‑1 (MQ=255)
                                              tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:4509804/54‑1 (MQ=255)
                                              tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:105533/54‑1 (MQ=255)
                                              tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                                        <  1:3481780/54‑1 (MQ=255)
                                               gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          <  1:946158/51‑1 (MQ=255)
                                               gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                                          <  1:1960161/51‑1 (MQ=255)
                                               gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                         >  1:3048195/1‑68 (MQ=255)
                                                ctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGt                      <  1:4689544/70‑1 (MQ=255)
                                                  gcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc                         >  1:5064551/1‑65 (MQ=255)
                                                   cGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtg                   <  1:4764864/70‑1 (MQ=255)
                                                     gTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgatg                <  1:3982140/71‑1 (MQ=255)
                                                      ttAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgat                 <  1:3604419/69‑1 (MQ=255)
                                                      ttAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgat                 <  1:5478436/69‑1 (MQ=255)
                                                       tAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa              >  1:610744/1‑71 (MQ=255)
                                                       tAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa              >  1:587680/1‑71 (MQ=255)
                                                        aaTGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATactact                                 >  1:2881159/1‑51 (MQ=255)
                                                        aaTGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa              >  1:1630218/1‑70 (MQ=255)
                                                          tGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtga                  >  1:1852637/1‑64 (MQ=255)
                                                             tACCGATAATAGTTATAGTACAGCCCTTATTCATCATCAta                                      >  1:5017773/1‑41 (MQ=255)
                                                             tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                      >  1:143625/1‑41 (MQ=255)
                                                             tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                      >  1:379287/1‑41 (MQ=255)
                                                             tACCGATAATAGTTATAGTACAGCCCTGATTCATCATCAta                                      >  1:2378101/1‑41 (MQ=255)
                                                               ccGATAATAGTTATAGTACAGCCCTGATTCATCATCAtact                                    <  1:4970106/41‑1 (MQ=255)
                                                               ccGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTa           >  1:1186802/1‑66 (MQ=255)
                                                               ccGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCt      >  1:678455/1‑71 (MQ=255)
                                                                   taataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttta  >  1:5567465/1‑68 (MQ=255)
                                                                   taataGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCtttta  <  1:3806569/71‑4 (MQ=255)
                                                                    |                                                                     
CCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAACTATTGCTTCTG  >  Exported/2572702‑2572839

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: