Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,572,727 G→T 43.9% intergenic (+1495/‑284) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,572,7270GT43.9% 131.1 / 160.9 157intergenic (+1495/‑284)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (45/43);  new base T (44/25);  total (89/68)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

GATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTC  >  minE/2572659‑2572794
                                                                    |                                                                   
gATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGcc                                                                   >  1:1086747/1‑71 (MQ=255)
    gATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGa                                                                <  1:1099265/70‑1 (MQ=255)
        ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                            <  1:1882051/70‑1 (MQ=255)
        ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                            <  1:4341378/70‑1 (MQ=255)
         cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                          >  1:28989/1‑71 (MQ=255)
         cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                          <  1:2017498/71‑1 (MQ=255)
           gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat                                                         <  1:3519381/67‑1 (MQ=255)
           gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                                        <  1:1731037/71‑1 (MQ=255)
            aaaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                                 >  1:2153154/1‑61 (MQ=255)
              aTCTGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt                                                              <  1:877913/62‑1 (MQ=255)
              aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt                                                              <  1:4897259/62‑1 (MQ=255)
               ccGGATTCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGc                                                    >  1:1231848/2‑71 (MQ=255)
               ccGGATTCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGc                                                    >  1:1660448/2‑71 (MQ=255)
                  ggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGcccc                                                 <  1:2673676/71‑1 (MQ=255)
                   attCAGCGGATACGCATACAGACTCCATCCCCCCTCCAGCCCTATCGGGTCCTGAGTGATATACCGCCCct                                                <  1:2205840/70‑1 (MQ=255)
                   attCAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGccc                                                  >  1:3231085/2‑69 (MQ=255)
                    ttCAACGTATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                          >  1:2577749/1‑60 (MQ=255)
                    ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAgt                                                              <  1:5182467/56‑1 (MQ=255)
                    ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                               >  1:1136220/1‑69 (MQ=255)
                    ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                          <  1:4411233/60‑1 (MQ=255)
                    ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                          >  1:2413257/1‑60 (MQ=255)
                     tCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCct                                                >  1:4019555/1‑69 (MQ=255)
                      cAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGCCCctgc                                              <  1:5249456/70‑1 (MQ=255)
                      cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                               >  1:3173499/1‑67 (MQ=255)
                      cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                             >  1:3994771/1‑71 (MQ=255)
                      aaaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                            >  1:711242/2‑56 (MQ=255)
                       aaCGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGCCCctgcaa                                            >  1:5397850/1‑69 (MQ=255)
                       aaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                                            >  1:3513883/1‑71 (MQ=255)
                        aCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                                            >  1:554642/1‑68 (MQ=255)
                        aCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                                            >  1:3815513/1‑68 (MQ=255)
                        aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                       >  1:610003/1‑59 (MQ=255)
                         cGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGCCCCTGCAGCg                                          <  1:4503678/71‑1 (MQ=255)
                         cGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCg                                          <  1:3995756/71‑1 (MQ=255)
                         cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                            <  1:1588314/53‑1 (MQ=255)
                          ggATACTGATAAAAATTCCATCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCTCCCCTGCAGCggg                                        >  1:2878568/1‑70 (MQ=255)
                          ggATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga                                                            <  1:5222836/52‑2 (MQ=255)
                          ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                       >  1:1156018/1‑57 (MQ=255)
                          ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                       >  1:5633958/1‑57 (MQ=255)
                           gATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGccc                                                  <  1:175360/61‑1 (MQ=255)
                           gATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtgat                                                           >  1:1654182/1‑50 (MQ=255)
                            aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                        >  1:2729850/1‑70 (MQ=255)
                             tACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAg                                                               >  1:2754358/1‑46 (MQ=255)
                             tACTGATAAAAATTCCATCCCCACTTCAGCCCAATCGGATCCTGAg                                                               >  1:4576623/1‑46 (MQ=255)
                              aCTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGa                                                                >  1:240013/1‑44 (MQ=255)
                               cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                                  >  1:4294827/1‑57 (MQ=255)
                                tGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                                  <  1:3630668/40‑1 (MQ=255)
                                tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCt                                                                  <  1:2534791/40‑1 (MQ=255)
                                tGATAAAAATTCCACCTACCTCCCAGCCCAATCGGAGCCt                                                                  <  1:432037/40‑1 (MQ=255)
                                 gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                                         >  1:1934447/1‑64 (MQ=255)
                                 gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                                         >  1:4979797/1‑64 (MQ=255)
                                 gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                                  <  1:5332975/39‑1 (MQ=16)
                                  aTAAAAATTCCATCCCGCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                                       >  1:3379299/1‑48 (MQ=255)
                                  aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                                       >  1:2503737/1‑48 (MQ=255)
                                  aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtata                                                        >  1:3729263/1‑48 (MQ=255)
                                  aTAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGa                                                                >  1:4320881/1‑40 (MQ=255)
                                  aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                                 <  1:1445588/71‑1 (MQ=255)
                                   tAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCagcagc                                           <  1:1807372/60‑1 (MQ=255)
                                    aaaaaTTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTAc                               >  1:1167399/1‑71 (MQ=255)
                                     aaaaTTCCACCTACCTTCCAGCCCAATCGGAGCATGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTAcc                              >  1:4425874/1‑71 (MQ=255)
                                       aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                                          >  1:1699164/1‑41 (MQ=255)
                                       aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAcctc                                                    <  1:4132582/47‑3 (MQ=255)
                                           ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa                         >  1:4211728/1‑70 (MQ=255)
                                            cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                        <  1:967637/70‑1 (MQ=255)
                                            cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                        <  1:2640986/70‑1 (MQ=255)
                                             acctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                           <  1:514564/50‑1 (MQ=255)
                                             acctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                           <  1:309814/50‑1 (MQ=255)
                                               ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                           <  1:399991/48‑1 (MQ=255)
                                               ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCgg                                         >  1:2143794/1‑50 (MQ=255)
                                               ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCg                                          >  1:2150441/1‑49 (MQ=255)
                                                tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                              <  1:5145287/44‑1 (MQ=255)
                                                tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                              <  1:3548074/44‑1 (MQ=255)
                                                tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                              <  1:158107/44‑1 (MQ=255)
                                                tacctTCCAGCCCAATCGGAGCCTGAGAGTTATACCGCCCCAGCTGCGGTTAAAGTTACCGAAAATAGAta                   >  1:2432791/1‑71 (MQ=255)
                                                  cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTTCGGTTAATGTTACCGATAATAGTta                   <  1:3063471/69‑1 (MQ=255)
                                                  cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGa                            <  1:3796863/60‑1 (MQ=255)
                                                  cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta                   <  1:4763221/69‑1 (MQ=255)
                                                   ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                <  1:5513053/71‑1 (MQ=255)
                                                   ctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAg                >  1:1943673/1‑71 (MQ=255)
                                                    ttCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTa                                >  1:815782/1‑54 (MQ=255)
                                                      tcAGCCCAATCGGATCCTGAGTGATATACCGcccc                                                 <  1:1132602/34‑1 (MQ=37)
                                                      tcAGCCCAATCGGATCCTGAGTGATATACCGcccc                                                 <  1:4242852/34‑1 (MQ=37)
                                                      tcAGCCCAATCGGATCCTGAGTGATATACCGcccc                                                 <  1:3332834/34‑1 (MQ=37)
                                                      ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                              >  1:274979/1‑38 (MQ=255)
                                                      ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTTCTGCGGTTAATGTTACCGATAATAGTTATAGt               <  1:5113808/69‑1 (MQ=255)
                                                      ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa              <  1:5045414/70‑1 (MQ=255)
                                                       cAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCgg                                         >  1:701797/1‑42 (MQ=255)
                                                       cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa              >  1:1811317/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)
                                                              aTCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGtt                                 >  1:2463427/1‑43 (MQ=255)
                                                               tCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAt    >  1:4835131/1‑71 (MQ=255)
                                                                 ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataata                       >  1:4172927/1‑50 (MQ=255)
                                                                 ggAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGAtt   >  1:5074099/1‑70 (MQ=255)
                                                                    |                                                                   
GATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCCCTGATTC  >  minE/2572659‑2572794

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: