Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,758 A→C 19.8% intergenic (+1526/‑253) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7580AC19.8% 202.3 / 64.2 121intergenic (+1526/‑253)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/54);  new base C (14/10);  total (57/64)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-02

Only 100 of 125 total aligned reads displayed.

ATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAA  >  minE/2572693‑2572827
                                                                 |                                                                     
aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATg                                                                    >  1:943676/1‑69 (MQ=255)
   aaaaTTCCACCTACCTTCCAGCCCAATCGGAGCATGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTAcc                                                               >  1:4425874/1‑71 (MQ=255)
       ttCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGat                                                            >  1:3171790/1‑70 (MQ=255)
         ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa                                                          >  1:4211728/1‑70 (MQ=255)
         ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa                                                          <  1:5065269/70‑1 (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)
                cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGa                                                             <  1:3796863/60‑1 (MQ=255)
                cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                                                    <  1:4763221/69‑1 (MQ=255)
                cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                                                    <  1:5310379/69‑1 (MQ=255)
                 ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                                                 <  1:5513053/71‑1 (MQ=255)
                 ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                                                 >  1:1943673/1‑71 (MQ=255)
                   tCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                                                 >  1:1695969/1‑69 (MQ=255)
                    ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTTCTGCGGTTAATGTTACCGATAATAGTTATAGt                                                <  1:5113808/69‑1 (MQ=255)
                     cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                                               >  1:1811317/1‑69 (MQ=255)
                     cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa                                               >  1:4094487/1‑69 (MQ=255)
                     cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTAc                                              <  1:2124118/70‑1 (MQ=255)
                        ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                          >  1:2273757/1‑71 (MQ=255)
                        cccAATCGGGTCCGGGGTGATATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGcc                                          >  1:64164/1‑71 (MQ=255)
                         cccaTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGccc                                         >  1:1898110/4‑71 (MQ=255)
                         cccaTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGccc                                         >  1:1181916/4‑71 (MQ=255)
                         ccAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGccc                                         >  1:2761371/1‑71 (MQ=255)
                          cAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGata                                                           >  1:3425088/1‑52 (MQ=255)
                          cAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                        <  1:28268/71‑1 (MQ=255)
                           caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg                                       <  1:2543541/70‑1 (MQ=255)
                           caTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTg                                       <  1:3753246/70‑1 (MQ=255)
                           aaTCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTaa                                                                      <  1:4299080/40‑1 (MQ=255)
                            aTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGa                                      >  1:3363095/1‑71 (MQ=255)
                            aTCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCACTGa                                      >  1:3547067/1‑71 (MQ=255)
                            aTCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                                                                  >  1:2463427/1‑43 (MQ=255)
                             tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt                                     >  1:4835131/1‑71 (MQ=255)
                             tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt                                     >  1:2551937/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)
                                  tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                <  1:2827269/70‑1 (MQ=255)
                                  tccTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATtcctc                                <  1:5675327/70‑1 (MQ=255)
                                   ccTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGt                                                <  1:5209863/54‑1 (MQ=255)
                                   ccTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatca                               >  1:1458927/1‑71 (MQ=255)
                                    cTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                              <  1:1215386/71‑1 (MQ=255)
                                     tGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGt                                                <  1:1550200/52‑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:4311450/1‑61 (MQ=255)
                                        gtgtTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt                                    >  1:1850727/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)
                                          aatataCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGCCCTGATTCCTCATCAtgct                        >  1:3625482/3‑71 (MQ=255)
                                           ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                        <  1:459614/54‑1 (MQ=255)
                                           ttATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCt                                        <  1:4608315/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)
                                             ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTtat                                                   >  1:4843679/1‑41 (MQ=255)
                                             ataCCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATagtact                     >  1:311533/1‑71 (MQ=255)
                                              taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                          >  1:3381715/1‑49 (MQ=255)
                                              taCCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc                                          >  1:3598182/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)
                                                   ccccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCg               >  1:1675088/1‑71 (MQ=255)
                                                     ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGa                                      >  1:4672802/1‑46 (MQ=255)
                                                     ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGa                                      >  1:707945/1‑46 (MQ=255)
                                                     ccTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGa                                      >  1:4879658/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:878015/55‑1 (MQ=255)
                                                      ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                            <  1:1131489/55‑1 (MQ=255)
                                                      ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                            <  1:2079703/55‑1 (MQ=255)
                                                      ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                            <  1:3793395/55‑1 (MQ=255)
                                                      ctgctgCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCa            <  1:3059692/71‑1 (MQ=255)
                                                       tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                            <  1:1512071/54‑1 (MQ=255)
                                                       tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                            <  1:303116/54‑1 (MQ=255)
                                                       tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                            <  1:4509804/54‑1 (MQ=255)
                                                       tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                            <  1:3481780/54‑1 (MQ=255)
                                                       tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcatca                            <  1:105533/54‑1 (MQ=255)
                                                       tgctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAg           >  1:3132992/1‑71 (MQ=255)
                                                        gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATtcatcat                              <  1:946158/51‑1 (MQ=255)
                                                        gctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGc             >  1:3048195/1‑68 (MQ=255)
                                                         ctgcGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTc         >  1:67307/1‑71 (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:2148182/69‑1 (MQ=255)
                                                               ttAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCtgat     <  1:3604419/69‑1 (MQ=255)
                                                                tAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGaa  >  1:587680/1‑71 (MQ=255)
                                                                 aaTGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATactact                     >  1:5012677/1‑51 (MQ=255)
                                                                 |                                                                     
ATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTCATCATCATACTACTGCTGCGGCAGTCTGATGAA  >  minE/2572693‑2572827

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: